RU2556290C1 - Способ определения теплофизических свойств твердых материалов - Google Patents
Способ определения теплофизических свойств твердых материалов Download PDFInfo
- Publication number
- RU2556290C1 RU2556290C1 RU2013154721/28A RU2013154721A RU2556290C1 RU 2556290 C1 RU2556290 C1 RU 2556290C1 RU 2013154721/28 A RU2013154721/28 A RU 2013154721/28A RU 2013154721 A RU2013154721 A RU 2013154721A RU 2556290 C1 RU2556290 C1 RU 2556290C1
- Authority
- RU
- Russia
- Prior art keywords
- temperature
- sample
- heat
- change
- time
- Prior art date
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Изобретение относится к теплофизическим измерениям и может быть использовано в теплофизическом приборостроении. Способ заключается в нагревании одной из поверхностей образца или ее участка до максимальной температуры, которую поддерживают до момента времени, когда измеряемая разность температур на границах исследуемого участка уменьшится до заданного значения. Затем снижают мощность нагрева до наступления стадии остывания образца с измеряемой скоростью изменения температуры, не превышающей установленного значения. Измеряют удельное количество тепла, приращения температуры на границах за два последовательных интервала времени, первый из которых определяется моментами времени: подачи тепла и достижения установленной скорости изменения температуры. Теплопроводность и объемную теплоемкость определяют по формулам. Технический результат - увеличение точности определения теплофизических свойств и уменьшение времени измерения. 4 ил.
Description
Настоящее изобретение относится к теплофизическим измерениям, в частности измерению теплопроводности и объемной теплоемкости твердых материалов, и может быть использовано в теплофизическом приборостроении.
Из существующего уровня техники известен способ определения твердых материалов, в котором измерение теплопроводности осуществляется в стационарном режиме, а теплоемкости в переходной области до наступления стационарного режима из уравнения теплового баланса. В начальной стадии опыта плоский образец находится в контакте с термостатируемой средой и имеет температуру t1. Затем противоположная грань образца приводится в контакт с другой термостатируемой средой с температурой t2. Измеряют тепловые потоки, действующие на противоположных гранях образца до установления стационарного режима τy, и среднеобъемную температуру в стационарном режиме 0,5(t1+t2), а теплоемкость определяют по известному уравнению теплового баланса. В стационарном режиме определяют теплопроводность по известному уравнению (см. Теплофизические измерения и приборы / под ред. Е.С. Платунова. Л.: Машиностроение, 1986. С.149).
Недостатками данного технического решения являются большое время измерения, определяемое наступлением стационарного режима, низкая точность и сложность технической реализации вследствие необходимости измерения двух нестационарных тепловых потоков в переходной области до наступления стационарного режима.
Наиболее близким к заявленному техническому решению является способ определения теплофизических свойств жидкообразных и мелкодисперсных сред, включающий воздействие одномерным тепловым потоком на одну из поверхностей образца, измерение разности температур, количества тепла, поступившего в образец, приращения температуры в точках образца за одни и те же интервалы времени, и по измеренным параметрам определяют искомые свойства (см. а.с. СССР №1017985, G01N 25/18, 1981 г.).
Недостатком данного технического решения является низкая точность, обусловленная большой погрешностью определения интегрального параметра
по двум точкам на границах образца или его участка в переходной области теплопередачи. Второй недостаток - это невыполнимость условия выбора длительности первого интервала [0, τ1] по всем значениям из диапазона измерения искомых величин, что позволяет найти два других интервала [τ1, τ2] и [τ1, τ3], вследствие недостижимости условия
, когда справедлива формула для определения теплопроводности, или условия
, когда справедлива формула для определения объемной теплоемкости, по причине наступления стационарного режима. Это приводит к тому, что измерение одной из величин невозможно, что, при технической реализации, приведет к неограниченному времени измерения.
Задачей, на решение которой направлено заявляемое изобретение, является повышение точности определения теплофизических свойств и уменьшение времени измерения.
Поставленный технический результат достигается тем, что в заявленном способе, включающем нагревание одной из поверхностей образца или ее участка, измерение температуры на границах исследуемого участка и их разности, удельного количества тепла, поступившего в образец, дополнительно измеряют скорость изменения температуры на тепловоспринимающей границе, нагрев осуществляют до максимальной температуры и поддерживают ее до момента времени, когда разность температур уменьшится до заданного значения, снижают мощность нагрева до наступления стадии остывания образца со скоростью, не превышающей установленного значения, измеряют изменение указанных величин за два последовательных интервала времени, первый из которых определяется моментами времени: подачи тепла и достижения установленной скорости изменения температуры, а искомые величины определяют по формулам
где λ, С - теплопроводность и объемная теплоемкость;
Δτ1, Δτ2 - соответственно, интервалы времени:[0, τ1] и [τ1, τ2];
t(0,τ), t(L,τ) - температура на границах участка исследования образца;
р0, р1 - весовые коэффициенты;
Уменьшение разности температур и уменьшение скорости изменения температуры до определенного значения позволяет увеличить точность определения интегрального параметра
. Нагрев до максимальной температуры и поддержание ее в трубке данного значения позволяет уменьшить время достижения установленной разности температур и, следовательно, уменьшить время измерения. Использование стадии нагревания и остывания образца обеспечивает получение системы линейных уравнений с малым числом обусловленности, при котором имеет место незначительное увеличение погрешности определения искомых тепловых величин по сравнению с погрешностью исходных данных.
Для обоснования повышения точности комплексного определения теплопроводности и объемной теплоемкости необходимо рассмотреть основные положения теории, на которой базируется заявленный способ. Он основан на математическом описании теплопередачи в объекте исследования в виде интегральной формы уравнения теплопроводности. Для одномерного варианта температурного поля ее можно записать в следующем виде:
где L, t(0,τ), t(L,τ) - длина исследуемого участка образца и температура на его границах;
- удельное количество тепла, поступившего в образец за интервал времени [0,τ]; λ, С - теплопроводность и объемная теплоемкость.
Для определения повторного интеграла используется приближенная формула
, основанная на интерполяционном многочлене Лагранжа и которая является точной для многочлена степени m, где
,
, (i≠k) - весовой коэффициент; t(xk,τ) - температура в точке xk; m+1 - количество точек на интервале [0, L].
Для упрощения технической реализации способа необходимо использовать две точки измерения температуры на границах интервала [0, L]. Тогда приближенная формула будет точной для линейного распределения температуры и имеет вид
. Погрешность данной формулы для других распределений можно оценить из выражения: Δ(τ)<kL4 a -1(мах dt(τ)/dτ), где k - коэффициент, L - длина исследуемого участка образца; а - температуропроводность образца, мах dt(τ)/dτ) - максимальная скорость изменения температуры. Поэтому для уменьшения относительной погрешности определения повторного интеграла необходимо уменьшить скорость изменения температуры и относительную разность температур [t(0,τ)-t(L,τ)]/[t(0,τ)+t(L,τ)].
Если зафиксировать два интервала времени [0, τ1] и [τ1, τ2], то получим систему двух уравнений с неизвестными λ и С. Использование стадии нагревания и остывания образца на интервале [0, τ1] и остывания на интервале [τ1, τ2] позволяет получить систему двух линейных уравнений с малым числом обусловленности, которое обеспечивает незначительное увеличение погрешности определения искомых тепловых величин по сравнению с погрешностью исходных данных.
Сущность изобретения поясняется чертежами, на которых изображено следующее.
На фиг.1 - схема измерительной ячейки для измерения комплекса: теплопроводность - объемная теплоемкость низкотеплопроводных материалов на образцах цилиндрической формы.
На фиг.2 - схема измерительной ячейки для измерения комплекса: теплопроводность - объемная теплоемкость высокотеплопроводных материалов на образцах в виде пластины.
На фиг.3 - изменение температуры на тепловоспринимающей поверхности образца (), разности температур () и плотности теплового потока () во времени, полученных при имитационном моделировании комплексного измерения λС низкотеплопроводных материалов
На фиг.4 - изменение температуры на тепловоспринимающем участке поверхности образца (), разности температур () и плотности теплового потока () во времени, полученных при имитационном моделировании комплексного измерения λС высокотеплопроводных материалов
Примеры реализации способа для измерения теплофизических свойств низкотеплопроводных и высокотеплопроводных материалов показаны соответственно на фиг.1 и фиг.2. На фиг.3 приведены графики, поясняющие заявленный способ, которые были получены на модели измерительной ячейки, приведенной на фиг.1: λ=0,05 Вт/(м·К), а=10-7 м2/с; L=4 мм; τ1, τ2 - моменты времени окончания, соответственно, первого и второго интервала интегрирования, измерения приращений температур и количества тепла, поступившего в образец. На фиг.4 приведены графики, поясняющие заявленный способ, полученные на модели измерительной ячейки, приведенной на фиг.2: λ=7 Вт/(м·К), а=3,3·10-6 м2/с; L=10 мм, толщина H=1 мм.
Способ осуществляется следующим образом. Образец 1 нагревают до максимальной установленной температуры под воздействием нагревателя 3 и поддерживают в трубке данного значения до установления заданной разности температур. Для низкотеплопроводного материала, как показано на фиг.3, заданная разность установилась после подачи второго импульса тепла. Для высокотеплопроводных материалов данное условие выполняется после достижения максимальной температуры. После уменьшения теплового потока до определенной величины, в частном случае до нуля, образец охлаждают со скоростью, удовлетворяющей заданному значению. Данный момент времени обозначен τ1 и является окончанием первого интервала измерения. Дальнейшее остывание образца происходит в течение заданного интервала [τ1, τ2] времени, который выбирается исходя из достаточной точности определения искомых величин. В реальных условиях измерения он составляет 30÷40 с. Имитационное измерение теплопроводности и объемной теплоемкости на ячейке, представленной на фиг.1, при погрешности в измерительной информации, соответствующей погрешности вычисления температурного поля в системе Маткад и методической погрешности определения количества тепла, показало, что максимальная погрешность определения теплоемкости не превышает 2%, а теплопроводности - 1% для начальной температуры 0°С, максимальной температуры 25°С, разности температур меньше 10°С, интервала [τ1, τ2] больше 30 с.
Погрешности имитационного измерения теплопроводности и объемной теплоемкости на ячейке, показанной на фиг.2, при аналогичных условиях и наличии контактного теплового сопротивления Rк=1,9·10-4 м2·К/Вт между образцом: нагревателем и опорами, на торцах которых закреплены спаи термопар, составили: δλ<±2%, δС<±8%. Относительно большая погрешность, по сравнению с предыдущим результатом, объясняется присутствием погрешности измерения температур вследствие наличия теплового сопротивления в зоне контакта между образцом и термочувствительными элементами.
Заявленный способ прошел теоретические исследования на различных моделях измерительных ячеек методом имитационного моделирования с целью использования в приборе для комплексного измерения теплопроводности и объемной теплоемкости.
Claims (1)
- Способ определения теплофизических свойств твердых материалов, включающий нагревание одной из поверхностей образца или ее участка, измерение температуры на границах исследуемого участка и их разности, удельного количества тепла, поступившего в образец, отличающийся тем, что дополнительно измеряют скорость изменения температуры на тепловоспринимающей границе, нагрев осуществляют до максимальной температуры и поддерживают ее до момента времени, когда разность температур уменьшится до заданного значения, снижают мощность нагрева до наступления стадии остывания образца со скоростью, не превышающей установленного значения, измеряют изменение указанных величин за два последовательных интервала времени, первый из которых определяется моментами времени: подачи тепла и достижения установленной скорости изменения температуры, а искомые величины определяют по формулам
где λ, С - теплопроводность и объемная теплоемкость;
Δt1, Δτ2 - соответственно, интервалы времени: [0, τ1] и [τ1, τ2];
t(0,τ), t(L,τ) - температура на границах участка исследования образца;
p0, p1 - весовые коэффициенты;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013154721/28A RU2556290C1 (ru) | 2013-12-11 | 2013-12-11 | Способ определения теплофизических свойств твердых материалов |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013154721/28A RU2556290C1 (ru) | 2013-12-11 | 2013-12-11 | Способ определения теплофизических свойств твердых материалов |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2013154721A RU2013154721A (ru) | 2015-06-20 |
RU2556290C1 true RU2556290C1 (ru) | 2015-07-10 |
Family
ID=53433455
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2013154721/28A RU2556290C1 (ru) | 2013-12-11 | 2013-12-11 | Способ определения теплофизических свойств твердых материалов |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2556290C1 (ru) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2758414C1 (ru) * | 2020-12-10 | 2021-10-28 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский авиационный институт (национальный исследовательский университет)" | Устройство для определения комплекса теплофизических характеристик композиционных материалов |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1017985A1 (ru) * | 1981-12-29 | 1983-05-15 | Новомосковский Филиал Московского Ордена Ленина И Ордена Трудового Красного Знамени Химико- Технологического Института Им.Д.И.Менделеева | Способ определени теплофизических свойств жидкообразных и мелкодисперсных сред и устройство дл его осуществлени |
SU1163232A1 (ru) * | 1982-07-28 | 1985-06-23 | Гомельский Государственный Университет | Способ определени коэффициента температуропроводности материалов |
RU2149387C1 (ru) * | 1996-10-08 | 2000-05-20 | Клебанов Михаил Геннадиевич | Способ неразрушающего контроля теплофизических характеристик материалов |
RU2255329C1 (ru) * | 2003-12-22 | 2005-06-27 | Тамбовский военный авиационный инженерный институт | Способ определения теплофизических свойств материалов |
JP2005227010A (ja) * | 2004-02-10 | 2005-08-25 | Nichias Corp | 熱伝導率測定装置及び熱伝導率測定方法 |
RU2328725C1 (ru) * | 2006-11-17 | 2008-07-10 | Государственное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" (ГОУ ВПО "ТГТУ") | Способ неразрушающего определения теплофизических свойств твердых материалов |
-
2013
- 2013-12-11 RU RU2013154721/28A patent/RU2556290C1/ru not_active IP Right Cessation
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1017985A1 (ru) * | 1981-12-29 | 1983-05-15 | Новомосковский Филиал Московского Ордена Ленина И Ордена Трудового Красного Знамени Химико- Технологического Института Им.Д.И.Менделеева | Способ определени теплофизических свойств жидкообразных и мелкодисперсных сред и устройство дл его осуществлени |
SU1163232A1 (ru) * | 1982-07-28 | 1985-06-23 | Гомельский Государственный Университет | Способ определени коэффициента температуропроводности материалов |
RU2149387C1 (ru) * | 1996-10-08 | 2000-05-20 | Клебанов Михаил Геннадиевич | Способ неразрушающего контроля теплофизических характеристик материалов |
RU2255329C1 (ru) * | 2003-12-22 | 2005-06-27 | Тамбовский военный авиационный инженерный институт | Способ определения теплофизических свойств материалов |
JP2005227010A (ja) * | 2004-02-10 | 2005-08-25 | Nichias Corp | 熱伝導率測定装置及び熱伝導率測定方法 |
RU2328725C1 (ru) * | 2006-11-17 | 2008-07-10 | Государственное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" (ГОУ ВПО "ТГТУ") | Способ неразрушающего определения теплофизических свойств твердых материалов |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2758414C1 (ru) * | 2020-12-10 | 2021-10-28 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский авиационный институт (национальный исследовательский университет)" | Устройство для определения комплекса теплофизических характеристик композиционных материалов |
RU2758414C9 (ru) * | 2020-12-10 | 2022-02-03 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский авиационный институт (национальный исследовательский университет)" | Устройство для определения комплекса теплофизических характеристик композиционных материалов |
Also Published As
Publication number | Publication date |
---|---|
RU2013154721A (ru) | 2015-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103411996B (zh) | 固体材料导热系数测量装置及测量方法 | |
Warzoha et al. | Determining the thermal conductivity of liquids using the transient hot disk method. Part II: Establishing an accurate and repeatable experimental methodology | |
CN104833695B (zh) | 基于红外热成像技术的金属薄板热导率测量方法 | |
Moreira et al. | Heat transfer coefficient: a review of measurement techniques | |
CN103954650B (zh) | 一种固体材料热扩散系数测试方法与系统 | |
Hubble et al. | A hybrid method for measuring heat flux | |
RU2387981C1 (ru) | Способ комплексного определения теплофизических характеристик материалов | |
CN106124078A (zh) | 一种采用双热电偶测量强瞬变流体温度的方法 | |
Mahvi et al. | Accurate and inexpensive thermal time-of-flight sensor for measuring refrigerant flow in minichannels | |
Lin et al. | Inverse problem of unsteady conjugated forced convection in parallel plate channels | |
CN104749214A (zh) | 一种基于瞬态平面热源法测量液体导热系数的恒温热浴装置 | |
Pan et al. | Identification of temperature-dependent thermal conductivity and experimental verification | |
Bhattacharya et al. | Characterization of the temperature oscillation technique to measure the thermal conductivity of fluids | |
RU2556290C1 (ru) | Способ определения теплофизических свойств твердых материалов | |
Zandt et al. | Capabilities for dielectric-constant gas thermometry in a special large-volume liquid-bath thermostat | |
Wu et al. | Quantitative evaluation of the natural convection effect on thermal conductivity measurement with transient plane source method | |
CN100437098C (zh) | 热传导性能检测装置及检测方法 | |
CN103245692B (zh) | 基于稳态分析的半球向全发射率与导热系数的测量方法 | |
CN102998333B (zh) | 一种表面传热系数的测量装置及方法 | |
CN202947983U (zh) | 表面传热系数的测量装置 | |
Silin et al. | Heat transfer in a short parallel-plate channel in the transition regime | |
Van der Westhuizen et al. | Implementation of liquid crystal thermography to determine wall temperatures and heat transfer coefficients in a tube-in-tube heat exchanger | |
CN103267772B (zh) | 基于瞬态分析的大温差样品的半球向全发射率测量方法 | |
RU2523090C1 (ru) | Способ определения удельной теплоемкости материалов | |
RU2631007C1 (ru) | Теплосчетчик на основе накладных датчиков |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20151212 |