RU2322662C2 - Способ измерения коэффициента температуропроводности (варианты) и устройство для его осуществления - Google Patents
Способ измерения коэффициента температуропроводности (варианты) и устройство для его осуществления Download PDFInfo
- Publication number
- RU2322662C2 RU2322662C2 RU2006116343/28A RU2006116343A RU2322662C2 RU 2322662 C2 RU2322662 C2 RU 2322662C2 RU 2006116343/28 A RU2006116343/28 A RU 2006116343/28A RU 2006116343 A RU2006116343 A RU 2006116343A RU 2322662 C2 RU2322662 C2 RU 2322662C2
- Authority
- RU
- Russia
- Prior art keywords
- temperature
- heat exchanger
- measuring
- meter
- measured
- Prior art date
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Изобретение относится к строительной технике. Сущность изобретения заключается в формировании одномерного теплового потока, измерении зависимости от времени температуры теплоносителя на входе теплообменника и измерении зависимости от времени температуры между теплообменником и тепловой изоляцией. Техническим результатом изобретения является повышение точности и достоверности измерений. 3 н. и 3 з.п. ф-лы, 1 табл, 5 ил.
Description
Изобретение относится к строительной технике и может быть преимущественно использовано для измерения теплофизических характеристик различных строительных конструкций, например стен, потолков, полов, переборок, подволоков и др.
Известен способ бурения скважин и устройство для его осуществления [1], позволяющие получать образцы материалов с различной глубины. Измеряя параметры этих образцов, можно получить информацию о физических и химических свойствах и конфигурации глубинных слоев. Недостаток известных способа и устройства для его осуществления заключается в том, что они не обеспечивают неразрушающего контроля исследуемого объекта.
Известны многочисленные варианты способов ультразвуковой дефектоскопии и устройств, их реализующих, например [2, 3], позволяющие определить наличие неоднородностей в различных конструкциях и конфигурацию этих неоднородностей, однако приборы такого рода не позволяют провести измерение теплофизических характеристик исследуемых материалов, в частности теплового сопротивления.
Известны многочисленные варианты способов для измерения теплофизических характеристик различных радиоэлектронных приборов, например, описанный в [4] способ определения теплового сопротивления переход-корпус полупроводниковых диодов. Недостаток известного технического решения заключается в узкой области применения: его можно использовать для измерения теплофизических характеристик только радио-электронных приборов, причем только одного их класса - полупроводниковых диодов.
Известны многочисленные варианты устройств для измерения теплофизических характеристик различных радиоэлектронных приборов, например устройство для измерения теплового сопротивления транзисторов [5]. Недостаток известного технического решения заключается в узкой области применения: его можно использовать для измерения теплового сопротивления только радиоэлектронных приборов, причем только одного их класса - транзисторов.
Известно описанное в [6] устройство для определения характеристик материалов, содержащее источник импульсного нагрева, термопару и электронный блок обработки. Термопара расположена на поверхности исследуемого образца. Выход термопары подключен к входу электронного блока обработки. Главный недостаток известного устройства заключается в том, что при использовании импульсного нагрева необходима сложная обработка результатов измерений, для чего требуется сложная аппаратура. Это приводит к значительному удорожанию проведения измерений. Кроме того, большая сложность обработки результатов измерений приводит к снижению их точности и достоверности.
Наиболее близким по технической сущности к заявленному способу является способ бесконтактного неразрушающего контроля теплофизических свойств материалов [7]. Известный способ заключается в том, что измеряют двумя термоприемниками температуру в заданных точках поверхности образца и температуру окружающей среды, по полученным результатам определяют поправочный коэффициент, затем воздействуют на поверхность образца неподвижным точечным источником тепла. В заданный момент времени измеряют двумя термоприемниками избыточные температуры нагреваемой поверхности в заданных точках. Продолжают нагрев и измеряют момент времени, когда температура более удаленного от пятна нагрева термоприемника увеличится на заданную величину. По измеренным величинам определяют коэффициенты температуропроводности и теплопроводности.
Наиболее близким по технической сущности к заявленному устройству является термозонд для неразрушающего контроля теплофизических свойств материалов [8], содержащий линейный нагреватель и две термопары, расположенные симметрично относительно линейного нагревателя по обе стороны от него. Главный недостаток известного устройства заключается в том, что линейный нагреватель формирует тепловой поток, расходящийся в плоскости, перпендикулярной оси линейного нагревателя.
Известные технические решения (способ и устройство) имеют низкие потребительские свойства за счет узкой области применения, низкой точности и низкой достоверности измерений. Наличие указанных недостатков обусловлено следующими факторами. В известных технических решениях используется точечный или линейный источник тепловой энергии, поэтому его можно применять только при измерении теплофизических характеристик однородных объектов. Если исследуемый объект имеет различные неоднородности (например, железобетонная стена), то результаты измерений будут различаться при воздействии точечного или линейного источника тепловой энергии на различные точки поверхности исследуемого объекта. Низкие точность и достоверность известных технических решений обусловлены высокой сложностью модели, описывающей их работу. Так, в известном способе, необходимо провести предварительные измерения для определения поправочного коэффициента. Работа известных способа и устройства описываются уравнением теплового баланса, в котором необходимо учитывать потери тепловой мощности в окружающую среду за счет конвективного и лучистого теплообмена. Кроме того, в уравнении теплового баланса известного способа необходимо учитывать потери тепловой мощности из-за частичного поглощения лазерного излучения окружающей средой и частичного отражения лазерного излучения поверхностью исследуемого объекта. Эти компоненты учитываются приближенно расчетным путем. При обработке результатов измерений в известном способе используется функция ошибок, вычисляемая путем ее разложения в ряд Тейлора. В описании известного способа указано, что вычисление этой функции в аналитическом виде весьма затруднительно. Там же указано, что только для материалов с коэффициентом температуропроводности а≥10-7 м2/сек можно ограничиться первым членом ряда Тейлора, и только тогда можно использовать рабочую формулу для обработки результатов измерений. Все эти факторы в совокупности приводят к низким потребительским свойствам известных способа и устройства за счет узкой области применения и низких точности и достоверности.
Задачей изобретения является повышение потребительских свойств путем расширения области применения и повышения точности и достоверности.
Решение поставленной задачи в соответствии с п.1 формулы изобретения обеспечивается тем, что в известный способ, заключающийся в нагревании поверхности исследуемого объекта, измерении зависимости от времени температуры первым измерителем температуры, измерении зависимости от времени температуры вторым измерителем температуры, внесены следующие усовершенствования: нагревают участок поверхности исследуемого объекта посредством использования теплоносителя, поступающего в теплообменник, внешняя поверхность которого, кроме поверхности, примыкающей к участку поверхности исследуемого объекта, снабжена тепловой изоляцией, при этом измеряют зависимость от времени температуры теплоносителя на входе теплообменника первым измерителем температуры, измеряют зависимость температуры от времени вторым измерителем температуры, размещенным между теплообменником и тепловой изоляцией, и определяют коэффициент температуропроводности - а из соотношения:
а=4Fo*τ{V0(T0-Tw)/[2(Tw-Tн)Fo*+Vwτ]}2,
где Fo* - характеристическое значение критериального числа Фурье, τ - текущее время, V0=CM/(cγF), С - удельная теплоемкость теплоносителя, М - массовый расход теплоносителя, с - удельная теплоемкость материала исследуемого объекта, γ - плотность материала исследуемого объекта, F - площадь контакта теплообменника и исследуемого объекта, Т0 - температура, измеряемая первым измерителем температуры, Tw - температура, измеряемая вторым измерителем температуры, Тн - температура окружающей среды, Vw - скорость изменения температуры, измеряемой вторым измерителем температуры.
В частном случае в соответствии с п.2 формулы изобретения измеряют зависимость от времени температуры теплоносителя на выходе теплообменника третьим измерителем температуры и обеспечивают условия равенства показаний второго и третьего измерителей температуры.
Решение поставленной задачи в соответствии с п.3 формулы изобретения обеспечивается тем, что в известный способ, заключающийся в нагревании поверхности исследуемого объекта, измерении зависимости от времени температуры первым измерителем температуры, измерении зависимости от времени температуры вторым измерителем температуры, внесены следующие усовершенствования: нагревают участок поверхности исследуемого объекта посредством использования теплоносителя, поступающего в теплообменник, внешняя поверхность которого, кроме поверхности, примыкающей к участку поверхности исследуемого объекта, снабжена тепловой изоляцией, при этом измеряют зависимость от времени температуры теплоносителя на входе теплообменника первым измерителем температуры, измеряют зависимость температуры от времени вторым измерителем температуры, размещенным между теплообменником и тепловой изоляцией, причем коэффициент температуропроводности а определяют из соотношения:
a=Naэ,
N={cоэ(T0-Tw)[2(Twэ-Tн)Fo*+Vwэτ]}2/{co(T0-Twэ)[2(Twэ-Tн)Fo*+Vwτ]}2,
где аэ - коэффициент температуропроводности эталонного объекта соэ - объемная теплоемкость материала эталонного объекта, Тwэ - температура, измеряемая вторым измерителем температуры при исследовании эталонного объекта, Vwэ - скорость изменения температуры, измеряемой вторым измерителем температуры при исследовании эталонного объекта, c0 - объемная теплоемкость материала исследуемого объекта, Fo* - характеристическое значение критериального числа Фурье, τ - текущее время, Т0 - температура, измеряемая первым измерителем температуры, Tw - температура, измеряемая вторым измерителем температуры, Тн - температура окружающей среды.
В частном случае в соответствии с п.4 формулы изобретения измеряют зависимость от времени температуры теплоносителя на выходе теплообменника третьим измерителем температуры и обеспечивают условия равенства показаний второго и третьего измерителей температуры.
Решение поставленной задачи в соответствии с п.5 формулы изобретения обеспечивается тем, что в известное устройство, содержащее нагреватель, первый измеритель температуры, второй измеритель температуры, внесены следующие усовершенствования: оно дополнительно содержит входной трубопровод, соединительный трубопровод, теплообменник, выходной трубопровод, при этом выход входного трубопровода соединен с входом нагревателя, выполненного с возможностью нагрева теплоносителя, а выход нагревателя соединен с входом соединительного трубопровода, выход которого соединен с входом теплообменника, выход теплообменника соединен с входом выходного трубопровода, внешняя поверхность теплообменника, кроме поверхности, примыкающей к участку поверхности исследуемого объекта, снабжена тепловой изоляцией, первый измеритель температуры размещен в соединительном трубопроводе, а второй измеритель температуры размещен между теплообменником и тепловой изоляцией.
В частном случае в соответствии с п.6 формулы изобретения устройство дополнительно содержит третий измеритель температуры, размещенный в выходном трубопроводе.
Покажем, что поставленная задача изобретения действительно решается в заявленном техническом решении.
Первичной определяемой по данным измеренных параметров обычно является температуропроводность, а теплопроводность вычисляется по известному коэффициенту температуропроводности с помощью соотношения [7]:
где λ - коэффициент теплопроводности материала исследуемого объекта, [Вт/(м·К)]; а - коэффициент температуропроводности материала исследуемого объекта [м2/с]; c0 - объемная теплоемкость материала исследуемого объекта [Дж/(м3·К)], с - удельная теплоемкость материала исследуемого объекта [Дж/(кг·К)]; γ - плотность материала исследуемого объекта [кг/м3].
Следует отметить, что за рубежом температуропроводность принята в качестве эталонной величины, а в России традиционно используется эталон теплопроводности.
Заметим, что в ряде случаев конечной величиной, подлежащей определению, является удельное тепловое сопротивление, определяемое по формуле:
где r - удельное тепловое сопротивление [м2К/Вт], обратное величине коэффициента теплопередачи αс [Вт/м2К]; L - толщина прогреваемого слоя исследуемого объекта [м].
На практике величина удельного теплового сопротивления г должна определяться экспериментально, в частности, по экспериментально определенным температуропроводности или теплопроводности.
Физическая модель исследуемого объекта может быть представлена в виде однородной неограниченной пластины толщиной L. В начальный момент времени одна из поверхностей пластины вводится в соприкосновение с источником тепловой мощности - теплообменником, имеющим ту же температуру, что и исследуемая пластина.
Введем следующие допущения и ограничения.
1. Протяженность источника тепловой мощности вдоль поверхности исследуемого объекта достаточно велика для того, чтобы можно было считать, что тепловой поток через стенку одномерный и однородный по сечению.
2. Теплофизические параметры материала исследуемого объекта одинаковы по всем направлениям и постоянны во времени.
3. Теплообменник теплоизолирован от окружающей среды, кроме части его поверхности, соприкасающейся с поверхностью исследуемого объекта; таким образом, внешней средой для теплообменника является контактирующий с ним исследуемый объект, вглубь которого (в поперечном направлении) распространяется тепловой поток, выделяемый в теплообменнике за счет передачи тепловой энергии от нагретого теплоносителя.
4. Теплообменник является изотермическим, вследствие чего можно ограничиться анализом его среднеобъемной температуры.
Среднеобъемная температура проточного теплообменника в рамках принятых допущений будет описываться уравнением вида:
где Cw - полная теплоемкость теплообменника [Дж/К]; Tw - температура теплообменника [К]; Тc - температура поверхности исследуемого объекта, примыкающей к рабочей поверхности теплообменника [К]; Tf - температура теплоносителя [К]; σ - тепловая проводимость между поверхностями теплообменника и исследуемого объекта [Вт/К]; αк - контактный коэффициент теплообмена на границе соприкосновения поверхностей теплообменника и исследуемого объекта [Вт/м2К]; F - площадь контакта теплообменника и исследуемого объекта [м2]; α - конвективный коэффициент теплообмена в теплообменнике [Вт/м2К]; S - площадь внутренней поверхности теплообменника [м2], τ - текущее время.
В дальнейшем рассматривается среднее по теплообменнику значение коэффициента теплоотдачи α.
В реальных случаях термической инерцией теплоносителя можно пренебречь, тогда распределение температур в теплоносителе вдоль оси z, совпадающей с усредненным направлением движения теплоносителя в теплообменнике, будет описываться уравнением:
где С - удельная теплоемкость теплоносителя [Дж/(кг·К)]; М - его массовый расход [кг/с]; l - его длина [м]; Т0 - температура горячего теплоносителя на входе.
Из (4), считая теплообменник изотермическим, а распределение α по оси z постоянным, можно получить распределение температуры теплоносителя по оси z:
где 
- средняя по оси z температура теплоносителя.
Из (3) с учетом (5) получим:
Рассматривая параметры, входящие в (6) и определяющие тепловой режим теплообменника, необходимо отметить сложности априорной оценки φ - важной характеристики, определяющей эффективность охлаждения. Оценка φ требует расчета коэффициента конвективного теплообмена по критериальным соотношениям.
Безразмерный параметр φ является весьма общей характеристикой, которая присутствует в математическом описании среднеобъемных температур теплообменников. При φ>3...4 задача прогнозирования теплового режима предельно упрощается, поскольку в этом случае Е≈1 и кроме существенного упрощения расчетных соотношений отпадает необходимость расчета коэффициента конвективного теплообмена.
Основную формулу, описывающую процесс нагрева или охлаждения, удобно записать после преобразования соотношения (6) при условии Е=1 в виде:
В формулу (7) входит темп нагрева m0, который складывается из собственного темпа нагрева mн и дополнительной компоненты, возрастающей с увеличением массового расхода теплоносителя.
Заметим, что при выполнении условия φ>3...4 из (5) следует, что температура теплоносителя на выходе теплообменника Тв равна температуре теплообменника Tw, то есть выполняется условие:
Отдаваемая теплообменником тепловая мощность при условии Е=1 может быть определена из соотношения:
Подставив в (9) выражение для Tw из (7), получим:
Из (10) следует, что в начальный момент времени и в конце процесса выхода на стационарный тепловой режим значение отдаваемой тепловой мощности определяется соотношением:
Как показывает анализ, для эффективной работы рассматриваемой измерительной системы необходимо обеспечение выполнения условия:
откуда
Однако в процессе измерений с учетом допущения 3 теплопринимающей средой для теплообменника является исследуемый объект, вследствие конечной (и обычно малой) теплопроводности его материала особенно быстро нагревается поверхность исследуемого объекта, контактирующая с рабочей поверхностью теплообменника.
Таким образом, величина Тc в соотношениях (7), (10)-(12) не что иное, как температура нагреваемой поверхности исследуемого объекта.
Вследствие возрастания величины Тc во времени отдаваемая теплообменником тепловая мощность падает в процессе измерений от начального значения
до конечного
Отсюда следует, что соотношения (7) и (10) необходимо дополнить зависимостью Тс(τ), для чего необходимо определить нестационарное температурное поле в стенке, описываемое следующей краевой задачей:
В краевой задаче введены следующие обозначения: x - координата вдоль оси, ортогональной плоскости контакта теплообменника и исследуемого объекта [м]; q - удельный тепловой поток [Вт/м2]; Р - полный тепловой поток, отдаваемый теплообменником и определяемый соотношением (10).
Граничное условие (17в) введено из тех соображений, что процесс измерений длится гораздо меньше времени, чем требуется для прогрева всей толщи стенки.
Преобразуем (10), с учетом обозначения σΣ в (12), к следующему виду:
Заметим, что
где σм - новое введенное обозначение параметра, имеющего размерность тепловой проводимости.
Подставив (19) в (18), получим:
После подстановки (20) в (17б) можно получить граничное условие в следующем виде:
Получим оценку нестационарной температуры поверхности исследуемого объекта, проинтегрировав обе части уравнения (17а) по x в пределах от 0 до L:
Подставив (22) с учетом (21) в (17а), получим:
Для упрощения уравнения (23) удобно перейти от температур к перегреву ϑ=Т0-Тс0:
при начальном условии, определяемом с учетом (17 г):
Умножим все члены уравнения (24) на L и перейдем к безразмерным параметрам
где Fo - безразмерное критериальное число Фурье.
Решение уравнения (26) имеет вид:
В результате интегрирования в (27) с учетом начального условия (25) после преобразований получим окончательное решение уравнения (26) в виде:
Из (28) можно определить нестационарную температуру нагреваемой поверхности исследуемого объекта:
Подставив (29) в (7) с учетом тождества Тc≡Тсо, получим выражение для температуры теплообменника Тw, которая, как следует из (5), при φ>4 совпадает с температурой теплоносителя на выходе теплообменника:
Из (10) с учетом (29) можно также определить отдаваемый теплообменником тепловой поток:
Начальная стадия нагрева по определению будет ограничиваться таким отрезком времени от начала процесса нагрева, в пределах которого выполняются следующие неравенства:
При весьма малых значениях аргумента после разложения экспоненты в ряд можно ограничиться линейным приближением, то есть:
поэтому функции f1 и f2 с высокой степенью точности можно представить в виде:
где fin - приближенное значение функции fi на начальном временном участке (i=1, 2); Сс - полная теплоемкость, нагреваемого слоя стенки в пределах площади F контакта с теплообменником [Дж/К].
При выполнении условий (34) выражения (29) и (30) можно привести к виду:
Из (35) можно определить скорость роста температуры поверхности исследуемого объекта Vc:
Полученные соотношения для температур Тсо и Tw, а также для скорости роста температуры поверхности исследуемого объекта непригодны для определения теплофизических свойств (теплопроводности, температуропроводности) материала исследуемого объекта, поскольку рассмотренная модель физического процесса относится только к поверхности, оставляя без рассмотрения процесс распространения температурного поля вглубь исследуемого объекта. Поэтому необходимо получить решение для нестационарного температурного поля в прилегающем к поверхности слое материала исследуемого объекта.
Введем еще одно допущение (в дополнение к четырем ранее сформулированным).
Допущение 5:
что позволяет перейти к рассмотрению средней скорости роста температуры в нагреваемом слое.
Подставим Vc из (38) в левую часть уравнения (17а), полагая с учетом (37), что Vc=const. В результате получим новое уравнение:
Решение уравнения (39) имеет вид:
где С1 и С2 - постоянные интегрирования.
Из граничного условия (17б) получим:
Постоянную интегрирования С2 найдем из дополнительного граничного условия:
С учетом граничных условий (41) и (42) решение (40) примет вид:
Из (43) нетрудно определить температуру на поверхности Тсо=Тс(х=0); с учетом очевидного соотношения λ=асγ, получим:
Из (44) нетрудно определить значение коэффициента температуропроводности:
В полученное соотношение (45) входит в качестве параметра толщина прогреваемого слоя L, значение которой следует признать неопределенным, поскольку с течением времени эта толщина непрерывно возрастает.
При хорошем тепловом контакте между поверхностями теплообменника и исследуемого объекта, характеризуемом значениями критерия Bi>50, всегда сохраняется постоянное характеристическое значение критериального числа Фурье Fo*, характеризующее связь глубины прогрева с временем от начала подвода тепловой энергии [7]. Вследствие этого толщина прогретого слоя L связана с текущим временем τ следующим соотношением [7]:
Подставив (46) в (45), после преобразований можно получить новое уравнение:
Из уравнения (47) можно определить коэффициент температуропроводности:
В соотношение (48) входят трудноопределимые в эксперименте величины температуры поверхности исследуемого объекта Тсо и скорости изменения температуры этой поверхности Vc. Это обусловлено трудностями в установке измерителя температуры между поверхностями теплообменника и исследуемого объекта таким образом, чтобы он регистрировал только температуру исследуемого объекта, исключая влияние теплообменника на показания измерителя температуры.
Задача значительно упрощается, если учесть результаты численного математического моделирования рассматриваемого процесса нагрева, в соответствии с которыми перепад температур между поверхностями теплообменника и исследуемого объекта в течение всего процесса мал по сравнению с уровнями температур. Это позволяет принять следующее допущение:
Допущение 6:
что в принципе дает возможность ограничиться измерением только температуры теплообменника. Скорость роста температуры при этом определяется из соотношения:
или
В формуле (50) фигурирует среднее значение скорости 
за все время от начала прогрева, а в формуле (51) - мгновенное Vwм, при этом ΔTw - приращение температуры теплообменника в течение выбранного промежутка времени Δτ; индекс «н» соответствует начальному значению; Тн - начальная температура всей измерительной схемы и окружающей среды в самом распространенном случае, когда все начальные температуры совпадают.
Представим формулу для определения плотности теплового потока с учетом (9) и (17б) в виде:
Введем также обозначение
где V0 - имеет размерность скорости [м/с].
С учетом (49), (50), (52) и (53) формулу (48) можно представить в виде:
С целью проверки работоспособности формулы для определения температуропроводности (54) было проведено математическое моделирование, в процессе которого конечно-разностным методом рассчитывалось нестационарное температурное поле в исследуемом объекте, контактирующем с теплообменником, нестационарная температура которого вычислялась по формуле (10).
Расчет проводился в нелинейной постановке с использованием итераций.
Численные исследования проводились при следующих значениях параметров:
В результате численных исследований определено, в частности, что Fo*=0,0665 при ΔТ=0,5 К. С учетом принятых в (55) параметров значение V0=8,2·10-6 м/с.
Подставив в (54) все указанные значения параметров, получим рабочую расчетную формулу:
В формулу (56) входит мгновенное значение Vw=Vwм.
Результаты исследований представлены в таблице 1.
| Таблица 1. | ||||||||||
| Результаты численных расчетов Tw от текущего времени τ, мгновенных скоростей роста температуры теплообменника Vwм и значений коэффициента температуропроводности а, вычисленных по формуле (56). | ||||||||||
| τ, мин | 3 | 5 | 8 | 10 | 13 | 20 | 30 | 40 | 50 | |
| I | τ, с | 180 | 300 | 480 | 600 | 780 | 1200 | 1800 | 2400 | 3000 |
| Тwм, °С | 27,5 | 33 | 38,7 | 41,3 | 44 | 48,3 | 52 | 54,5 | 56,4 | |
| Δτ, с | 180 | 120 | 180 | 120 | 180 | 420 | 600 | 600 | 600 | |
| II | ΔТ, °С | 7,5 | 5,5 | 5,7 | 2,6 | 2,7 | 4,3 | 3,7 | 2,5 | 1,9 |
Vw,  |
0,0417 | 0,0458 | 0,0317 | 0,0217 | 0,015 | 0,01 | 0,0062 | 0,004 | 0,0032 | |
| III | а, 10-7  |
1,23 | 0,5 | 0,47 | 0,64 | 0,82 | 0,87 | 1,0 | 1,39 | 1,43 |
Данные таблицы 1 разделены на три группы (I, II, III). Первую группу составляют данные первичных исследований (расчета, эксперимента). Во вторую группу входят величины, определенные на основании первичных данных и служащие исходными данными для определения мгновенной скорости; при этом Δτn=τn-τn-1; ΔТn=Тn-Тn-1, где n - номер данной колонки таблицы 1. Третья строка группы II содержит данные о вычисленной мгновенной скорости. В третьей группе приведены значения конечной определяемой физической величины (а).
Представленные в таблице 1 значения мгновенной скорости (для простоты в таблице и в формулах (54) и (56) индекс «м» при Vw опущен) определялись по формуле (51).
Из таблицы 1 можно сделать следующие выводы.
1. Скорость роста температуры теплообменника (и, соответственно, нагреваемой поверхности стенки) со временем снижается - от начала процесса до пятидесятой минуты - на порядок.
2. Наибольшее отличие величины коэффициента температуропроводности, определенного по полученной аналитической формуле (54), а в данном конкретном случае - по формуле (56), от исходного значения а=1,28·10-7 м2/с (55), наблюдается на временном отрезке от пятой до двадцатой минуты. Именно в этом временном промежутке, как показали результаты численных расчетов, имеет место максимальное отличие температур теплообменника и прилегающей поверхности исследуемого объекта, и допущение 6, описываемое соотношениями (49), выполняется наименее строго.
3. В целом предложенная методика определения коэффициента температуропроводности, основанная на использовании формулы (54) позволяет получить достаточно точные значения искомой величины, даже несмотря на то, что значения мгновенных скоростей определялись весьма приближенно. Отсюда можно сделать вывод о том, что разработанная физическая модель адекватна физическому процессу.
Повысить точность определения коэффициента температуропроводности можно, вводя необходимые поправки на нелинейность теплового процесса на основе использования результатов численного анализа, а также используя для определения мгновенных скоростей значения производных температуры по времени в каждый заданный момент времени.
Самым эффективным методом компенсации погрешности является применение метода сравнения с эталонным образцом. Этот метод основан на предварительном измерении величины а образца из материала с хорошо известными теплофизическими свойствами.
В последнем случае обработка результатов измерений проводится по формуле:
В (57) индекс «э» относится к эталонному образцу.
При использовании соотношения (57) погрешность определения коэффициента температуропроводности должна существенно снижаться за счет взаимной компенсации систематических погрешностей для эталонного и исследуемого образцов материалов.
Величина коэффициента теплопроводности λ определяется по измеренному значению коэффициента температуропроводности а с помощью соотношения (1). Это возможно потому, что параметры с и γ большинства материалов хорошо известны, в отличие от λ и а.
Таким образом, заявленное устройство обеспечивает быстрое, простое и достоверное определение коэффициента температуропроводности.
Проведем анализ основных закономерностей протекания тепловых процессов в рамках рассмотренной физической модели.
Прежде всего отметим, что благодаря введенному допущению 5, выраженному соотношением (38), удается прийти к координатной зависимости (43), содержащей скорость нагрева Vc. Использование в расчетах величины мгновенной скорости Vсм позволяет исключить необходимость учета существенной нелинейности процесса нагрева, а введение допущения 6 позволяет существенно упростить методику измерений и обработки экспериментальных данных.
Нелинейность процесса проявляется в непрерывном снижении скорости нагрева, что видно из данных таблицы 1. Указанная нелинейность может быть проанализирована с использованием соотношений (37) и (46).
Запишем соотношение (37) в виде:
откуда можно определить толщину прогретого слоя в приближении бесконечно большой теплопроводности:
Подставив в (59) все значения параметров из (55), получим, например, для τ=3000 с: Vw=0,0032 К/с (табл.1), откуда L=0,15 м=15 см.
С другой стороны, тот же расчет по формуле (46) приводит для τ=3000 с к значению L=0,058 м≈6 см, что в 2,5 раза меньше, чем значение, определенное по формуле (59). Такой результат вполне понятен, учитывая низкую теплопроводность и температуропроводность рассматриваемого материала стенки.
Из проведенного анализа можно сделать следующий вывод: для исследуемого объекта толщиной 30 см предложенный способ измерения обеспечивает результат на той стадии, когда исследуемый объект прогревается всего на 1/5 толщины.
Завершим описание сущности заявленного технического решения анализом результатов численных исследований.
На фиг.1 представлены расчетные зависимости перепадов температур между теплообменником и исследуемым объектом от времени от начала процесса нагрева при величинах коэффициента теплопроводности исследуемого объекта в Вт/(м·К): 0,05 (линия L1); 0,1 (линия L2); 0,3 (линия L3); 0,5 (линия L4); 1 (линия L5); 2 (линия L6); 3 (линия L7). Все приведенные линии соответствуют тем же значениям коэффициента температуропроводности, умноженным на 106. Максимальные значения перепада температур реализуются как раз на том временном интервале, в пределах которого погрешности расчетов коэффициента температуропроводности были максимальны (таблица 1).
На фиг.2 представлены начальные участки зависимостей температуры теплообменника от времени нагрева при величинах коэффициента теплопроводности исследуемого объекта в Вт/(м·К): 0,05 (линия L8); 0,1 (линия L9); 0,3 (линия L10); 0,5 (линия L11); 1 (линия L12); 2 (линия L13); 3 (линия L14). Как видно из фиг.2, на начальном временном участке эти зависимости нелинейны, что обусловлено наличием второй степени параметра τ, как это видно из формулы (36). Из фиг.2 следует, что по динамике нагрева можно определять теплофизические свойства материала стенки.
На фиг.3 представлены обратные зависимости - a(Tw) для разных моментов времени от начала поступления теплоносителя в теплообменник, в минутах: 10 (линия L15); 20 (линия L16); 30 (линия L17); 40 (линия L18); 50 (линия L19); 60 (линия L20), а на фиг.4 - аналогичные зависимости а от мощности тепловыделения в теплообменнике при разных моментах времени от начала поступления теплоносителя в теплообменник в минутах: 10 (линия L21); 20 (линия L22); 30 (линия L23); 40 (линия L24); 50 (линия L25); 60 (линия L26). Данные этих фигур могут служить руководством для экспериментального определения коэффициента температуропроводности материала.
Таким образом, заявленное техническое решение может реализовываться как с использованием полученных аналитических формул, так и на основе результатов численных расчетов.
Заявленные способ и устройство по сравнению со способом-прототипом и устройством-прототипом обладают более высокими потребительскими свойствами за счет расширения области применения и повышения точности и достоверности измерений. В отличие от прототипа-способа и прототипа-устройства, в описываемом техническом решении используется протяженный в двух взаимно перпендикулярных направлениях источник тепловой энергии, поэтому измерения дают эффективные (усредненные) значения теплофизических характеристик, а в большинстве случаев на практике нужны не точечные, а именно эффективные значения тепловых характеристик, например при исследовании стен зданий. В заявленных способе и устройстве обеспечивается формирование одномерного теплового потока (тепловой поток распространяется прямолинейно от внутренней поверхности исследуемого объекта к внешней поверхности исследуемого объекта), который описывается значительно более простой математической моделью, чем в прототипе. Таким образом, по сравнению с известными техническими решениями, заявленные способ и устройство обладают более широкой областью применения. В описываемых способе и устройстве используется более простая и более достоверная модель. В известном способе нужно проводить предварительные измерения для определения поправочного коэффициента, а в описываемых технических решениях этого делать не нужно. В уравнении теплового баланса, описывающего работу заявленных способа и устройства, не надо учитывать потери тепловой энергии в окружающую среду за счет конвективного и лучистого теплообмена вследствие того, что применяется протяженный источник тепловой энергии, снабженный тепловой изоляцией. В уравнении теплового баланса описываемого способа отсутствуют необходимые в уравнении теплового баланса способа-прототипа члены, описывающие частичное поглощение тепловой мощности окружающей средой при передаче тепловой энергии от источника энергии к поверхности исследуемого объекта и при частичном отражении от поверхности исследуемого объекта. В модели, описывающей работу заявленного способа, отсутствует необходимость разложения каких-либо функций в ряд Тейлора. Таким образом, описываемые технические решения обладают более высокими точностью и достоверностью. Расширение области применения и повышение точности и достоверности измерений приводит к повышению потребительских свойств заявленных технических решений по сравнению с прототипами.
Сущность изобретения поясняется описанием варианта выполнения заявленного устройства и чертежами, на которых:
- на фиг.1-4 приведены графики, поясняющие сущность изобретения;
- на фиг.5 приведена схема устройства, соответствующая п.5 формулы изобретения.
Устройство для измерения коэффициента температуропроводности содержит (фиг.5) нагреватель 1, первый измеритель температуры 2, второй измеритель температуры 3, входной трубопровод 4, соединительный трубопровод 5, теплообменник 6, выходной трубопровод 7, при этом выход входного трубопровода 4 соединен с входом нагревателя 1, выполненного с возможностью нагрева теплоносителя, выход нагревателя 1 соединен с входом соединительного трубопровода 5, выход которого соединен с входом теплообменника 6, выход теплообменника 6 соединен с входом выходного трубопровода 7, внешняя поверхность теплообменника 6, кроме поверхности, примыкающей к участку поверхности исследуемого объекта 8, снабжена тепловой изоляцией 9, первый измеритель температуры 2 размещен в соединительном трубопроводе 5, а второй измеритель температуры 3 размещен между теплообменником 6 и тепловой изоляцией 9. Третий измеритель температуры 10 размещен в выходном трубопроводе 7. На фиг.5 стрелками 11 и 12 показано направление движения теплоносителя.
Устройство для измерения температуропроводности работает следующим образом. Теплоноситель под давлением по входному трубопроводу 4 поступает в нагреватель 1, где он нагревается, после чего по соединительному трубопроводу 5 поступает в теплообменник 6. Часть своей тепловой энергии теплоноситель передает исследуемому объекту 8. Из теплообменника 6 теплоноситель поступает в выходной трубопровод 7. Первый измеритель температуры 2 измеряет зависимость от времени температуры теплоносителя на входе теплообменника 6. Второй измеритель температуры 3 измеряет зависимость от времени температуры между теплообменником 6 и тепловой изоляцией 9. Температуропроводность определяют из соотношения (54).
Если известны данные по исследованию эталонного объекта, то есть объекта, параметры которого известны с достаточной точностью, то температуропроводность исследуемого объекта может быть определена из соотношения (57).
В расчетные формулы подставляются измеренные значения Tw, Vw, Twэ, Vwэ, измеренные вторым 3 или третьим 10 измерителем температуры. Измерения являются достоверными только при условии равенства показаний второго 3 и третьего 10 измерителей температуры. Если в процессе измерений это условие не выполняется с первого же измерения, необходимо добиться этого равенства регулировкой массового расхода теплоносителя.
ИСТОЧНИКИ ИНФОРМАЦИИ
1. Сухов Р.И., Лебедкин Ю.М., Кузнецов В.Г. и др. Способ бурения скважин и устройство для его осуществления. Патент РФ на изобретение №2237148, приор. 1999.10.06, опубл. 2001.07.20, МПК7 Е21В 6/02, Е21В 7/00, Е21В 10/36.
2. Пилин Б.П., Марков А.А., Молотков С.Л. Способ ультразвуковой дефектоскопии и устройство, его реализующее. Патент РФ на изобретение №2131123, приор. 1996.01.12, опубл. 1999.05.27, МПК6 G01N 29/04.
3. Бобров В.Т., Тарабрин В.Ф., Ордынец С.А., Кулешов Р.В. Ультразвуковой дефектоскоп «Ласточка». Патент РФ на изобретение №2231783, приор. 2001.08.09., опубл. 2003.07.10, МПК7 G01N 29/04.
4. Сергеев В.А. Способ определения теплового сопротивления переход-корпус полупроводниковых приборов. Патент РФ на изобретение №2178893, приор. 2001.03.13, опубл. 2002.01.27, МПК7 G01R 31/26.
5. Сергеев В.А. Устройство для измерения теплового сопротивления транзисторов. Заявка на патент РФ на изобретение №2000127414/09, приор. 2000.10.31, опубл. 2002.10.10, МПК7 G01R 31/26.
6. Медведев В.В., Троицкий О.Ю. Устройство для определения характеристик материалов. Патент РФ на изобретение №2212653, приор. 2002.05.28, опубл. 2003.09.20, МПК7 G01N 25/18.
7. Чернышев В.П., Сысоев Э.В., Павлов Р.В. Способ бесконтактного неразрушающего контроля теплофизических свойств материалов. Патент РФ на изобретение №2251098, приор. 2003.11.17, опубл. 2005.04.27, MПK7 G01N 25/18.
8. Чудинов Ю.В., Ишук И.Н., Фесенко А.И. Термозонд для неразрушающего контроля теплофизических свойств материалов. Заявка на патент РФ на изобретение №2003129494, приор. 02.10.2003, опубл. 27.03.2005, МПК7 С01N 25/18.
9. Лыков А.В. Теория теплопроводности. Высшая школа. - М.: 1967. - 599 с.
Claims (6)
1. Способ измерения коэффициента температуропроводности, заключающийся в нагревании поверхности исследуемого объекта, измерении зависимости от времени температуры первым измерителем температуры, измерении зависимости от времени температуры вторым измерителем температуры, отличающийся тем, что нагревают участок поверхности исследуемого объекта посредством использования теплоносителя, поступающего в теплообменник, внешняя поверхность которого, кроме поверхности, примыкающей к участку поверхности исследуемого объекта, снабжена тепловой изоляцией, при этом измеряют зависимость от времени температуры теплоносителя на входе теплообменника первым измерителем температуры, измеряют зависимость температуры от времени вторым измерителем температуры, размещенным между теплообменником и тепловой изоляцией, и определяют коэффициент температуропроводности а из соотношения
a=4Fo*τ{V0(T0-Tw)/[2(Tw-Tн)Fo*+Vwτ]}2,
где Fo* - характеристическое значение критериального числа Фурье, τ - текущее время, V0=CM/(cγF), С - удельная теплоемкость теплоносителя, М - массовый расход теплоносителя, с - удельная теплоемкость материала исследуемого объекта, γ - плотность материала исследуемого объекта, F - площадь контакта теплообменника и исследуемого объекта, Т0 - температура, измеряемая первым измерителем температуры, Tw - температура, измеряемая вторым измерителем температуры, Тн - температура окружающей среды, Vw - скорость изменения температуры, измеряемой вторым измерителем температуры.
2. Способ измерения коэффициента температуропроводности по п.1, отличающийся тем, что измеряют зависимость от времени температуры теплоносителя на выходе теплообменника третьим измерителем температуры и обеспечивают условия равенства показаний второго и третьего измерителей температуры.
3. Способ измерения коэффициента температуропроводности, заключающийся в нагревании поверхности исследуемого объекта, измерении зависимости от времени температуры первым измерителем температуры, измерении зависимости от времени температуры вторым измерителем температуры, отличающийся тем, что нагревают участок поверхности исследуемого объекта посредством использования теплоносителя, поступающего в теплообменник, внешняя поверхность которого, кроме поверхности, примыкающей к участку поверхности исследуемого объекта, снабжена тепловой изоляцией, при этом измеряют зависимость от времени температуры теплоносителя на входе теплообменника первым измерителем температуры, измеряют зависимость температуры от времени вторым измерителем температуры, размещенным между теплообменником и тепловой изоляцией, причем коэффициент температуропроводности а определяют из соотношения
a=Naэ,
N={cоэ(T0-Tw)[2(Twэ-Tн)Fo*+Vwэτ]}2/{co(T0-Twэ)[2(Twэ-Tн)Fo*+Vwτ]}2,
где аэ - коэффициент температуропроводности эталонного объекта, Ооэ - объемная теплоемкость материала эталонного объекта, Тwэ - температура, измеряемая вторым измерителем температуры при исследовании эталонного объекта, Vwэ - скорость изменения температуры, измеряемой вторым измерителем температуры при исследовании эталонного объекта, cо - объемная теплоемкость материала исследуемого объекта, Fo* - характеристическое значение критериального числа Фурье, τ - текущее время, Т0 - температура, измеряемая первым измерителем температуры, Tw - температура, измеряемая вторым измерителем температуры, Тн - температура окружающей среды.
4. Способ измерения коэффициента температуропроводности по п.3, отличающийся тем, что измеряют зависимость от времени температуры теплоносителя на выходе теплообменника третьим измерителем температуры и обеспечивают условия равенства показаний второго и третьего измерителей температуры.
5. Устройство для измерения коэффициента температуропроводности, содержащее нагреватель, первый измеритель температуры, второй измеритель температуры, отличающееся тем, что оно дополнительно содержит входной трубопровод, соединительный трубопровод, теплообменник, выходной трубопровод, при этом выход входного трубопровода соединен с входом нагревателя, выполненного с возможностью нагрева теплоносителя, а выход нагревателя соединен с входом соединительного трубопровода, выход которого соединен с входом теплообменника, выход теплообменника соединен с входом выходного трубопровода, внешняя поверхность теплообменника, кроме поверхности, примыкающей к участку поверхности исследуемого объекта, снабжена тепловой изоляцией, первый измеритель температуры размещен в соединительном трубопроводе, а второй измеритель температуры размещен между теплообменником и тепловой изоляцией.
6. Устройство для измерения коэффициента температуропроводности по п.5, отличающееся тем, что оно дополнительно содержит третий измеритель температуры, размещенный в выходном трубопроводе.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2006116343/28A RU2322662C2 (ru) | 2006-05-06 | 2006-05-06 | Способ измерения коэффициента температуропроводности (варианты) и устройство для его осуществления |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2006116343/28A RU2322662C2 (ru) | 2006-05-06 | 2006-05-06 | Способ измерения коэффициента температуропроводности (варианты) и устройство для его осуществления |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2006116343A RU2006116343A (ru) | 2007-12-10 |
| RU2322662C2 true RU2322662C2 (ru) | 2008-04-20 |
Family
ID=38903186
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2006116343/28A RU2322662C2 (ru) | 2006-05-06 | 2006-05-06 | Способ измерения коэффициента температуропроводности (варианты) и устройство для его осуществления |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2322662C2 (ru) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2722088C1 (ru) * | 2019-05-06 | 2020-05-26 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Самарский государственный университет путей сообщения" (СамГУПС) | Способ измерения удельного теплового сопротивления и устройство для его осуществления |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113848231B (zh) * | 2020-06-28 | 2024-03-08 | 宝山钢铁股份有限公司 | 基于炼焦煤热解过程中热扩散率的结焦性判断方法 |
-
2006
- 2006-05-06 RU RU2006116343/28A patent/RU2322662C2/ru not_active IP Right Cessation
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2722088C1 (ru) * | 2019-05-06 | 2020-05-26 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Самарский государственный университет путей сообщения" (СамГУПС) | Способ измерения удельного теплового сопротивления и устройство для его осуществления |
Also Published As
| Publication number | Publication date |
|---|---|
| RU2006116343A (ru) | 2007-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Park et al. | Experimental study on condensation heat transfer in vertical minichannels for new refrigerant R1234ze (E) versus R134a and R236fa | |
| Rose | Heat-transfer coefficients, Wilson plots and accuracy of thermal measurements | |
| Del Col et al. | Effect of cross sectional shape during condensation in a single square minichannel | |
| Lu et al. | Inverse estimation of the inner wall temperature fluctuations in a pipe elbow | |
| Kim et al. | Experimental and numerical analysis of heat transfer phenomena in a sensor tube of a mass flow controller | |
| Piccolo et al. | Experimental study of heat transfer characteristics of finned-tube and circular-pore heat exchangers in oscillatory flow | |
| Pătrăşcioiu et al. | Prediction of the outlet temperatures in triple concentric—tube heat exchangers in laminar flow regime: case study | |
| Bertsche et al. | Experimental investigation on heat transfer in laminar, transitional and turbulent circular pipe flow with respect to flow regime boundaries | |
| Du et al. | Tikhonov regularization stabilizes multi-parameter estimation of geothermal heat exchangers | |
| Liu et al. | Experimental study on forced convective heat transfer characteristics in quartz microtube | |
| RU2322662C2 (ru) | Способ измерения коэффициента температуропроводности (варианты) и устройство для его осуществления | |
| Shah et al. | Estimation of local heat flux with CFD and enhanced conjugate gradient method for laminar and turbulent flow in a helical coil tube heat exchanger | |
| RU2344338C1 (ru) | Способ определения толщины отложений на внутренней поверхности трубопроводов | |
| Kamsanam et al. | Development of experimental techniques for measurement of heat transfer rates in heat exchangers in oscillatory flows | |
| JP3926571B2 (ja) | 固体の温度伝導率等測定方法および装置 | |
| RU60730U1 (ru) | Устройство для измерения коэффициента температуропроводности | |
| RU2568983C1 (ru) | Способ определения коэффициента теплопроводности жидкой тепловой изоляции в лабораторных условиях | |
| Stadnyk et al. | Features of heat transfer in the environment when it is sprayed with rotary rollers. | |
| RU2330270C2 (ru) | Способ измерения удельного теплового сопротивления и устройство для его осуществления | |
| Samadi et al. | Analytical solution for partial heating on the exterior of the pipe with application to measuring fluid flow rate | |
| Hedayati-Dezfooli et al. | A design of experimental apparatus for studying coupled heat and moisture transfer in soils at high-temperature conditions | |
| Minakov et al. | Experiment-calculated investigation of the forced convection of nanofluids using single fluid approach | |
| RU54193U1 (ru) | Устройство для измерения теплофизических характеристик (варианты) | |
| RU2646437C1 (ru) | Способ определения коэффициента теплопроводности жидкой тепловой изоляции при нестационарном тепловом режиме | |
| Hanby et al. | Modeling the dynamic response of conduits |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20110507 |