RU2374631C2 - Способ определения комплекса теплофизических свойств твердых материалов - Google Patents
Способ определения комплекса теплофизических свойств твердых материалов Download PDFInfo
- Publication number
- RU2374631C2 RU2374631C2 RU2008102328/28A RU2008102328A RU2374631C2 RU 2374631 C2 RU2374631 C2 RU 2374631C2 RU 2008102328/28 A RU2008102328/28 A RU 2008102328/28A RU 2008102328 A RU2008102328 A RU 2008102328A RU 2374631 C2 RU2374631 C2 RU 2374631C2
- Authority
- RU
- Russia
- Prior art keywords
- excess temperature
- temperature
- time
- sample
- thermal
- Prior art date
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Изобретение относится к теплофизике. Способ включает тепловое импульсное воздействие на плоскую поверхность исследуемого образца и измерение избыточной температуры образца в одной точке в заданном интервале времени. Тепловое импульсное воздействие осуществляют лучистым тепловым потоком известной плотности и длительности, а измерение избыточной температуры с момента подачи теплового импульса проводят в центральной части нагреваемой поверхности образца, при этом регистрируют значение максимальной избыточной температуры и время ее достижения. Для данного теплового процесса используют математическую модель прямой задачи теплопроводности на основе метода конечных разностей и решают обратную коэффициентную задачу теплопроводности вариационным методом на интервале наблюдения избыточной температуры. Технический результат - повышение точности определения комплекса теплофизических свойств твердых материалов. 4 ил.
Description
Изобретение относится к теплофизическим измерениям. Область применения - определение теплофизических свойств (ТФС) материалов и изделий неразрушающим методом.
Известен способ идентификации комплекса теплофизических свойств твердых материалов, включающий воздействие тепловыми импульсами от линейного источника на плоскую поверхность исследуемого и эталонного образцов, измерение избыточных температур в момент подачи тепловых импульсов в точках, расположенных на фиксированных расстояниях от линии нагрева на поверхности образцов, по идентифицированным параметрам теплофизических свойств образцов и действительным значениям теплофизических свойств эталона находят искомый комплекс теплофизических свойств (патент РФ №2125258, 1999 г., МПК: G01N 25/18 (2006.01)). Под избыточной температурой понимается температура, отсчитываемая от начальной температуры, при которой находился образец в момент подачи первого теплового импульса. Под эталонным образцом понимается образец материала с известными теплофизическими свойствами.
Известен также способ идентификации комплекса теплофизических свойств твердых материалов, включающий тепловое импульсное воздействие от линейного источника нагрева на плоскую поверхность исследуемого и эталонного образцов, измерение избыточной температуры на плоской поверхности образцов на фиксированном расстоянии от линии нагрева с момента подачи теплового импульса, причем тепловое импульсное воздействие и измерение избыточной температуры осуществляют в плоскости контакта исследуемого и эталонного образцов, а измерение избыточной температуры производят в одной точке в заданном интервале времени, используя математическую модель прямой задачи теплопроводности на основе метода конечных разностей (решение о выдаче патента на изобретение по заявке 2005117325/28 от 06.06.2005, МПК: G01N 25/00 (2006.01)). Этот способ принят за ближайший аналог.
Недостатками способов являются низкая точность измерений и длительное время проведения измерений.
Техническим результатом изобретения является повышение точности определения комплекса теплофизических свойств твердых материалов и упрощение теплофизических измерений.
Технический результат достигается тем, что в способе определения комплекса теплофизических свойств твердых материалов, заключающемся в тепловом импульсном воздействии на плоскую поверхность исследуемого образца, измерении избыточной температуры на плоской поверхности образца с момента подачи теплового импульса в одной точке в заданном интервале времени, причем тепловое импульсное воздействие на плоскую поверхность образца осуществляют лучистым тепловым потоком известной плотности и длительности, а измерение избыточной температуры с момента подачи теплового импульса проводят в центральной части нагреваемой поверхности образца, при этом регистрируют значение максимальной избыточной температуры и время ее достижения, используя математическую модель прямой задачи теплопроводности на основе метода конечных разностей, применяют промежуточную сеточную функцию источника нагрева излучением и конвективного теплообмена:
для учета нелинейной зависимости теплофизических свойств материала от температуры используются кусочно-линейные функции (принято условное разбиение шкалы избыточной температуры на три участка [0;Т1],(Т1;Т2],(Т2,∞)):
решают обратную коэффициентную задачу теплопроводности вариационным методом на интервале наблюдения избыточной температуры [τ1,τ2], при выборе параметра регуляризации по невязке:
искомые теплофизические свойства а, λ находят из минимума функции Jξ(а,λ),
где Тk m,n,p - сеточная функция;
Т - избыточная температура, регистрируемая в эксперименте;
Tmax - максимальное значение избыточной температуры, регистрируемое в эксперименте;
Т1 - избыточная температура, регистрируемая в момент времени τ1,
Т2 - избыточная температура, регистрируемая в момент времени τ2;
τmax - момент времени регистрации избыточной температуры Тmax;
а - температуропроводность;
λ - теплопроводность;
x, y, z - оси правой декартовой прямоугольной системы координат;
k - номер отсчета по времени;
h - шаг сетки по расстоянию;
Δτ - шаг сетки по времени;
q - плотность теплового потока;
α - коэффициент теплоотдачи;
К1а,К2а,К3а - линейные коэффициенты зависимости температуропроводности от температуры;
K1λ,K2λ,K3λ - линейные коэффициенты зависимости теплопроводности от температуры;
τ - текущее время, отсчитываемое от момента подачи импульса;
τ* - длительность теплового импульса;
ξ1,ξ2 - параметры регуляризации;
J - функционал невязки.
Сущность изобретения поясняется чертежами.
На фиг.1 представлена пространственная сетка (M,N,P) математической модели в плоскости (x,y,z).
На фиг.2 представлены относительные графики, аппроксимирующие зависимость: 1 - для температуропроводности в виде а(T)/а(0) при К1а=-0,001·10-3, K2а=-0,01·10-3;
K3а=-0,05·10-3; 2 - для теплопроводности в виде λ(T)/λ(0) при К1λ=0,1·10-3, K2λ=4·10-3, K3λ=8·10-3.
На фиг.3 приведен вариант определения теплофизических свойств исследуемого материала с минимально возможным значением невязки, где кривая 1 - график температуры, построенный математической моделью, кривая 2 - график температуры реального испытания, где τ1 - время начала измерения температуры, отсчитываемое от момента подачи теплового импульса; τ2 - время окончания измерения температуры, отсчитываемое от момента подачи теплового импульса; Δτ - шаг сетки по времени.
На фиг.4 показана схема устройства, реализующая предлагаемый способ определения комплекса ТФС твердых материалов.
Устройство (Фиг.4) содержит инфракрасный излучатель 1, инфракрасный измеритель температуры 2 и исследуемый материал 3.
Инфракрасный излучатель 1 в течение заданного времени τ* равномерно воздействует по всей площади поверхности исследуемого материала радиусом R, а с помощью инфракрасного измерителя температуры в интервале времени [τ1,τ2] B центре материала измеряют избыточные значения температур. Для соблюдения условия полуограниченности исследуемого материала толщиной Н необходимо, чтобы H/R≥20 (Карслоу Г., Егер Д. Теплопроводность твердых тел. - М.: Наука, 1964. - 114 с.).
Сущность изобретения заключается в следующем.
По всей площади плоской поверхности исследуемого образца 3 осуществляют тепловое импульсное воздействие источником тепла 1 (инфракрасным излучателем), формирующим лучистый тепловой поток известной плотности и длительности, при этом измеряют инфракрасным измерителем температуры 2 избыточную температуру в центре образца в заданном интервале времени [τ1,τ2] (фиг.3, фиг.4). Охлаждение поверхности осуществляется путем естественного конвективного теплообмена с окружающей средой. Для определения теплофизических свойств исследуемого материала используют математическую модель прямой задачи теплопроводности на основе метода конечных разностей (фиг.1), полученную решением нелинейной задачи теплопроводности с учетом теплоотдачи с окружающей средой:
где Tn - температура поверхности исследуемого образца; h (τ)- ступенчатая функция.
Нелинейную задачу теплопроводности (7) при граничных условиях (8) решают методом конечных разностей. Разностная схема для трехмерного уравнения теплопроводности имеет вид (3). Граничные условия аппроксимируют на основании выражения:
полученного с помощью метода теплового баланса. При этом q и α отнесены к единице площади поперечного сечения пространственной сетки и выражены в Вт/м2 и Дж/(м2·К).
Выражение (9) приводят к виду (1), (2) и рассчитывают изменение значения избыточной температуры во времени в заданном интервале наблюдения избыточной температуры [τ1,τ2] с учетом зависимости температуропроводности и теплопроводности от температуры согласно кусочно-линейным функциям (4), (5).
Для расчета искомой сеточной функции «прогонками вдоль осей x,y,z», соответствующей температуре на k+1 временном слое, сначала применяют выражения (1), (2) при далее (3) при
Для получения при численном решении физически оправданных результатов следует соблюдать условие сходимости вычислений:
При минимальном значении функционала невязки значения избыточных температур, рассчитанные математической моделью и полученные при реальном испытании, совпадают с минимальной погрешностью, следовательно, значениям теплофизических свойств исследуемого материала ставятся в соответствие значения λ, и а математической модели (фиг.2, фиг.3). Расчет невязки с использованием дополнительных слагаемых и в формуле (6) позволяет достичь наилучшего приближения дискретной модели к экспериментальной термограмме, так как оно учитывает условие существования слабого экстремума. Решение задачи определения ТФС производится каким-либо из вариационных методов исчисления.
Результаты предварительных экспериментов показали, что при мощности теплового потока 150 Вт/м2 и времени τ*=60 сек максимальное значение избыточной температуры для класса теплоизоляционных материалов составляет 50…120°C, что требует необходимость учета нелинейной зависимости изменения ТФС материла от температуры. Так с ростом температуры теплопроводность увеличивается, а температуропроводности уменьшается, а доля изменения ТФС может составлять до 10-30% на 100°С избыточной температуры, при этом наилучшая аппроксимация зависимости ТФС от температуры достигается параболической, экспоненциальной или кусочно-линейной функциями (Платунов Е.С. и др. Теплофизические измерения и приборы - Л.: Машиностроение, 1986).
Для обеспечения высокой воспроизводимости при каждом измерении на основании закона Ламберта угол наблюдения для металлов должен быть 0…40°, а для диэлектриков в интервале 0…60° (Инфракрасная термография. Основы, техника, применение: Пер. с франц. - М.: Мир, 1988. - С.47). Инфракрасный излучатель должен быть расположен параллельно исследуемым материалам строго на определенной высоте для обеспечения заданного значения мощности теплового потока.
Для учета возникновения вынужденной конвекции при проведении измерений в полевых условиях коэффициент теплоотдачи может быть записан как функция от скорости ветра (Справочник физических величин / Под ред. проф. Г.А.Рябинина. - СПб., Лениздат; Издательство «Союз», 2001. - С.99).
Применение предлагаемого способа позволяет повысить точность определения комплекса ТФС материалов по сравнению с прототипом за счет исключения составляющих частных погрешностей измерения расстояния от источника тепла до датчика температуры, теплофизических свойств материала эталона, а также за счет учета зависимости теплофизических свойств исследуемого материла от температуры в виде кусочно-линейных функций, более точно аппроксимирующих реальный физический процесс.
В связи с тем, что предлагаемый способ позволяет определять комплекс ТФС без использования материала эталона, на котором размещают линейный источник на фиксированном расстоянии от датчика температуры, а измерение температуры производят в любой точке в центре исследуемого материала, проведение теплофизических измерений упрощается.
Claims (1)
- Способ определения комплекса теплофизических свойств твердых материалов, включающий тепловое импульсное воздействие на плоскую поверхность исследуемого образца, измерение избыточной температуры на плоской поверхности образца с момента подачи теплового импульса в одной точке в заданном интервале времени, отличающийся тем, что тепловое импульсное воздействие на плоскую поверхность образца осуществляют лучистым тепловым потоком известной плотности и длительности, а измерение избыточной температуры с момента подачи теплового импульса проводят в центральной части нагреваемой поверхности образца, при этом регистрируют значение максимальной избыточной температуры и время ее достижения, используя математическую модель прямой задачи теплопроводности на основе метода конечных разностей, применяют промежуточную сеточную функцию источника нагрева излучением и конвективного теплообмена:
, при τ≤τ*;
, при τ>τ*,
при , ;
сеточную функцию, учитывающую теплоперенос по всему объему образца:
,
при , , ;
для учета нелинейной зависимости теплофизических свойств материала от температуры используют кусочно-линейные функции:
решают обратную коэффициентную задачу теплопроводности вариационным методом на интервале наблюдения избыточной температуры [τ1,τ2], при выборе параметра регуляризации по невязке:
,
искомые теплофизические свойства а, 1 находят из минимума функции Jξ(a,λ), где - сеточная функция;
Т - избыточная температура, регистрируемая в эксперименте;
- избыточная температура, рассчитанная математической моделью;
Tmax - максимальное значение избыточной температуры, регистрируемое в эксперименте;
- максимальное значение избыточной температуры, рассчитанное математической моделью;
Т1 - избыточная температура, регистрируемая в момент времени τ1,
Т2 - избыточная температура, регистрируемая в момент времени τ2;
τmax - момент времени регистрации избыточной температуры Тmax;
- момент времени, рассчитанный математической моделью и соответствующий температуре ;
a - температуропроводность;
λ - теплопроводность;
m - номер отсчета по координате x, ;
n - номер отсчета по координате y, ;
p - номер отсчета по координате z, ,
x, y, z - оси правой декартовой прямоугольной системы координат;
k - номер отсчета по времени;
h - шаг сетки по расстоянию;
Δτ - шаг сетки по времени;
q - плотность теплового потока;
α - коэффициент теплоотдачи;
K1α, K2α, K3α - линейные коэффициенты зависимости температуропроводности от температуры;
K1λ, K2λ, K3λ - линейные коэффициенты зависимости теплопроводности от температуры;
τ - текущее время, отсчитываемое от момента подачи импульса;
τ* - длительность теплового импульса;
ξ1,ξ2 - параметры регуляризации;
J - функционал невязки.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2008102328/28A RU2374631C2 (ru) | 2008-01-21 | 2008-01-21 | Способ определения комплекса теплофизических свойств твердых материалов |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2008102328/28A RU2374631C2 (ru) | 2008-01-21 | 2008-01-21 | Способ определения комплекса теплофизических свойств твердых материалов |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2008102328A RU2008102328A (ru) | 2009-07-27 |
RU2374631C2 true RU2374631C2 (ru) | 2009-11-27 |
Family
ID=41048049
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2008102328/28A RU2374631C2 (ru) | 2008-01-21 | 2008-01-21 | Способ определения комплекса теплофизических свойств твердых материалов |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2374631C2 (ru) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2534429C1 (ru) * | 2013-06-13 | 2014-11-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО "ТГТУ" | Способ измерения теплофизических свойств твердых материалов методом плоского мгновенного источника тепла |
CN105675646A (zh) * | 2016-03-23 | 2016-06-15 | 哈尔滨工业大学 | 基于本征光热信息同时测量高温半透明介质热导率及吸收系数的方法 |
RU2601234C1 (ru) * | 2015-07-13 | 2016-10-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО "ТГТУ" | Способ измерения теплофизических свойств теплоизоляционных материалов методом плоского импульсного источника теплоты |
RU2611080C1 (ru) * | 2016-01-22 | 2017-02-21 | Роман Владимирович Мироненко | Установка по определению критического значения лучистого теплового потока для различных материалов и веществ |
RU2613194C1 (ru) * | 2015-11-02 | 2017-03-15 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тамбовский государственный технический университет" (ФГБОУ ВО "ТГТУ") | Способ измерения теплофизических свойств анизотропных материалов методом линейного импульсного источника теплоты |
-
2008
- 2008-01-21 RU RU2008102328/28A patent/RU2374631C2/ru not_active IP Right Cessation
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2534429C1 (ru) * | 2013-06-13 | 2014-11-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО "ТГТУ" | Способ измерения теплофизических свойств твердых материалов методом плоского мгновенного источника тепла |
RU2601234C1 (ru) * | 2015-07-13 | 2016-10-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО "ТГТУ" | Способ измерения теплофизических свойств теплоизоляционных материалов методом плоского импульсного источника теплоты |
RU2613194C1 (ru) * | 2015-11-02 | 2017-03-15 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тамбовский государственный технический университет" (ФГБОУ ВО "ТГТУ") | Способ измерения теплофизических свойств анизотропных материалов методом линейного импульсного источника теплоты |
RU2611080C1 (ru) * | 2016-01-22 | 2017-02-21 | Роман Владимирович Мироненко | Установка по определению критического значения лучистого теплового потока для различных материалов и веществ |
CN105675646A (zh) * | 2016-03-23 | 2016-06-15 | 哈尔滨工业大学 | 基于本征光热信息同时测量高温半透明介质热导率及吸收系数的方法 |
CN105675646B (zh) * | 2016-03-23 | 2018-04-24 | 哈尔滨工业大学 | 基于本征光热信息同时测量高温半透明介质热导率及吸收系数的方法 |
Also Published As
Publication number | Publication date |
---|---|
RU2008102328A (ru) | 2009-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2374631C2 (ru) | Способ определения комплекса теплофизических свойств твердых материалов | |
RU2701775C1 (ru) | Способ определения кинетических теплофизических свойств твердых материалов | |
Peeters et al. | Determining directional emissivity: Numerical estimation and experimental validation by using infrared thermography | |
JP3568271B2 (ja) | レーザフラッシュ法を用いた熱定数の測定方法及びその装置 | |
Kirsanov et al. | Measurement of thermal relaxation and temperature damping time in a solid | |
RU2534429C1 (ru) | Способ измерения теплофизических свойств твердых материалов методом плоского мгновенного источника тепла | |
CN106546353B (zh) | 一种均匀材质构件内部温度场的超声测量方法 | |
RU2568983C1 (ru) | Способ определения коэффициента теплопроводности жидкой тепловой изоляции в лабораторных условиях | |
RU2284030C2 (ru) | Способ определения комплекса теплофизических свойств твердых материалов | |
RU2610115C1 (ru) | Устройство для определения температуры газа в полых высокотемпературных элементах газотурбинных двигателей | |
RU2324165C1 (ru) | Способ идентификации комплекса теплофизических свойств твердых материалов | |
RU2613194C1 (ru) | Способ измерения теплофизических свойств анизотропных материалов методом линейного импульсного источника теплоты | |
RU2328724C1 (ru) | Способ идентификации комплекса теплофизических свойств твердых материалов | |
RU2701881C1 (ru) | Устройство для бесконтактного определения теплофизических свойств твердых тел | |
RU2478939C1 (ru) | Способ измерения коэффициента температуропроводности теплоизоляционных материалов методом регулярного режима третьего рода | |
RU2303777C2 (ru) | Способ идентификации комплекса теплофизических свойств твердых материалов | |
RU2324164C1 (ru) | Способ идентификации комплекса теплофизических свойств твердых материалов | |
RU2287807C1 (ru) | Способ определения теплофизических свойств многослойных строительных конструкций и изделий | |
RU2686859C1 (ru) | Способ измерения теплового сопротивления между корпусом полупроводникового прибора и радиатором охлаждения | |
RU2625599C1 (ru) | Способ определения теплопроводности твердых тел | |
Nenarokomov et al. | Investigation of radiative-conductive heat transfer by noncontact measurements. Identification of mathematical models of heat transfer for nontransparent materials | |
RU2255329C1 (ru) | Способ определения теплофизических свойств материалов | |
Chudzik | Measuring system with a dual needle probe for testing the parameters of heat-insulating materials | |
RU2324166C1 (ru) | Способ идентификации комплекса теплофизических свойств твердых материалов | |
RU2327148C1 (ru) | Способ неразрушающего контроля теплофизических свойств многослойных строительных конструкций |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20200122 |