RU2819561C1 - Способ определения коэффициента диффузии в массивных изделиях из капиллярно-пористых материалов - Google Patents

Abstract

Изобретение относится к измерительной технике и может быть использовано при исследовании процессов массопереноса для определения коэффициента диффузии в строительных изделиях из капиллярно-пористых материалов, а также в пищевой, химической и других отраслях промышленности. Способ определения коэффициента диффузии в массивных изделиях из капиллярно-пористых материалов заключается в том, что в исследуемом изделии создают равномерное начальное содержание распределенного в твердой фазе растворителя, гидроизолируют верхнюю плоскую поверхность образца, в начальный момент времени осуществляют импульсное точечное увлажнение верхней поверхности исследуемого изделия. Измеряют изменение во времени сигналов менее удаленного первого и более удаленного второго от источника гальванических датчиков, расположенных на расстояниях соответственно r ₁ и r ₂ от точки нанесения импульса дозой растворителя. Фиксируют моменты времени τ₁ и τ₂, при которых достигаются одинаковые значения сигналов соответственно первого E ₁ и второго E ₂ датчиков из диапазона (0,7 – 0,9)E _e на нисходящих ветвях кривых изменения сигналов во времени этих двух датчиков и рассчитывают коэффициент диффузии. При этом вначале делают пробное импульсное точечное увлажнение верхней поверхности исследуемого изделия дозой растворителя. Если доза пробного импульсного точечного увлажнения оказывается недостаточной, при которой максимумы сигналов обоих датчиков после нанесения этого воздействия оказываются меньше (0,7)E _e или не наблюдается вовсе, то осуществляют новое точечное импульсное воздействие увеличенной в (r ₂/r ₁)³ раза дозой растворителя. Процедуру повторяют до момента идентификации выраженного максимума сигнала любого из датчиков в диапазоне 1>E _max/E _e≥0,7. Если доза пробного импульсного точечного увлажнения оказывается избыточной, при которой после нанесения этого воздействия максимумы сигналов обоих датчиков находятся в зоне плато насыщения датчиков E _max/E _e≈1 и не идентифицируются, то осуществляют новое импульсное точечное воздействие уменьшенной в (r ₂/r ₁)³ раза дозой растворителя. Процедуру повторяют до момента наблюдения выраженного максимума сигнала любого из датчиков в диапазоне 1>E _max/E _e≥0,7. Если относительное значение максимума сигнала E _max2/E _e более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в диапазоне 0,75>E _max2/E _e ≥0,7, то осуществляют новое импульсное воздействие, превышающее предыдущее приблизительно в 1,2 раза, а затем проводят расчет искомого коэффициента диффузии. Если относительное значение максимума сигнала E _max1/E _e менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 0,75>E _max1/E _e ≥0,7, то осуществляют новое импульсное воздействие, превышающее предыдущее приблизительно в 1,2×(r ₂/r ₁)³ раза. Затем проводят расчет искомого коэффициента диффузии. Если относительное значение максимума сигнала более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в пределах 1>E _max2/E _e≥0,95, то осуществляют новое импульсное воздействие уменьшенной в 2,25 раза дозой растворителя, затем производят расчет искомого коэффициента диффузии. Если относительное значение максимума сигнала E _max1/E _e менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 1>E _max1/E _e ≥0,95, то дозу следующего импульса увеличивают в (r ₂/r ₁)³/2,25 раза, после чего производят расчет искомого коэффициента диффузии. Если относительное значение максимума сигнала E _max2/E _e более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в диапазоне 0,95>E _max2/E _e ≥0,75, то производят расчет искомого коэффициента диффузии. Если относительное значение максимума сигнала менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 0,95>E _max1/E _e≥0,75, то осуществляют следующее воздействие увеличенной в (r ₂/r ₁)³ раза дозой, а затем производят расчет искомого коэффициента диффузии. Техническим результатом является повышение точности измерения коэффициента диффузии и сокращение длительности эксперимента. 3 з.п. ф-лы, 3 ил., 4 табл.

Description

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

Известен способ определения коэффициента массопроводности и потенциалопроводности массопереноса (А.С. 174005, кл. G 01 k N 421, 9₅₁, 1965), заключающийся в импульсном увлажнении слоя материала и измерении на заданном расстоянии от этого слоя изменения влагосодержания материала во времени. Коэффициент массопроводности вычисляется по установленной зависимости. Недостатком этого способа являются осуществление разрушающего контроля опытного образца при размещении датчиков во внутренних слоях исследуемого тела, невозможность определения коэффициента диффузии других растворителей, кроме воды, большая трудоемкость метода при подготовке образцов, необходимость индивидуальной градуировки датчиков по каждому материалу.

Наиболее близким является способ определения коэффициента диффузии растворителей в массивных изделиях из капиллярно-пористых материалов (патент РФ на изобретение № 2784198, G 01 N 13/00, G 01 N 27/26, 23.11.2022, Бюл. № 33), заключающийся в том, что в исследуемом изделии создают равномерное начальное содержание распределенного в твердой фазе растворителя, гидроизолируют верхнюю плоскую поверхность образца, в начальный момент времени осуществляют импульсное точечное увлажнение верхней поверхности исследуемого изделия, затем измеряют изменение во времени сигналов двух гальванических преобразователей, расположенных на разных расстояниях r ₁ и r ₂ от точки нанесения импульса дозой растворителя, фиксируют моменты времени τ₁ и τ₂, при которых достигаются одинаковые значения сигналов соответственно первого датчика E ₁ и второго датчика E ₂ из диапазона (0,7 – 0,9)E _e на нисходящих ветвях кривых изменения сигналов во времени этих двух датчиков и рассчитывают коэффициент диффузии, где E _e - максимально возможное значение сигнала датчиков, соответствующее переходу растворителя из области связанного с твердой фазой исследуемого материала в область свободного состояния.

Недостатками этого способа являются:

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

2. Низкая чувствительность работы применяемых гальванических преобразователей при завышенной дозе вносимого растворителя при импульсном воздействии по сравнению с требуемым (заранее неизвестным), что делает невозможным применение данного метода. При этом может возникать ситуация, когда в течение значительного времени эксперимента сигналы одного или обоих датчиков практически не изменяются и наблюдаются на уровне E ₁/E _e = E ₂/E _e ≈1 (фигура 2, кривая 1), т.к. находятся в зоне плато насыщения статической характеристики датчиков. Причем негативные последствия превышения вносимой дозы увеличиваются по мере отклонения в большую сторону величины вносимой дозы.

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

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

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

Техническая задача достигается тем, что в отличие от прототипа (патент РФ на изобретение № 2784198, G 01 N 13/00, G 01 N 27/26, 23.11.2022, Бюл. № 33) вначале делают пробное импульсное точечное увлажнение верхней поверхности исследуемого изделия дозой растворителя, затем, если доза пробного импульсного точечного увлажнения оказывается недостаточной, при которой максимумы сигналов обоих датчиков после нанесения этого воздействия оказываются меньше (0,7)E _e или не наблюдается вовсе, то осуществляют новое точечное импульсное воздействие увеличенной в (r ₂/r ₁)³ раза дозой растворителя, причем эту процедуру повторяют до момента идентификации выраженного максимума сигнала любого из датчиков в диапазоне 1˃E _max/E _e ≥0,7.

Если доза пробного импульсного точечного увлажнения оказывается избыточной, при которой после нанесения этого воздействия максимумы сигналов обоих датчиков находятся в зоне плато насыщения датчиков E _max/E _e≈1 и не идентифицируются, то осуществляют новое импульсное точечное воздействие уменьшенной в (r ₂/r ₁)³ раза дозой растворителя, причем эту процедуру повторяют до момента наблюдения выраженного максимума сигнала любого из датчиков в диапазоне 1˃E _max/E _e≥0,7.

После этого, если относительное значение максимума сигнала E _max2/E _e более удаленного от точки нанесения импульсного воздействия второго датчика после этих воздействий идентифицируется в диапазоне 0,75˃E _max2/E _e≥0,7, то осуществляют новое импульсное воздействие, превышающее предыдущее приблизительно в 1,2 раза, а затем проводят расчет искомого коэффициента диффузии.

Если относительное значение максимума сигнала E _max1/E _e менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 0,75˃E _max1/E _e ≥0,7, то осуществляют новое импульсное воздействие, превышающее предыдущее приблизительно в 1,2×(r ₂/r ₁)³, а затем проводят расчет искомого коэффициента диффузии. Если относительное значение максимума сигнала более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в пределах 1˃E _max2/ E _e ≥0,95, то осуществляют новое импульсное воздействие уменьшенной в 2,25 раза дозой растворителя и производят расчет искомого коэффициента диффузии.

Если относительное значение максимума сигнала E _max1/E _e менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 1˃E _max1/ E _e ≥0,95, то дозу следующего импульса увеличивают в (r ₂/r ₁)³/2,25 раза, после чего производят расчет искомого коэффициента диффузии.

Если относительное значение максимума сигнала E _max2/E _e более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в диапазоне 0,95˃E _max2/E _e ≥0,75, то производят расчет искомого коэффициента диффузии.

Если относительное значение максимума сигнала менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 0,95˃E _max1/E _e≥0,75, то осуществляют следующее воздействие увеличенной в (r ₂/r ₁)³ раза дозой, а затем производят расчет искомого коэффициента диффузии.

Причем после каждого нанесенного импульсного воздействия следующее импульсное воздействие производят после снижения сигналов преобразователей до начального значения, а расчет искомого коэффициента диффузии производят при относительном значении максимума сигнала второго датчика E _max2/E _e≥0,75 и при значениях сигналов обоих датчиков E ₁ и E ₂, приблизительно равных (E _max2 - 0,05E _e).

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

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

Применяемые гальванические преобразователи обладают монотонной зависимостью выходной характеристики - электродвижущей силы (ЭДС) от концентрации распределенного в твердой фазе растворителя. Предпочтительным для использования при определении коэффициента диффузии является рабочий участок статической характеристики гальванических преобразователей, который находится в диапазоне (0,7 – 0,9)E _e. При уменьшении содержания растворителя за пределы, соответствующего 0,7E _e, наблюдается нестабильная работа гальванических преобразователей с низкой помехозащищенностью сигнала или сигнал вообще отсутствует. При увеличении содержания растворителя за пределы, соответствующего 0,9E _e, наблюдается интенсивное снижение чувствительности ЭДС преобразователя к изменению содержания растворителя (плато насыщения). При появлении свободных молекул растворителя, не связанных с твердой фазой, сигнал преобразователя стабилизируется на уровне E _e. Поэтому равные значения сигналов обоих датчиков E ₁ и E ₂, при которых фиксируются моменты времени и , входящие в расчетное выражение для определения искомого коэффициента диффузии, необходимо иметь возможность выбирать из диапазона (0,7 – 0,9)E _e. Причем, концентрация растворителя, соответствующая сигналу датчика равного 0,9E _e, приблизительно в 2 раза выше концентрации растворителя, соответствующей сигналу датчика равного 0,7E _e. При сигналах датчика больших 0,9E _e за счет интенсивного снижения чувствительности ЭДС преобразователя к увеличению содержания растворителя это соотношение имеет тенденцию к существенному росту. Так, концентрация растворителя, соответствующая сигналу датчика равного 0,95E _e, выше концентрации растворителя, соответствующей сигналу датчика равного 0,75E _e, уже приблизительно в 2,25 раза (фигура 3, кривые 2, 4).

Процесс распространения растворителя в массивном изделии из капиллярно-пористых материалов (при условии, что минимальные размеры изделия относительно точки импульсного воздействия превышают 10 r ₂, где r ₂ - расстояние от точки импульсного воздействия до электродов наиболее удаленного от нее гальванического преобразователя) после подачи импульса в виде дозы растворителя описывается краевой задачей массопереноса в неограниченной среде при нанесении импульсного воздействия от точечного источника массы.

В этом случае изменение концентрации растворителя в капиллярно-пористом материале в зоне действия источника описывается функцией:

где - концентрация растворителя на поверхности сферы радиусом r относительно точки импульсного подвода дозы растворителя к образцу в момент времени τ; D - коэффициент диффузии растворителя; – плотность абсолютно сухого исследуемого материала; Q – количество жидкой фазы, подведенной из дозатора к плоской поверхности изделия исследуемого капиллярно-пористого материала.

Коэффициент диффузии растворителя D при организации данного процесса массопереноса в изделии связан соотношением:

где τ_{max i} – время, соответствующее максимуму на кривой U(r _i ,τ) изменения концентрации на расстоянии r _i от источника.

Учитывая (2), уравнение (1) для заданной точки контроля ЭДС гальванического преобразователя r=r _i можно преобразовать к виду:

Из (3) можно получить значение достигаемого максимума U _{max i} при τ=τ_{max i} :

Из уравнения (4) с учетом (2), получено:

Из (5) следует, что для обеспечения значения максимума E _max2 кривой на более удаленном втором датчике, равном максимуму E _max1 кривой на ближнем первом датчике, необходимо увеличивать дозу Q ₁ импульсного точечного увлажнения в (r ₂/r ₁)³ раза. И, наоборот, для обеспечения значения максимума E _max1 кривой на менее удаленном первом датчике, равном максимуму E _max2 кривой на более удаленном втором датчике, необходимо уменьшать дозу Q ₂ импульсного точечного увлажнения в (r ₂/r ₁)³ раза.

Расчетное соотношение для определения искомого коэффициента диффузии имеет вид:

где и - моменты времени, при которых фиксируются одинаковые значения сигналов соответственно первого датчика E ₁ и второго датчика E ₂ из диапазона (0,7 – 0,9)E _e на нисходящих ветвях кривых изменения сигналов во времени этих двух датчиков.

Среднеквадратическая оценка относительной погрешности измерения искомого коэффициента диффузии по расчетному соотношению (6), имеет вид:

где и – относительная погрешность определения моментов времени соответственно и (при условии равенства абсолютных погрешностей определения моментов времени);

В формулах (8) и (9) символами ∆ обозначены абсолютные погрешности определения разности и логарифма .

При фиксированных значениях r ₁ и r ₂, реализованных в устройстве, погрешности разности и логарифма , определяемые по формулам (8), (9), являются доминантами результирующей погрешности измерения искомого коэффициента диффузии (7), т.к. в них также присутствуют погрешности и .

При увеличении вносимой дозы растворителя уменьшается разница между значениями и , входящих в расчетное выражение (7) (фигура 3, кривые 1 и 2), что приводит к росту относительных погрешностей (8),(9), т.к. при уменьшении разницы между значениями и уменьшаются оба знаменателя в выражениях (8) и (9), а стремится к нулю. Поэтому измерение искомого коэффициента диффузии необходимо проводить не только в области стабильной работы применяемых преобразователей в диапазоне (0,7 – 0,9)E _e, но и при возможно большей разнице значений и .

Если доза пробного импульсного точечного увлажнения оказывается недостаточной, при которой максимумы сигналов обоих датчиков после нанесения этого воздействия оказываются меньше (0,7)E _e или не наблюдаются вовсе, то осуществляют новое точечное импульсное воздействие увеличенной в (r ₂/r ₁)³ раза дозой растворителя, причем эту процедуру повторяют до момента идентификации выраженного максимума сигнала любого из датчиков в диапазоне 1˃E _max/E _e ≥0,7 (фигура 1, кривые 1 - 4).

Если доза пробного импульсного точечного увлажнения оказывается избыточной, при которой после нанесения этого воздействия максимумы сигналов обоих датчиков находятся в зоне плато насыщения датчиков E _max/E _e≈1 и не идентифицируются, то осуществляют новое импульсное точечное воздействие уменьшенной в (r ₂/r ₁)³ раза дозой растворителя, причем эту процедуру повторяют до момента наблюдения выраженного максимума сигнала любого из датчиков в диапазоне 1˃E _max/E _e ≥0,7 (фигура 2, кривые 2 или 3 и 4).

Обеспечение получения выраженного максимума сигнала любого из датчиков в диапазоне 1˃E _max/E _e ≥0,7 позволяет находить рациональные значения доз импульсного точечного увлажнения уже на следующем шаге исследования.

Так, если в результате предыдущих итераций на более удаленном от источника втором датчике наблюдается значение максимума на нижней границе рабочего участка статической характеристики (0,7 – 0,9)E _e, т.е. E _max2≈ 0,7E _e (фигура 1, кривая 3), то определение искомого коэффициента диффузии связано с существенной погрешностью измерения момента времени τ₂ вследствие низкой чувствительности датчика вблизи максимума кривой, где производная сигнала по времени стремится к нулю. Как показала практика, в этом случае целесообразно получение максимума на более удаленном втором датчике на уровне не ниже E _max2 ≈ 0,75 E _e, а расчет искомого коэффициента диффузии проводить на основании данных о моментах времени и , при которых фиксируются значения сигналов обоих датчиков E ₁ и E ₂, равных (E _max2 - 0,05E _e). Это обеспечивается применением дополнительно увеличенной в 1,2 раза дозы, при которой достигалось значение E _max2 ≈ 0,7E _e на предыдущем шаге итераций. После этого рассчитывают искомый коэффициент диффузии по формуле с использованием кривых 1 и 2 (фигура 1) по формуле (6). В этом случае при значении максимума на более удаленном втором датчике даже на нижнем уровне E _max2 ≈ 0,75 E _e, значение момента времени надежно фиксируется на нижнем уровне стабильного и помехозащищенного сигнала гальванического преобразователя E ₂/E _e ≈ 0,7. Аналогичная картина наблюдается и при относительном значении E _max2/E _e, находящемся не на границе, а внутри диапазона 0,75˃E _max2/E _e˃0,7, с той лишь разницей, что получаемое в результате нового импульсного воздействия значение E _max2/E _e будет выше 0,75. Моменты времени и , входящих в расчетное выражение (6) рекомендуется фиксировать при значениях сигналов обоих датчиков E ₁ и E ₂, равных (E _max2 - 0,05E _e).

Если в результате предыдущих итераций наблюдается минимально допустимое значение E _max1 ≈ 0,7E _e на ближнем первом датчике, (фигура 1, кривая 4), а стабильного сигнала на более удаленном втором датчике вообще не наблюдается (на фигуре 1 нет соответствующей кривой), то при увеличении дозы в (r ₂/r ₁)³ раза, получается то же значение максимума E _max2 ≈ 0,7E _e, но на более удаленном втором датчике (фигура 1, кривая 3). Это значит, что для достижения максимума на более удаленном втором датчике на уровне не ниже E _max2 ≈ 0,75 E _e, необходимо дополнительно увеличить наносимую в процессе импульсного воздействия дозу в 1,2 раза (рассмотрено выше). Таким образом, если в результате предыдущих итераций наблюдается минимально допустимое значение E _max1 ≈ 0,7E _e на ближнем первом датчике, для повышения точности необходимо увеличить предыдущую дозу приблизительно в 1,2×(r ₂/r ₁)² раза и рассчитывать искомый коэффициент диффузии по формуле (6) на основании данных кривых 1 и 2 (фигура 1). Аналогичная картина наблюдается и при относительном значении E _max1/E _e, находящемся не на границе, а внутри диапазона 0,75˃E _max1/E _e˃0,7, с той лишь разницей, что получаемое в результате нового импульсного воздействия значение E _max2/E _e будет выше 0,75. Моменты времени и , входящих в расчетное выражение (6) рекомендуется фиксировать при значениях сигналов обоих датчиков E ₁ и E ₂, равных (E _max2 - 0,05E _e).

Рассмотрим теперь ситуации, когда в результате предыдущих итераций наблюдается значение максимума сигналов E _max1/E _e менее удаленного от точки нанесения импульсного воздействия первого датчика или E _max2/E _e более удаленного второго датчика находятся в диапазонах соответственно 1˃E _max1/E _e ≥0,95 и 1˃E _max2/E _e ≥0,95.

Если относительное значение максимума сигнала более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в диапазоне 1˃E _max2/E _e≥0,95, то осуществляют новое импульсное воздействие уменьшенной в 2,25 раза дозой растворителя, после чего производят расчет искомого коэффициента диффузии по формуле (6) на основании данных кривых 3 и 4 (фигура 3). Если при этом относительное значение максимума сигнала более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется на нижней границе диапазона 1˃E _max2/E _e≥0,95, т.е. E _max2 ≈ 0,95E _e (фигура 3, кривая 2), то после уменьшения дозы в 2,25 раза E _max2 достигает уровня E _max2 ≈ 0,75E _e. Это случай уже рассмотренный выше. Если E _max2 надежно идентифицируется внутри диапазона 1˃E _max2/E _e˃0,95, то после уменьшения дозы в 2,25 раза E _max2 будет находиться в рабочем диапазоне 0,95˃E _max2/E _e˃0,75. Верхнее ограничение 0,95˃E _max2/E _e подтверждается кривой 4 (фигура 2) для случая завышенного значения дозы, при которой максимум сигнала менее удаленного от точки нанесения импульсного воздействия первого датчика не идентифицируется и равен E _max1/E _e ≈1. Необходимость применения на последнем шаге итерации уменьшенной дозы в 2,25 раза по сравнению с предыдущей итерацией объясняется стремлением повысить точность определения искомого коэффициента диффузии. В самом деле, расчет искомого коэффициента диффузии может быть произведен на предыдущем шаге итерации по кривым 1, 2 (фигура 3). Однако, при увеличенной дозе растворителя на предыдущем шаге итерации существенно увеличивается длительность эксперимента, а также уменьшается разница между значениями и , входящих в расчетное выражение (6) (фигура 3, кривые 1 и 2), что приводит к росту относительных погрешностей (8),(9), т.к. при уменьшении разницы между значениями и уменьшаются оба знаменателя в выражениях (8) и (9), а стремится к нулю. Поэтому измерение искомого коэффициента диффузии необходимо проводить не только в области стабильной работы применяемых преобразователей в диапазоне (0,7 – 0,9)E _e, но и при возможно большей разнице значений и .

Если в результате предыдущих итераций наблюдается значение максимума сигнала E _max1/E _e менее удаленного от точки нанесения импульсного воздействия первого датчика в диапазоне 1˃E _max1/E _e ≥0,95, то дозу следующего импульса осуществляют с коэффициентом (r ₂/r ₁)³/2,25. после чего производят расчет искомого коэффициента диффузии по формуле (6) на основании данных кривых 3 и 4 (фигура 3). При этом, если увеличить предыдущую дозу в (r ₂/r ₁)³ раза, то на более удаленном втором датчике будет наблюдаться максимум сигнала в диапазоне 1˃E _max2/E _e ≥0,95. Если еще уменьшить полученную дозу в 2,25 раза, то это приведет к случаю, рассмотренном в предыдущем абзаце. Таким образом, дозу в целом нужно увеличивать в (r ₂/r ₁)³/2,25 раза.

Если в результате предыдущих итераций относительное значение максимума сигнала более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в диапазоне 0,95˃E _max2/E _e≥0,75, то производят расчет искомого коэффициента диффузии по формуле (6) на основании данных кривых 3 и 4 (фигура 2).

Если в результате предыдущих итераций относительное значение максимума сигнала менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 0,95˃E _max1/E _e≥0,75, то осуществляют воздействие увеличенной в (r ₂/r ₁)³ раза дозой, чтобы получить относительное значение максимума сигнала более удаленного от точки нанесения импульсного воздействия второго датчика в диапазоне 0,95˃E _max2/E _e≥0,75, после чего производят расчет искомого коэффициента диффузии по формуле (6).

Примеры. Были проведены исследования коэффициента диффузии этанола в плитах, отформованных из пеногипсобетона, толщиной 50 мм, плотностью в сухом состоянии 600 кг/м³. Расстояние от точки нанесения дозы растворителя до расположения электродов гальванических преобразователей: r ₁ = 4 мм и r ₂ = 6 мм. Количество внесенного растворителя определялось по мерной емкости. Исследования проводились при комнатной температуре.

Пример 1. На фигуре 1 представлены графики изменения ЭДС гальванических преобразователей при различных дозах импульсного точечного увлажнения при поиске рациональной дозы для расчета искомого коэффициента диффузии.

В результате предыдущей итерации при дозе ≈3.3×10^-6 кг наблюдается минимально допустимое значение E _max1 ≈ 0,7E _e на ближнем первом датчике, (фигура 1, кривая 4), а стабильного сигнала на более удаленном втором датчике вообще не наблюдается (на фигуре 1 нет соответствующей кривой), поэтому произвести расчет искомого коэффициента диффузии не представляется возможным. После увеличения в (6/4)³ раза наносимой при импульсном воздействии дозы до уровня ≈1.1×10^-5 кг наблюдался максимум сигнала более удаленного второго датчика на уровне E _max2/E _e ≈ 0,7 (фигура 1, кривая 3). Однако определение искомого коэффициента диффузии в этом случае связано с существенной погрешностью измерения момента времени τ₂ вследствие низкой чувствительности датчика вблизи максимума кривой, где производная сигнала по времени стремится к нулю. При дополнительном увеличении дозы приблизительно в 1,2 раза до уровня порядка 1,34 ×10^-5 кг наблюдается значение максимума более удаленного второго датчика на уровне E _max2/E _e ≈0,76. Таким образом, суммарное увеличение дозы произошло в 1,2×(r ₂/r ₁)³ раза.

Если предыдущую итерацию сразу проводили при дозе ≈1.1×10^-5 кг, то максимум сигнала более удаленного второго датчика наблюдался на уровне E _max2/E _e ≈0,7 (фигура 1, кривая 3). Поэтому следующую дозу импульсного воздействия необходимо увеличивать приблизительно в 1,2 раза до уровня порядка 1,34 ×10^-5 кг, что приводит к значению максимума более удаленного второго датчика на уровне E _max2/E _e ≈0,76.

И в первом, и во втором случаях расчет искомого коэффициента диффузии осуществлялся при значениях сигналов обоих датчиков E ₁ и E ₂, приблизительно равных (E _max2/E _e - 0,05)≈0,71 на основании данных кривых 1 и 2 (фигура 1). Результаты расчета по формуле (6) представлены в таблице 1.

Пример 2. В результате предыдущей итерации при дозе ≈1.68×10^-4 кг максимумы сигналов обоих датчиков находились в зоне плато насыщения применяемых датчиков и не идентифицировались, в течение практически всего эксперимента наблюдались сигналы E _max1/E _e=E _max2/E _e≈1 (фигура 2, кривые не представлены).

Таблица 1.

Значение ЭДС E 1/E e = E 2/E e	,с	,с	,с		D×109, м2/с.
0,71	4195	2717	1478	0,652	3,62

При уменьшении дозы в (6/4)³ раза до ≈4.98×10^-5 кг, идентифицировался максимум сигнала более удаленного второго датчика на уровне E _max2/E _e≈0,995, на ближнем первом датчике максимум сигнала не идентифицировался, т.к. находился в зоне плато насыщения применяемых датчиков E _max1/E _e ≈1 (фигура 2, кривые 1 и 2).

После уменьшения в (6/4)³ раза наносимой при импульсном воздействии дозы до уровня ≈1.47×10^-5 кг максимум сигнала более удаленного второго датчика наблюдался на уровне E _max2/E _e≈0,79 (фигура 2, кривая 4).

Расчет коэффициента диффузии можно проводить при различных равных значениях сигналов обоих датчиков E ₁ и E ₂. Однако определение искомого коэффициента диффузии при E ₁=E ₂=0,79E _e связано с существенной погрешностью измерения момента времени τ₂ вследствие низкой чувствительности датчика вблизи максимума кривой, где производная сигнала по времени стремится к нулю. Поэтому целесообразно расчеты проводить на участках кривой E ₂/E _e(τ) за пределами области расположения максимума при равных значениях ЭДС преобразователей E ₁ и E ₂ из всего рационального диапазона (0,7 – 0,9)E _e. В таблице 2 представлены результаты измерения при различных значениях ЭДС преобразователей.

Анализ данных таблицы 2 показывает, что при выборе для расчета искомого коэффициента диффузии более низкого уровня приравниваемых значений E ₁ и E ₂ закономерно снижаются значения и особенно , что приводит к росту составляющих (8) и (9) результирующей погрешности измерения искомого коэффициента диффузии (7).

Таблица 2.

Значение ЭДС E 1/E e = E 2/E e	,с	,с	,с		D×109, м2/с.
0,74	4219	2729	1490	0,653	3,60
0,7	4644	3156	1488	0,579	3,44

Например, фиксирование моментов времени и при E ₁=E ₂=0.7E _e увеличивает погрешность (9) по сравнению с фиксированием при E ₁=E ₂=0.74E _e за счет уменьшения в знаменателе (9). Это происходит вследствие приближения кривых друг к другу при одновременном увеличении длительности эксперимента. Для повышения точности измерения искомого коэффициента диффузии целесообразно фиксировать значения моментов времени и при максимально возможных одинаковых значениях E ₁ и E ₂. С целью снижения негативного влияния повышения погрешности в окрестности максимума кривой 4 фиксирование моментов времени τ₁ и τ₂ целесообразно проводить при одинаковых значениях E ₁ и E ₂, меньших максимума E _max2 приблизительно на 0.05E _e:

Кроме того, при данном значении сокращается длительность эксперимента по сравнению с использованием значения E ₁=E ₂=0.7E _e.

Пример 3. В результате предыдущей итерации при дозе ≈9.74×10^-5 кг максимумы сигналов обоих датчиков не идентифицировались, т.к. находились в зоне насыщения применяемых датчиков, и в течение практически всего эксперимента наблюдались сигналы E ₁/E _e=E ₂/E _e≈1 (фигура 3, кривые не представлены). Произвели уменьшение дозы в (6/4)³ раза до уровня ≈2.88×10^-5 кг. После этого наблюдалось изменение сигналов датчиков в виде кривых 1 и 2 (фигура 3) с максимумом кривой 2 на уровне E _max2/E _e ≈0,95. На основании данных этого эксперимента производился расчет искомого коэффициента диффузии, результаты которого представлены в таблице 3. Использование значений E ₁ и E ₂ вблизи максимума кривой 2 приводит к увеличению погрешности измерения момента времени за счет снижения чувствительности изменения ЭДС более удаленного от источника преобразователя от времени, где производная сигнала по времени стремится к нулю (фигура 3, кривая 2). Поэтому с целью снижения негативного влияния повышения погрешности в окрестности максимума кривой 2 максимальные равные значения E ₁ и E ₂, при которых производят фиксирование моментов времени τ₁ и τ₂, входящих в формулу (6), выбирали меньше максимума E _max2 приблизительно на 0.05E _e.

Таблица 3.

Значение ЭДС E 1/E e = E 2/E e	,с	,с	,с		D×109, м2/с.
0,9	4601	3235	1366	0,528	3,62
0,85	5381	4173	1208	0,381	3,71
0,8	6159	4960	1199	0,325	3,59

Затем дополнительно снизили дозу в 2,25 раза до уровня ≈1.28×10^-5 кг и получили изменение ЭДС на двух датчиках в виде кривых 3 и 4 (фигура 3) с максимумом E _max2/E _e≈0,75. Расчет искомого коэффициента диффузии осуществлялся при значениях сигналов обоих датчиков E ₁ и E ₂, приблизительно равных (E _max2/E _e - 0,05)≈0,7. Результаты расчета по формуле (6) представлены в таблице 4.

Таблица 4.

Значение ЭДС E 1/E e = E 2/E e	,с	,с	,с		D×109, м2/с.
0,7	4195	2671	1524	0,677	3,57

Анализ данных, приведенных в таблицах 3, 4 показывает, что с увеличением дозы снижаются значения и особенно , входящих в знаменатели составляющих (8) и (9), что приводит к росту составляющих (8) и (9) результирующей погрешности измерения искомого коэффициента диффузии (7).

Применение дополнительно уменьшенной в 2,25 раза дозы позволяет существенно сократить длительность эксперимента (таблица 4) и снизить погрешность составляющей (9) общей погрешности определения искомого коэффициента диффузии (7) приблизительно в (1,28 – 2,08) раза по сравнению с данными таблицы 3 за счет увеличения в знаменателе (9). Также обеспечивается снижение погрешности составляющей (8) общей погрешности определения искомого коэффициента диффузии (7) приблизительно в (1,12 – 1,27) раза за счет увеличения в знаменателе (8).

Claims (5)

1. Способ определения коэффициента диффузии в массивных изделиях из капиллярно-пористых материалов, заключающийся в том, что в исследуемом изделии создают равномерное начальное содержание распределенного в твердой фазе растворителя, гидроизолируют верхнюю плоскую поверхность образца, в начальный момент времени осуществляют импульсное точечное увлажнение верхней поверхности исследуемого изделия, затем измеряют изменение во времени сигналов менее удаленного первого и более удаленного второго от источника гальванических датчиков, расположенных на расстояниях соответственно r₁ и r₂ от точки нанесения импульса дозой растворителя, фиксируют моменты времени τ₁ и τ₂, при которых достигаются одинаковые значения сигналов соответственно первого E₁ и второго E₂ датчиков из диапазона (0,7 – 0,9)E_e на нисходящих ветвях кривых изменения сигналов во времени этих двух датчиков и рассчитывают коэффициент диффузии, отличающийся тем, что вначале делают пробное импульсное точечное увлажнение верхней поверхности исследуемого изделия дозой растворителя, затем, если доза пробного импульсного точечного увлажнения оказывается недостаточной, при которой максимумы сигналов обоих датчиков после нанесения этого воздействия оказываются меньше (0,7)E_e или не наблюдается вовсе, то осуществляют новое точечное импульсное воздействие увеличенной в (r₂/r₁)³ раза дозой растворителя, причем эту процедуру повторяют до момента идентификации выраженного максимума сигнала любого из датчиков в диапазоне 1˃E_max/E_e≥0,7, если доза пробного импульсного точечного увлажнения оказывается избыточной, при которой после нанесения этого воздействия максимумы сигналов обоих датчиков находятся в зоне плато насыщения датчиков E_max/E_e≈1 и не идентифицируются, то осуществляют новое импульсное точечное воздействие уменьшенной в (r₂/r₁)³ раза дозой растворителя, причем эту процедуру повторяют до момента наблюдения выраженного максимума сигнала любого из датчиков в диапазоне 1˃E_max/E_e≥0,7, после этого, если относительное значение максимума сигнала E_max2/E_e более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в диапазоне 0,75˃E_max2/E_e ≥0,7, то осуществляют новое импульсное воздействие, превышающее предыдущее приблизительно в 1,2 раза, а затем проводят расчет искомого коэффициента диффузии, если относительное значение максимума сигнала E_max1/E_e менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 0,75˃E_max1/E_e ≥0,7, то осуществляют новое импульсное воздействие, превышающее предыдущее приблизительно в 1,2×(r₂/r₁)³ раза, а затем проводят расчет искомого коэффициента диффузии, если относительное значение максимума сигнала более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в пределах 1˃E_max2/E_e≥0,95, то осуществляют новое импульсное воздействие уменьшенной в 2,25 раза дозой растворителя, а затем производят расчет искомого коэффициента диффузии, если относительное значение максимума сигнала E_max1/E_e менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 1˃E_max1/E_e ≥0,95, то дозу следующего импульса увеличивают в (r₂/r₁)³/2,25 раза, после чего производят расчет искомого коэффициента диффузии, если относительное значение максимума сигнала E_max2/E_e более удаленного от точки нанесения импульсного воздействия второго датчика идентифицируется в диапазоне 0,95˃E_max2/E_e ≥0,75, то производят расчет искомого коэффициента диффузии, если относительное значение максимума сигнала менее удаленного от точки нанесения импульсного воздействия первого датчика идентифицируется в диапазоне 0,95˃E_max1/E_e≥0,75, то осуществляют следующее воздействие увеличенной в (r₂/r₁)³ раза дозой, а затем производят расчет искомого коэффициента диффузии,

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

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

3. Способ по п. 1, отличающийся тем, что расчет искомого коэффициента диффузии производят при относительном значении максимума сигнала второго датчика E_max2/E_e≥0,75.

4. Способ по п. 1, отличающийся тем, что расчет искомого коэффициента диффузии производят при значениях сигналов обоих датчиков E₁ и E₂, приблизительно равных (E_max2/E_e - 0,05).