Погрешности и условия применения импульсных методов определения теплофизических характеристик материалов тема автореферата и диссертации по физике, 01.04.14 ВАК РФ
Кац, Марк Давыдович
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Томск
МЕСТО ЗАЩИТЫ
|
||||
2009
ГОД ЗАЩИТЫ
|
|
01.04.14
КОД ВАК РФ
|
||
|
На правах, рукописи
я
Кац Марк Давидович
ПОГРЕШНОСТИ И УСЛОВИЯ ПРИМЕНЕНИЯ ИМПУЛЬСНЫХ МЕТОДОВ ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК МАТЕРИАЛОВ
Специальность 01.04.14 - теплофизика и теоретическая теплотехника
Автореферат диссертации на соискание ученой степени кандидата физико-математических наук
Томск -2009
003476460
Работа выполнена на кафедре автоматизации теплоэнергетических процессов ГОУ-ВПО «Томский политехнический университет»
Научный руководитель: доктор физико-математических наук,
профессор
Кузнецов Гений Владимирович
Официальные оппоненты: доктор физико-математических наук,
профессор
Борзых Владимир Эрнестович
доктор физико-математических наук, профессор
Крайнов Алексей Юрьевич
Ведущая организация: Московский институт теплотехники
Защита состоится Л2 Сл ч TJ&jpSL 2009 года в 14 часов 30 минут на заседании Совета по защите докторских и кандидатских диссертаций ДС 212.025.01 при Томском политехническом университете по адресу: 634050, г. Томск, пр. Ленина, 30, ауд. 228 (10 корпус ТПУ)
С диссертацией можно ознакомиться в Научно-технической библиотеке Томского политехнического университета
Автореферат разослан 2009 года.
Ученый секретарь совета по защите докторских и кандидатских диссерт;
О.Ю. Долматов
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы.
В настоящее время ежегодно создаются новые конструкционные материалы, важнейшими характеристиками которых являются теплофизические (ТФХ): удельная теплоемкость с, коэффициенты теплопроводности Я и температуропроводности а. В связи с растущими объемами испытательных, исследовательских и производственных задач в энергетике, машиностроении и строительстве постоянно разрабатываются новые и совершенствуются известные высокотехнологичные методы определения ТФХ.
В последние десятилетия большой интерес вызывает метод лазерного импульса (Parkers method, Flash method, метод лазерного импульсного нагрева), сущность которого состоит в поглощении в тонком слое нагреваемой («горячей») поверхности образца импульса энергии и регистрации изменения во времени температуры его обратной («холодной») поверхности. Полученная по итогам экспериментов информация позволяет рассчитать теплофизические характеристики исследуемого материала.
Метод лазерного импульсного нагрева постоянно совершенствуется, а также расширяется перечень материалов, для которых с его использованием определяются ТФХ (металлы и сплавы, изоляционные, биологические, анизотропные, пленочные, частично прозрачные материалы, а также жидкие металлы, жидкости и расплавы солей, двухслойные системы). В японском институте метрологии изучается возможность применения этого метода для эталонного измерения коэффициента температуропроводности.
Однако практическая реализация метода лазерного импульса не соответствует его математической модели, в основе которой лежит приближенное решение нестационарной одномерной задачи теплопроводности для пластины при бесконечно малой продолжительности воздействия на нее лазерного луча с однородной энергетической плотностью и при отсутствии теплообмена с окружающей средой. Анализ литературы по рассматриваемой проблеме показал, что математическая постановка задачи, лежащая в основе метода лазерного импульсного нагрева, является источником методических погрешностей, которые необходимо учитывать в реальном эксперименте.
Основные погрешности определения ТФХ материалов импульсным методом лазерного воздействия обусловлены явными и неявными предположениями, сделанными при постановке задачи. Это ограничения по условиям теплообмена на границах; допущения о неограниченной величине теплового потока лазерного излучения к нагреваемой поверхности, бесконечно малой длительности теплового импульса, пространственной однородности теплового потока; одномерности процесса теплопереноса в образце.
Влияние длительности импульса и теплоотвода с «холодной» и «горячей» поверхностей образца исследовались и ранее, но при этом использовались, в основном, одномерные модели, в которых не учитывались неодномерность процесса теплопереноса в образце, локальность импульсного воздействия и j другие значимые факторы. Математическое моделирование температурных /
полей в образце в условиях реальных экспериментов может стать основным методом анализа как погрешности определения ТФХ методом импульсного лазерного нагрева, так и выделения реального диапазона изменения параметров эксперимента, при которых возможно достоверное определение ТФХ. Соответственно, анализ погрешностей этого метода с применением двухмерной и трехмерной моделей нестационарного процесса теплопроводности в образце материала при проведении эксперимента является актуальной задачей, имеющей как фундаментальное, так и прикладное значение.
Цель диссертационной работы.
Целью диссертационной работы является теоретический анализ погрешностей и условий применения импульсных методов определения ТФХ материалов при помощи математического моделирования нестационарных процессов теплопроводности в образцах материалов.
Для достижения поставленной цели необходимо решить следующие задачи:
1. Анализ влияния конвективного и лучистого теплообмена на нагреваемой и «холодной» поверхностях образца на результаты определения теплофизических характеристик.
2. Исследование влияния типа материала на достоверность определения ТФХ.
3. Оценка влияния длительности импульса на результаты определения теплофизических характеристик.
4. Анализ влияния локальности воздействия в условиях энергетической однородности лазерного луча на результаты определения ТФХ.
5. Исследование влияния локальности воздействия и энергетической неоднородности лазерного луча на результаты определения теплофизических характеристик.
6. Анализ влияния неодномерности процесса теплопереноса в материале на результаты определения ТФХ.
Научная новизна работы. Научная новизна диссертационной работы заключается в следующем:
- впервые сформулирован подход к анализу погрешностей и условий применения импульсных методов определения ТФХ материалов, в основе которого лежит численное моделирование реальных температурных полей, формирующихся в образцах в процессе экспериментов, из решения нестационарной задачи теплопроводности с нелинейными граничными условиями при высоких градиентах температур в приповерхностном слое образца, высоких темпах нагрева, локальности теплового воздействия и сравнение вычисленных по результатам численных экспериментов значений теплофизических характеристик материалов с истинными;
- численно исследовано влияние основных параметров (плотности теплового потока, времени воздействия лазерного луча, локальности нагрева, условий теплообмена с окружающей средой) на погрешности экспериментального определения ТФХ импульсным методом лазерного нагрева;
- установлены границы диапазонов возможного изменения величины импульсного теплового потока к поверхности материала, продолжительности
импульса, геометрических размеров образцов, обеспечивающих достоверность полученных результатов в зависимости от типа материала.
Практическая значимость работы.
Внесен вклад в развитие представлений о режимах проведения экспериментов по определению ТФХ методом лазерного импульсного нагрева. Проведен теоретический анализ погрешностей по определению теплофизических характеристик импульсными методами. Результаты диссертационной работы являются теоретической основой для определения параметров и границ применения импульсных методов.
Защищаемые положения:
1. Новый подход к анализу погрешностей и условий применения импульсных методов определения ТФХ материалов, в основе которого лежит численное моделирование реальных температурных полей, формирующихся в образцах в процессе экспериментов, из решения нестационарной задачи теплопроводности с нелинейными граничными условиями при высоких градиентах температур в приповерхностном слое образца, высоких темпах нагрева, локальности теплового воздействия и сравнение вычисленных по результатам численных экспериментов значений теплофизических характеристик материалов с истинными. Этот подход на этапе выбора методики экспериментального определения теплофизических характеристик материалов позволяет оценить как возможность применения импульсного метода для конкретного материала, так и величины методических погрешностей, характерные для этого метода и соответствующего материала.
2. Результаты численного анализа погрешностей метода лазерного импульсного нагрева в одномерной постановке, которые позволяют сделать обоснованное заключение о существенных ограничениях по толщинам образцов и видам материалов.
3. Результаты численного анализа погрешностей метода лазерного импульсного нагрева в двухмерной постановке, позволяющие сделать выводы о масштабах погрешностей в определении ТФХ материалов, обусловленных неодномерностью процесса реального теплопереноса в образце при проведении эксперимента, локальностью лазерного воздействия и неоднородностью плотности теплового потока лазерного импульса.
4. Результаты численного анализа погрешностей метода лазерного импульсного нагрева в трехмерной постановке, на основании которых можно сделать обоснованное заключение о том, что методические погрешности реально оцениваются только при решении нестационарной пространственной задачи теплопроводности и составляют в реальных условиях по коэффициентам теплопроводности и удельной теплоемкости от 4 % до 180 % от численного значения определяемых величин (в зависимости от отношения площадей зоны лазерного нагрева и поверхности образца).
Достоверность полученных результатов.
Обоснованность и достоверность численных исследований следует из проведенных проверок на сходимость и устойчивость решений на множестве сеток, а также подтверждается сравнением результатов, полученных при
решении задачи теплопроводности при импульсном лазерном нагреве аналитическим путем и численно. Полученные в диссертационной работе численные результаты моделирования для простейших постановок соответствуют аналитическому решению задачи импульсного нагрева.
Личный вклад автора.
Все работы по теме диссертации осуществлены автором или при его основном участии: постановка нестационарной задачи теплопроводности с нелинейными граничными условиями при высоких градиентах температур в приповерхностном слое образца, темпах нагрева и локальности теплового воздействия; выбор методов и разработка алгоритма её решения; определение численных значений величины и продолжительности теплового потока; проведение теоретических исследований влияния различных факторов на погрешности определения теплофизических характеристик методом импульсного лазерного нагрева; обработка и анализ полученных результатов, формулировка основных выводов диссертационной работы.
Апробация работы.
Основные результаты работы докладывались и обсуждались на следующих конференциях: 8-й конференции «Теплофизика и гидродинамика технологических процессов» (Томск, 1989), 4-й Всероссийской научно- технической конференции "Энергетика: Экология, Надежность, Безопасность" (Томск , 1998), 14-й Всероссийской научно-технической конференции "Энергетика: Экология, Надежность, Безопасность" (Томск, 2008), «Современные технические системы и технологии» (Севастополь, 2008), «Теплофизические основы энергетических технологий» (Томск, 2009).
Публикации.
Основные результаты, полученные при выполнении диссертации, опубликованы в 15 работах, в числе которых одно авторское свидетельство на изобретение, одно свидетельство о государственной регистрации программы для ЭВМ, пять статей в журналах, рекомендованных ВАК РФ, четыре статьи депонированы ВИНИТИ. Также опубликованы материалы в сборниках трудов трех Всероссийских конференций.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы диссертационной работы, сформулирована цель работы, показана новизна и практическая значимость полученных результатов, представлены положения, выносимые на защиту.
В первой главе проведен анализ современного состояния проблемы погрешностей и условий применения импульсных методов определения ТФХ материалов. Показано, что возможности аналитических или приближенных решений не позволяют учитывать такие важнейшие факторы, как нелинейные эффекты, связанные с теплоотводом с нагреваемой и «холодной» поверхностей образца материала за счет излучения, локальность теплового воздействия лазерного импульса, неоднородность процессов теплопроводности в образце.
Отмечено, что во всех публикациях по применению импульсных методов определения ТФХ материалов в качестве примеров рассматриваются задачи для
материалов с высокой теплопроводностью (металлы и сплавы). В реальной же практике наибольший интерес представляют теплофизические характеристики новых неметаллических материалов, для которых существуют жесткие ограничения по максимально допустимым температурам разогрева.
Во второй главе приведена методика исследований по сформулированной проблеме, сущность которой заключается в решении комплекса задач теплопроводности, соответствующих реальным условиям экспериментов (неоднородность теплопереноса в образце, локальность воздействия лазерного луча, ограниченность мощности лазерного воздействия, распределение энергии излучения по поперечной координате и др.) по определению ТФХ материалов импульсными лазерными методами. Учет факторов, соответствующих реальным условиям экспериментов, исключает возможность применения аналитических или приближенных методов решения задач теплопроводности. Поэтому все исследования проводились методом численного моделирования процессов теплопроводности в образцах материалов. Решение всех задач выполнено методом конечных разностей.
Анализ влияния группы факторов на погрешности определения теплофизических характеристик материалов проводился на основе одномерной, двухмерной и трехмерной моделей теплопроводности для области решения с неоднородными и нестационарными граничными условиями.
Представлены физическая и Математическая постановка задачи при лазерном нагреве поверхности образцу материала тепловым потоком конечной величины и продолжительности импульса. Рассматривается нестационарная
нелинейная задача теплопроводности материала при локальном воздействии на его поверхность теплового импульса конечной величины и временной протяженности. Температура исследуемого материала в начальный момент времени постоянна и равна Т0. Участок «горячей» поверхности (х = 0) нагревается импульсом энергии теплового потока д (рис. 1). При этом происходит нагрев приповерхностного слоя материала с дальнейшим процессом теплопереноса внутри образца. Учитывается теплообмен с окружающей средой за счет конвекции и излучения. Вклад радиационной теплопередачи с поверхностей материала определен законом Стефана-Больцмана, что приводит к рассмотрению нелинейных граничных условий. При постановке задачи приняты следующие допущения:
1) материал образца однородный и гомогенный;
2) теплофизические параметры материала не зависят от температуры. Учитывались основные физические процессы, протекающие при воздействии
ь
0 ^ ,2
Рис.1. Схема области решения: Ь — толщина, //-ширина, Z — высота образца; (у2,у1 - г2,г1) - размеры зоны лазерного нагрева; плотность теплового потока.
У< У О н У
9г
9х
импульсного лазерного излучения на материалы: теплообмен с окружающей средой, локальность и пространственная неоднородность лазерного воздействия, продолжительность импульса нагрева, неодномерность процесса теплопереноса в материале. Соответственно решения задач рассматривались в одномерной, двухмерной и трехмерной нестационарных постановках.
Проведен анализ погрешностей определения ТФХ материалов импульсным методом по результатам решения одномерной задачи теплопроводности (рис. 2) при воздействии теплового импульса на О _ поверхность (х=0) бесконечной пластины при
двух вариантах граничных условий: 1) отсутствие теплообмена с окружающей средой; 2) теплообмен за счет конвекции и излучения с окружающей средой. Особенностями решаемой задачи являются
Чг
А" А-
Чх
Рис. 2. Схема области решения.
Рис Л. Схема области решения.
конечные значения теплового потока нелинейность граничных условий.
Анализ влияния трех ранее не рассматривавшихся факторов (неодномерность процесса теплопереноса в материале, обусловленная локальным воздействием лазерного луча постоянной энергетической плотности; влияние теплообмена с поверхностей образца при локальном воздействии лазерного луча постоянной энергетической плотности; учет воздействия лазерного луча переменной по его сечению энергетической плотности) на погрешности определения ТФХ материалов импульсным методом сводится к решению двухмерной задачи теплопроводности при воздействии теплового импульса на поверхность (х=0) пластины. Область решения задачи представляет пластину толщиной Ь и шириной Н (рис. 3), участок (*=0) которой нагревается мощным импульсом энергии ц на участке у1...у2. Математическая постановка задачи, учитывающая теплообмен с поверхностей образца конвекцией и излучением, включает в себя двухмерное нестационарное уравнение теплопроводности с соответствующими начальным и граничными условиями.
При постановке задачи в трехмерной постановке учитывались не рассматривавшиеся ранее факторы: 1) неодномерность процесса теплопереноса в образце конечных размеров, обусловленная локальным воздействием лазерного луча постоянной энергетической плотности при отсутствии теплообмена с окружающей средой; 2) влияние теплообмена с поверхностей образца при локальном воздействии лазерного луча постоянной энергетической плотности; 3) учет воздействия лазерного луча переменной по его сечению
энергетической плотности. Область решения представляет пластину в форме прямоугольного параллелепипеда толщиной Ь, шириной Н и высотой Участок поверхности х=0 нагревается импульсом энергии я (рис. 1). Задача сводится к решению трехмерного нестационарного уравнения теплопроводности (1) с граничными (2)-{9) и начальным (10) условиями при учете конвективного теплообмена с поверхностей образца и интенсивности импульса энергии д(у,г)\
х = 0:
дТ , ср— = Я и 81
[дх
82Т &Т дх2 + &
2 &2/
= 4(у,г) + аг(Те-Т),
1>0,0<х<£, 0 <у<Н,
У^У^Уи г, < г < г2,
(1)
(2)
х = 0: -Л
х = Ь\
х = 0:
дГ дх
дТ_ дх
К®.
= ат{Те-Т),
= аг(Те-Т),
= ах(Те-Т),
У<У1-,У>У2,
0<у<Н, 0 <г<1¥,
{>гкип;
0<у<Н, 0<г<}У, ?>0;
(3)
(4)
(5)
>> = 0:
6Г
= 0,
0<х<1, 0<г<1Г, (>0;
(6)
У = Н:
ду
0<х<Ь, 0<г<Ж, 1>0;
(7)
0 <Х <Ь,
0 <у<Н, />0;
(8)
/ ^у1 ^
0<х<Ь, 0<у<Н; ' > 0,;
(9)
Г = 0: Г = Г0,
0 < л: 2 0<у<,Н,
(10)
Разностные аналоги исходных дифференциальных уравнений решены методом конечных разностей с применением неявной четырехточечной разностной схемы. Одномерные нелинейные разностные уравнения решены методом итераций с применением метода прогонки на каждом шаге итераций.
Особенность рассматриваемой задачи заключается в том, что в реальном возможном диапазоне изменения условий воздействия лазерного излучения на материалы в тонком приповерхностном слое образца возникают очень большие градиенты температур. При численном моделировании такого процесса возникает необходимость верификации разрабатываемых методик и алгоритмов численного решения задачи. С этой целью проводились специальные численные комплексные исследования для выбора сеточных параметров, исходя из условия баланса энергии в области решения. Анализ результатов показывает, что относительная погрешность вычисления поля температуры составляет не более 0,03 %. Поэтому имеется основание для заключения о корректности выбранного метода численного расчета.
В третьей главе представлены результаты численного исследования процессов теплопроводности в материале и анализ погрешностей, соответствующих реальным условиям проведения экспериментов по определению ТФХ материалов импульсными лазерными методами.
В качестве материалов исследуемых образцов рассматривались медь, сталь, карбид циркония и полимерные материалы. Выбор таких материалов не
накладывает каких-либо существенных ограничений на общность полученных результатов. Для установления масштабов влияния конвективного и лучистого теплообмена с нагреваемой и «холодной» поверхностей материала на результаты определения теплофизических характеристик проведено численное решение одномерной задачи теплопроводности при конечной величине интенсивности импульса лазерного нагрева.
На рис. 4 приведены типичные результаты численного моделирования для стальной пластины в виде распределения Т(х) для момента времени 1=0,3 с с учетом конвективного теплоотвода на «горячей» и «холодной» поверхностях при температурах Тп = 293 К и Те= 273 К . Из рис. 4 видно, что с увеличением коэффициента конвективного теплообмена, как и следовало ожидать, наблюдается уменьшение значения максимальной температуры «холодной» поверхности образца.
На рис. 5 приведена зависимость погрешности определения удельной
теплоемкости 8 от величины коэффициента конвективного теплообмена. Рис. 5 показывает, что с увеличением интенсивности естественной конвекции
Л
VI
VI
307,6 307.58 307,56 307,54 307,52 307.5 307,48
0 (и 0,1 0,6 0.В Д-Ю.М
Рис. 4. Распределение Т(х), 1=0,3 с, при а: 1)0; 2)10; 3)25 Вт/(м2К).
4 8 12 Рис. 5. Зависимость погрешности б определения удельной теплоемкости от коэффициента конвективного теплообмена а
величина погрешности возрастает и достигает 1,1 % при значении коэффициента конвективного теплообмена 20 Вт/(м2К).
Решение одномерной задачи теплопереноса проведено также на примере карбида циркония. Распространение тепла в карбиде циркония аналогично процессу теплопереноса стали. Влияние конвективного и лучистого теплообмена при малых величинах температур материала и среды приводят к незначительному росту погрешности определения удельной теплоемкости и температуропроводности. В реальных условиях при локальном разогреве поверхности до 1000 К и более при использовании более мощного источника нагрева интенсивность естественной конвекции и излучения в окрестности зоны нагрева может быть существенно выше, что приведет к увеличению погрешности определения теплофизических характеристик материалов.
В работе обоснованы условия применения метода лазерного импульса для многочисленных материалов с относительно низкой теплопроводностью (Я =0,2...0,5Вт/(м-К)) в зависимости от величины и продолжительности теплового импульса, толщины образца, типа материала, которые ранее отсутствовали. Для установления влияния величины импульсного теплового потока к поверхности материала и толщины экспериментальных образцов (в зависимости от типа материала), обеспечивающих достоверность определения теплофизических характеристик импульсным методом, была численно решена одномерная задача теплопроводности при конечной величине интенсивности импульса лазерного нагрева. Численный анализ температурных полей проводился для типичных полимерных материалов (полиэтилен, полиамид, полибутилентерефталат).
Распределения 7Хх) для образцов полимерных материалов при длительности импульса 0,1 с и величине теплового потока д = 108 Вт/м2 показали наличие жестких ограничений по величине теплового потока, длительности импульса, толщинам образцов. Не для всех материалов возможна реализация импульсных методов определения ТФХ вследствие превышения температуры нагреваемой поверхности образца выше температуры начала деструкции материала. Установлена необходимость существенного уменьшения величин предельных (максимальных) тепловых потоков для определения ТФХ материалов с низкой теплопроводностью импульсным методом. При толщине образца 5 мм величина максимального теплового потока составила: для полиэтилена высокой плотности 2,76 105 Вт/м2; полиамида ПА-610- 3,42Т05 Вт/м2; полибутилентерефталата - 5,06-Ю5 Вт/м2.
Численные исследования показали необходимость существенного уменьшения толщин образцов для определения ТФХ материалов с низкой теплопроводностью импульсным методом. Максимальные значения толщин образцов, при которых температура «холодной» стенки возрастает не менее чем на необходимое для расчетов ТФХ значение Т^ =3...6К(при предельных величинах теплового потока), составляют: для полиэтилена высокой плотности 1,8 мм; полиамида ПА-610 - 2,0 мм; для полибутилентерефталата -1 мм. Полученные результаты показывают, что при определении ТФХ материалов с
низкои теплопроводностью импульсными методами существуют ограничения по величине теплового потока и толщине образцов, что накладывает ограничения на типы исследуемых материалов, для которых возможно изготовление образцов толщиной 1-2 мм без изменения их характерной структуры.
Для установления влияния продолжительности импульса гям„ на результаты определения теплофизических характеристик проведено численное решение одномерной задачи теплопроводности при конечной величине теплового потока лазерного нагрева с учетом радиационного теплоотвода на «горячей» и «холодной» поверхности образца.
На рис. 6 приведены типичные результаты численного моделирования поля температур стали 1Х18Н9Т. Установлено, что, как и следовало ожидать, с увеличением длительности импульса растет температура на «горячей» поверхности образца. Проведено определение характеристик Ттвх и т0>5 при условии постоянства величины теплового потока q и изменения длительности импульса гамп = 0,1...12 мс. Установлено, что при увеличении продолжительности
Рис. 6. Распределение Т(ж), (=20 мс, при : 1) 0.1; 2) 5,0 ; 3) 10,0 мс.
импульса значения
Т
шах
и г05 растут с
/
у и
/
/
также аналитическими расчетами. Расчет удельной теплоемкости на примере стали 1Х18Н9Т показывает только
незначительное возрастание погрешности (до 0,2 %) с увеличением
продолжительности импульса, т. к. с увеличением ткиа возрастает значение поглощенной образцом теплоты £>. На рис. 7 приведена зависимость погрешности определения коэффициента
температуропроводности стали 1Х18Н9Т 3 от величины длительности импульса гимп. С увеличением гимп погрешность определения а возрастает существенно и достигает 14 % при значении времени импульса гимп = 0,012с. Аналогичные результаты получены и для других материалов. Проведенные численные исследования позволяют сделать вывод, что с ростом длительности импульса существенно увеличиваются погрешности определения коэффициентов температуропроводности и теплопроводности материалов.
О 0.004 0.008 0.012 г„,
Рнс.7. Зависимость погрешности определения коэффициента температуропроводности д от продолжительности импульса
Для анализа влияния локальности воздействия в условиях энергетической однородности лазерного луча проведены численные исследования процесса теплопроводности в двумерной постановке. На рис 8. показаны типичные распределения температуры Т{х,у) в образце стали Ст. 10 в моменты времени ( = 0,02 с и 1 = 0,06 с при гямв =0,00} 5 с.
»т. к .............................у................_ т. К ..........................г..............
Рис. 8. Распределение Т(х,у) в момент времени I: а) 0,02; б) 0,06 с.
При воздействии лазерного импульса на границе х = 0 распространение тепла осуществляется не только по координатному направлению х, но и по координатному направлению у. При этом температурный профиль имеет ярко выраженный максимум в центре действия лазерного луча. Температура в зоне лазерного нагрева достигает 304 К. В направлении от центра зоны лазерного нагрева к внешним граням пластины температура резко снижается (рис. 8а). В момент времени / = 0,06 с (рис. 86) происходит выравнивание температур «горячей» и «холодной» поверхностей пластины в зоне лазерного нагрева с продолжением перетока тепла по поперечной координате у. Полученные результаты позволяют сделать вывод о существенно неодномерном процессе распространения тепла в исследуемой системе. Этот вывод может показаться очевидным, если не принимать во внимание, что в подавляющем большинстве работ по моделированию лазерного импульсного воздействия на материалы рассматривались одномерные модели.
Типичными являются результаты численного анализа погрешностей определения ТФХ стали Сг. 10. При сравнении результатов численного моделирования с аналитическими результатами решения одномерной задачи импульсного лазерного нагрева установлено, что погрешность определения температуры Ттах составляет около 0,3 %, а погрешность определения времени Т05 — до 11 %, что сказывается при расчете ТФХ. Так погрешность определения удельной теплоемкости стали составляет до 23 %, коэффициента температуропроводности до 12 %, коэффициента теплопроводности до 35 %.
Учет конвективного и лучистого теплообмена на поверхностях пластины приводит, как и следовало ожидать, к увеличению погрешности. Погрешность
а
б
у,10н
определения теплоемкости при значении коэффициента конвективного теплообмена в режиме свободной конвекции при а - 20 Вт/(м2 К) и степени черноты поверхностей е = 0,5 составляет до 36 %.
Решение двухмерной задачи теплопереноса проведено также на примерах карбида циркония и полиамида ПА-610. Распространение тепла в карбиде циркония и полиамиде аналогично процессу теплопереноса в стали. Однако процесс теплопереноса в полиамиде происходит с существенно меньшей интенсивностью по сравнению со сталью. Выделенная особенность обусловлена различием коэффициентов теплопроводности полиамида и стали.
Для анализа влияния локальности воздействия и энергетической
неоднородности лазерного луча проведен численный анализ погрешностей определения ТФХ на базе двухмерной модели теплопроводности при конечной величине теплового потока лазерного нагрева.
Величина теплового потока лазерных установок не является постоянной по пространственными координатам ввиду наличия поперечных мод, определяемых условиями резонанса внутри установки. Пространственная неоднородность энергетической
характеристики источника лазерного излучения определяется выражением: = • ехр(—>>2 / Л2), где д0 - максимальная плотность теплового потока в центре лазерного луча, Вт/м ; к— коэффициент сосредоточенности потока.
На рис. 9 показано распределение температуры (положения характерных изотерм) стали Ст.10 в момент времени 1=0,02 с в плоскости оси лазерного луча. Распространение тепла осуществляется не только по координатному направлению оси х, но и по координатному направлению у. При этом температурный профиль имеет ярко выраженный максимум в точке, соответствующей наибольшему значению величины теплового, потока. Температура в центре зоны лазерного нагрева достигает 302 К при д0=2,25-Ю1 Дж/м2 . В направлении к границам пластины температура постепенно уменьшается. Полученные результаты позволяют сделать обоснованный вывод о существенно нестационарном и двухмерном характере процесса распространения тепла в исследуемом материале, что приводит уменьшению значения Ттах на «холодной» поверхности образца по сравнению с одномерной постановкой задачи, не учитывающей неоднородности теплового потока по сечению лазерного луча.
Сравнение результатов проведенного численного моделирования Т(х,у) при значении коэффициента сосредоточенности теплового потока к =0,06 с полученными аналитически по методике Паркера значениями температуры Ттах
Рис. 9, Распределение Т(х,у) стали в плоскости оси лазерного луча.
и времени г05 показало, что погрешность определения температуры составляет около 0,36 %, а отклонение в определении времени до 13,5 %. Последнее соответственно сказывается при вычислении ТФХ. Так погрешность определения удельной теплоемкости рассматриваемой стали составляет до 30 %, коэффициента температуропроводности до 14 %, коэффициента теплопроводности до 44 %. Полученные результаты позволяют сделать обоснованный вывод о существенно нестационарном и двухмерном характере процесса распространения тепла в исследуемом материале, что приводит уменьшению значения 7„ИЛ- на «холодной» поверхности образца по сравнению с постановкой задачи, не учитывающей неоднородности теплового потока по сечению лазерного луча.
Полученные результаты позволяют сделать вывод о том, что в реальной практике определение теплофизических характеристик материалов методом импульсного нагрева сопряжено с очень высокими (и неприемлемыми во многих практически важных приложениях) ошибками. Из приведенных результатов можно также сделать вывод, что двухмерная модель существенно более точно описывает температурное поле материала при воздействии импульсного лазерного источника. При этом масштабы погрешностей настолько велики, что можно сделать вывод о необходимости существенных ограничений размеров образца материала, как по толщине, так и по поперечным направлениям. Очевидно, что методами импульсного лазерного нагрева нецелесообразно определять теплофизические характеристики широко распространенных конструкционных материалов стекло- и углепластиков, а также резиноподобных армированных теплозащитных и теплоизоляционных материалов с наполнением в виде гранул, нитей, капсул, в виду невозможности изготовления образцов без нарушения целостности структуры материала.
Анализ влияния неодномерности процесса теплопереноса в материале
проведён путём численного решения трехмерной задачи теплопроводности при конечной величине теплового потока лазерного нагрева. На рис. 10 показано типичное положение характерных изотерм стали Ст. 10 в плоскости Ъ = 4,5 мм при I = 20 мс. Из рис. 10 следует, что температура изменяется не только по толщине пластины, но и по поперечным координатам, что приводит к существенному снижению значения Тщах на «холодной» поверхности образца по сравнению с одномерной постановкой задачи.
У.Ю"'«
Рис. 10. Положения характерных изотерм стали в плоскости I = 4,5 ММ при / = 20 мс
При сравнении результатов численного моделирования с одномерными аналитическими расчетами по методике Паркера установлено, что погрешность определения температуры Ттах и времени та 5 составляет соответственно около 0,06 % и 0,04 % при нагреве всей «горячей» поверхности образца лазерным импульсом, что сказывается при расчете ТФХ. Так погрешность определения удельной теплоемкости составляет около 4 %, а коэффициента температуропроводности 0,1%. При уменьшении площади зоны лазерного нагрева 8Н и неизменной «горячей» поверхности Б,, погрешности определения ТФХ возрастают. На рис. 11 приведена зависимость погрешности определения
удельной теплоемкости стали 6 от величины отношения площадей 5г/Би. Хорошо видно, что с уменьшением площади зоны лазерного нагрева Бн величина погрешности определения удельной теплоемкости достигает 140 % при 8г/Би= 9. Погрешность определения коэффициента теплопроводности
превышает погрешность определения удельной теплоемкости.
выводы
1.. Впервые сформулирован подход к анализу погрешностей и условий применения импульсных методов определения ТФХ материалов, в основе которого лежит численное моделирование реальных температурных полей, формирующихся в образцах в процессе экспериментов, из решения нестационарной задачи теплопроводности с нелинейными граничными условиями при высоких градиентах температур в приповерхностном слое образца, высоких темпах нагрева, локальности теплового воздействия и сравнение вычисленных по результатам численных экспериментов значений теплофизических характеристик материалов с истинными.
2. Проведен анализ влияния группы факторов на погрешности определения ТФХ материалов на основе одномерной, двухмерной и трехмерной моделей теплопроводности для области решения с неоднородными и нестационарными граничными условиями при высоких градиентах температур в приповерхностном слое образца и темпах нагрева, а также локальности теплового воздействия.
3. Результаты, полученные при решении одномерной задачи теплопроводности в условиях импульсного нагрева, позволяют сделать вывод об увеличении погрешности определения удельной теплоемкости с учетом теплообмена на всех открытых границах образца. Для материалов с низкой теплопроводностью при решении одномерной задачи установлены границы диапазонов возможного изменения интенсивности импульсного теплового потока, продолжительности импульса и толщины образца (в зависимости от типа материала), обеспечивающих достоверность определения ТФХ.
0 2 4 6 8 Эг/Бн Рис. 11. Зависимость погрешности определения удельной теплоемкости 5 отв^..
4. При решении двухмерной задачи теплопроводности импульсного нагрева показано влияние на погрешность определения ТФХ локального воздействия лазерного импульса и неоднородности распределения энергии по поперечной координате. Установлено, что неучет этих факторов приводит к погрешностям определения удельной теплоемкости стали до 25 %, коэффициента температуропроводности до 13%, коэффициента теплопроводности - до 38 %. Показано, что неоднородность пространственной энергетической плотности лазерного потока на поверхности пластины приводит к увеличению погрешностей. Так погрешность определения теплоемкости при значении коэффициента сосредоточенности теплового потока к = 0,06 составляет до 30 %, коэффициента температуропроводности до 14 %, коэффициента теплопроводности до 44 %.
5. Результаты, полученные при решении трехмерной задачи теплопроводности в условиях импульсного нагрева, позволяют сделать вывод о влиянии на погрешность определения ТФХ неодномерного теплопереноса в образце при локальности воздействия лазерного луча. Установлено, что на методику определения ТФХ импульсным методом накладываются жесткие ограничения по размерам образца. Площадь нагреваемой поверхности образца должна соответствовать площади лазерного нагрева.
6. Результаты проведенных исследований показали, что теоретическая оценка погрешностей метода лазерного импульса создает условия для достоверного определения ТФХ любых материалов с установлением методических ошибок, которые можно называть систематическими. Кроме того, знание масштабов этих ошибок создает предпосылки для постановки новых задач по снижению погрешностей определения ТФХ этим методом.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИОННОЙ РАБОТЫ ПРЕДСТАВЛЕНЫ В ПУБЛИКАЦИЯХ
1. Загромов Ю.А., Кац М.Д. Расчет теплофизических характеристик покрытия, нанесенного на металлическое основание, с использованием импульсного метода. ТПУ. - Томск. - 17с. Деп. в ВИНИТИ 09.03.87, №1686-В87.
2. Андык B.C., Загромов Ю.А., Кац М.Д. Использование импульсного метода для измерения теплоемкости материалов. ТПУ. - Томск. — 8 с. Деп. в ВИНИТИ 12.11.87, №17999-В87.
3. A.c. 1465751СССР: МКИ3 G 01 N 25/18. Устройство для определения коэффициента температуропроводности твердых материалов: / В.П. Вавилов, А.И. Иванов, В.В. Ширяев, А.И. Суппес, М.Д. Кац- заяв. 1.06.87; опубл. 15.03.89, Бюл. №10 - Зс.
4. Загромов Ю. А., Кац М. Д. Усилитель постоянного напряжения для регистрации малых изменений температуры //Материалы Всесоюзной конференции «Теплофизика и гидродинамика технологических процессов». -Томск: Изд-во ТПУ. - 1989-С. 128-131.
5. Кац М.Д., Карауш С.А., Загромов Ю.А., Бугаев И.В. Об учете теплообмена с поверхности образца при измерении теплофизических характеристик покрытий импульсным методом. // Материалы Всесоюзной
конференции «Вопросы теплообмена в строительстве». - Ростов н/Д: Ростовский инж. строит, ин-т. - 1990. -С. 136- 141.
6. Кац М.Д., Карауш С.А., Бугаев И.В. Влияние теплопотерь с поверхности двухслойного образца на измерение теплофизических характеристик импульсным методом // Инженерно-физический журнал. - 1991. - Т. 60. - №1. -С. 127-130.
7. Кац М.Д. Расчет температурных полей покрытия на металлическом основании при наличии теплообмена с окружающей средой // Материалы четвертой Всероссийской научно-технической конференции «Энергетика: экология, надежность, безопасность». - Томск: Изд-во ТПУ, 1998. - С. 96-97.
8. Кузнецов Г.В., Кац М.Д. Об условиях применения импульсных методов определения теплофизических характеристик конструкционных материалов // Известия Томского политехнического университета. - 2008. - Т. 312. - №4. -С. 10-13.
9. Кац М.Д. О погрешностях определения теплофизических характеристик конструкционных материалов импульсным методом // Материалы четырнадцатой Всероссийской научно-технической конференции «Энергетика: экология, надежность, безопасность». — Томск: Изд-во ТПУ, 2008.-С. 79-81.
10. Кузнецов Г.В., Кац М.Д. Анализ погрешностей определения импульсными методами теплофизических характеристик конструкционных материалов // Известия Томского политехнического университета. - 2008. - Т. 313.-№4.-С. 5-8.
11. Кузнецов Г.В., Кац М.Д. Теоретический анализ методических погрешностей определения теплофизических характеристик конструкционных материалов импульсным методом в образце конечных размеров // Измерительная техника. -2009. - №4. - С. 35-37.
12. Кузнецов Г.В., Кац М.Д. Анализ погрешности определения теплофизических характеристик конструкционных материалов импульсным методом при конечной длительности импульса лазерного луча // Инженерная физика. -2009. -№8. - С. 45-^8.
13. Кузнецов Г.В., Кац М.Д. Расчет нестационарного двухмерного температурного поля пластины при локальном импульсном тепловом воздействии малой временной протяженности с теплоотводом на границах // Свидетельство о государственной регистрации программы для ЭВМ №2009612866.- Бюлл. прогр. №6,2009.
14. Кузнецов Г.В., Кац М.Д. Погрешности определения теплофизических характеристик карбида циркония импульсным методом, обусловленные локальным воздействием лазерного луча. ТПУ. - Томск. - 12с. Деп. ВИНИТИ 18.06.09, № 382-В2009.
15. Кузнецов Г.В., Кац М.Д. Оценка погрешности определения теплофизических характеристик карбида циркония импульсным методом при пространственной неоднородности лазерного луча. ТПУ. - Томск. - 10с. Деп. ВИНИТИ 01.07.09, № 419 - В2009.
Подписано к печати 18.08.2009. Тираж 110 экз. Кол-во стр. 18. Заказ №34-09 Бумага офсетная. Формат А-5. Печать RISO. Отпечатано в типографии ООО «РауШмбх» Лицензия Серия ПД № 12-0092 от 03.05.2001г. 634034, г. Томск, ул. Усова 7, ком. 046 тел. (3822) 56-44-54
ВВЕДЕНИЕ.
1. СОВРЕМЕННОЕ СОСТОЯНИЕ ПРОБЛЕМЫ ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК МАТЕРИАЛОВ ИМПУЛЬСНЫМИ МЕТОДАМИ.
1.1. Усовершенствование техники эксперимента определения теплофизических характеристик материалов методом импульсного лазерного нагрева.
1.2. Аналитические и численные решения задач теплопереноса в материале при импульсном лазерном нагреве.
2. ПОСТАНОВКА ЗАДАЧИ ТЕПЛОПЕРЕНОСА ПРИ ЛАЗЕРНОМ ИМПУЛЬСНОМ НАГРЕВЕ И МЕТОД ЕЕ РЕШЕНИЯ.
2.1. Физическая постановка задачи.
2.2. Одномерная задача теплопроводности при воздействии импульса лазерного нагрева конечной длительности и интенсивности.
2.3. Метод решения.
2.4. Двухмерная задача теплопроводности при воздействии импульса лазерного нагрева конечной длительности и интенсивности.
2.5. Трехмерная задача теплопроводности при воздействии импульса лазерного нагрева конечной длительности и интенсивности.
2.6. Тестирование математических моделей и метода решения.
3. ЧИСЛЕННЫЙ АНАЛИЗ ПОГРЕШНОСТЕЙ И УСЛОВИЙ ПРИМЕНЕНИЯ- ИМПУЛЬСНЫХ МЕТОДОВ ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК МАТЕРИАЛОВ.
3.1. Анализ влияния конвективного и лучистого теплообмена на результаты определения теплофизических характеристик импульсными методами.
3.2. Влияние типа материала на результаты определения теплофизических характеристик импульсными методами.
3.3. Влияние длительности импульса лазерного луча на результаты определения теплофизических характеристик импульсными методами.
3.4. Анализ влияния локальности воздействия в условиях энергетической однородности лазерного луча на результаты определения теплофизических характеристик импульсными методами.
3.5. Влияние локальности воздействия и энергетической неоднородности лазерного луча на погрешности определения теплофизических характеристик импульсными методами.
3.6. Анализ влияния неодномерности процесса теплопереноса в материале на результаты определения теплофизических характеристик импульсными методами.
3.7. Обобщение полученных результатов. Рекомендации по применению методов импульсного нагрева для определения теплофизических характеристик материалов.
В настоящее время ежегодно создаются новые конструкционные материалы, такие как металлополимерные композиты, наноматериалы, аморфные магнитные сплавы, сплавы для изготовления дисков турбин, вольфраматы кальция* и бария для лазеров, механически легированные композиционные материалы для аэрокосмической техники, материалы специальных защитных покрытий и многие другие, которые обладают определенными преимуществами по сравнению с известными [1-10].
Важнейшими характеристиками любых конструкционных материалов являются теплофизические (ТФХ): удельная теплоемкость с, коэффициенты теплопроводности Л и температуропроводности а [11-24]. В связи с растущими объемами испытательных, исследовательских и производственных задач в энергетике, машиностроении, строительстве для достоверного определения ТФХ постоянно разрабатываются новые и совершенствуются известные методы их определения [17-24].
Существующая в настоящее время общепринятая классификация разделяет методы определения теплофизических характеристик материалов на стационарные и нестационарные [11-17, 20 21]. Эта классификация базируется на решении-уравнения теплопроводности для образца материала определенной геометрической формы с соответствующими граничными и начальными условиями.
Стационарные методы отличаются тем, что распределение температуры в образцах не зависит от времени и сложным образом связано с геометрией- тела, его ТФХ, граничными и начальными условиями. Недостатками стационарных методов являются длительность установления-необходимого теплового режима при заданной температуре и измерение только одного теплофизического параметра [11-17, 20 21].
Нестационарные методы определения ТФХ материалов основаны на временном изменении температурного поля в исследуемом материале и позволяют измерить несколько теплофизических характеристик за один эксперимент с минимальными временными затратами в условиях достаточно высоких температур [12-15. 17-24].
В последние десятилетия большой интерес вызывает метод лазерного импульса (Parkers method, Flash method, метод лазерного импульсного нагрева) [25-55], сущность которого состоит в поглощении в тонком слое нагреваемой «горячей» поверхности образца импульса энергии и регистрации изменения во времени температуры его обратной («холодной») поверхности. Полученная по итогам экспериментов информация позволяет рассчитать теплофизические характеристики исследуемого материала: температуропроводность, теплоемкость и теплопроводность.
В настоящее время продолжается совершенствование метода лазерного импульса и расширяется перечень материалов, для которых с его использованием определяются ТФХ [25-55]. При помощи этого метода находятся а, с, Я металлов и сплавов [30-37], изоляционных [38], биологических [39-40] материалов, жидких металлов [41], жидкостей и расплавов солей [42-43], двухслойных систем [44-47], анизотропных пленочных материалов [48], частично прозрачных материалов [49, 50], материалов ядерных реакторов [51]. В японском институте метрологии изучается возможность применения этого метода для эталонного измерения коэффициента температуропроводности [52-54]. При проведении теплового неразрушающего контроля также используется метод импульсного лазерного нагрева [19. 22, 55]. Все это говорит о неослабевающем интересе исследователей и инженеров к его применению.
Однако практическая реализация метода лазерного импульса не соответствует его математической модели. В основе метода лазерного импульса лежит приближенное решение нестационарной одномерной задачи теплопроводности для пластины при бесконечно малой продолжительности воздействия на нее лазерного луча с однородной энергетической плотностью и при отсутствии теплообмена с окружающей средой [25]. Анализ литературы по рассматриваемой проблеме [26, 27, 35] показал, что математическая постановка задачи, лежащая в основе метода лазерного импульса, является источником методических погрешностей, которые необходимо учитывать в реальном эксперименте. По оценке авторов [33] погрешность определения коэффициента температуропроводности материала методом импульсного лазерного нагрева, как правило, превышает 10 %.
Основные погрешности определения ТФХ материалов импульсным методом лазерного воздействия обусловлены явными и неявными предположениями, сделанными при постановке задачи [25]. Это ограничения по условиям теплообмена на границах; допущения о неограниченной величине теплового потока лазерного излучения к нагреваемой поверхности, бесконечно малой длительности теплового импульса, пространственной однородности теплового потока; одномерности процесса теплопереноса в образце [33].
Влияние отдельных факторов на изменение температурного поля образца исследовались и ранее [25-30, 33, 56-74]. Но в опубликованных к настоящему времени работах по анализу погрешностей определения ТФХ методом лазерного импульса использовались, в основном, одномерные модели, в которых не учитывались неодномерность процесса теплопереноса в образце, локальность импульсного воздействия и другие значимые факторы. Так, например, в [72] получено решение одномерной задачи теплопроводности для лазерного импульса произвольной временной формы при отсутствии теплообмена с поверхностей образца материала. В реальных же условиях процесс неодномерного теплопереноса в образце конечных размеров оказывает значительное влияние на основную значимую характеристику этого метода - температуру «холодной» поверхности образца.
Математическое моделирование температурных полей в образце в условиях реальных экспериментов может стать основным методом анализа как погрешности определения ТФХ методом импульсного лазерного нагрева, так и выделения реального диапазона изменения параметров эксперимента (интенсивности и продолжительности лазерного воздействия, размеров зоны нагрева, геометрических размеров образца, условий теплообмена с окружающей средой), при которых возможно достоверное определение ТФХ. Соответственно анализ погрешностей метода [25] с применением двухмерной и трехмерной моделей нестационарного процесса теплопереноса в образце материала при проведении лазерного импульсного эксперимента является актуальной задачей, имеющей как фундаментальное, так и прикладное значение.
Фундаментальность темы диссертации обусловлена перспективами применения импульсных методов определения ТФХ. Чем более полно будут проанализированы источники погрешностей вычисления Л, с и а при использовании метода лазерного импульсного нагрева и чем более точно будут оценены масштабы этих погрешностей, тем больше появится возможностей у экспериментаторов для создания эффективных методик определения ТФХ на базе метода [25] и анализа новых материалов.
Практическая значимость темы диссертации обусловлена тем, что теоретическая оценка погрешностей метода лазерного импульса создает условия для достоверного определения ТФХ любых материалов с установлением методических ошибок, которые можно называть систематическими. Кроме того, знание масштабов этих ошибок создает предпосылки для постановки новых задач по снижению погрешностей определения ТФХ.
Целью работы является теоретический анализ условий применения и погрешностей импульсных методов определения ТФХ материалов при помощи метода математического моделирования нестационарных процессов теплопроводности в образцах материалов.
Указанная цель достигается решением следующих задач:
1. Анализ влияния конвективного и лучистого теплообмена на результаты определения теплофизических характеристик импульсными методами.
2. Исследование влияния типа материала на результаты определения теплофизических характеристик импульсными методами.
3. Выявление воздействия длительности импульса лазерного луча на результаты определения теплофизических характеристик импульсными методами.
4. Анализ влияния локальности воздействия в условиях энергетической однородности лазерного луча на результаты определения теплофизических характеристик импульсными методами.
5. Исследование локальности воздействия и энергетической неоднородности лазерного луча на результаты определения теплофизических характеристик импульсными методами.
6. Анализ влияния неодномерности процесса теплопереноса в материале на результаты определения теплофизических характеристик импульсными методами.
Научная новизна работы.
1. Впервые сформулирован подход к анализу погрешностей и условий применения импульсных методов определения ТФХ материалов, в основе которого лежит численное моделирование реальных температурных полей, формирующихся в образцах в процессе экспериментов, из решения нестационарной задачи теплопроводности с нелинейными граничными условиями при высоких градиентах температур в приповерхностном слое образца, высоких темпах нагрева, локальности теплового воздействия и сравнение вычисленных по результатам численных экспериментов значений теплофизических характеристик материалов с истинными.
2. Численно исследовано влияние основных параметров (плотности теплового потока, времени воздействия' лазерного луча, локальности нагрева, условий теплообмена с окружающей- средой) на погрешности экспериментального определения, ТФХ импульсным методом лазерного нагрева.
3. Установлены границы диапазонов возможного изменения величины-импульсного теплового потока к поверхности1 материала, продолжительности импульса, геометрических размеров образцов, обеспечивающих достоверность полученных результатов в зависимости от типа материала.
Практическая значимость работы. Внесен вклад в развитие представлений о режимах проведения экспериментов по определению ТФХ методом лазерного импульсного нагрева. Проведен теоретический анализ погрешностей по определению теплофизических характеристик импульсными методами. Результаты диссертационной работы являются теоретической основой для определения- параметров и границ применения импульсных методов.
Достоверность полученных результатов. Обоснованность и достоверность численных исследований-следует из проведенных проверок на сходимость и-устойчивость решений на множестве сеток, подтверждается сравнением результатов, полученных при решении задачи Паркера [25] аналитическим путем и с помощью численного моделирования.
Автор защищает:
1. Новый подход к анализу погрешностей и условий применения импульсных методов определения ТФХ материалов, в основе которого лежит численное моделирование реальных температурных полей, формирующихся в образцах в процессе экспериментов, из решения нестационарной задачи теплопроводности с нелинейными граничными условиями при высоких градиентах температур в приповерхностном слое образца, высоких темпах нагрева и локальности теплового воздействия, и последующее сравнение вычисленных по результатам численных экспериментов значений теплофизических характеристик с истинными. Этот подход на этапе выбора методики экспериментального определения теплофизических характеристик материалов позволяет оценить как возможность применения импульсного метода для конкретного материала, так И'величины методических погрешностей, характерные для этого метода и соответствующего материала.
2. Результаты численного<■ анализа погрешностей метода лазерного импульсного нагрева в одномерной постановке, которые позволяют сделать обоснованное заключение о существенных ограничениях по толщинам образцов и видам материалов.
3. Результаты численного анализа погрешностей метода лазерного импульсного нагрева в двухмерной постановке, позволяющие сделать выводы о» масштабах погрешностей в определении ТФХ материалов, обусловленных неодномерностью процесса реального теплопереноса в образце при проведении эксперимента, локальностью лазерного воздействия и неоднородностью плотности теплового потока лазерного импульса.
4. Результаты численного анализа погрешностей метода лазерного импульсного нагрева в трехмерной постановке, на основании которых можно сделать обоснованное заключение о том, что методические погрешности реально оцениваются только при решении нестационарной пространственной задачи теплопроводности и составляют в реальных условиях по коэффициентам теплопроводности и удельной теплоемкости более от 4 % до 180 % от численного значения определяемых величин.
Аппробация работы: Основные результаты работы докладывались и обсуждались на следующих конференциях: 8-й конференции «Теплофизика и гидродинамика технологических процессов» (Томск, 1989), 4-й Всероссийской научно-технической конференции "Энергетика: Экология, Надежность, Безопасность" (Томск , 1998), 14-й Всероссийской научно-технической конференции "Энергетика: Экология, Надежность, Безопасность" (Томск, 2008), «Современные технические системы, и технологии». (Севастополь, 2008), «Теплофизические основы энергетических технологий» (Томск, 2009).
Публикации. Основные результаты диссертационной работы представлены в трудах вышеперечисленных конференций, а также в журналах: «Известия Томского > политехнического университета», «Инженерно-физический журнал», «Измерительная техника», «Инженерная физика», депонированы в ВИНИТИ. Получено' одно авторское свидетельство на изобретение, свидетельство о регистрации программы для ЭВМ. Всего по материалам диссертации опубликовано 15'работ. Из них в соавторстве 13 работ с к.т.н., доцентом Загромовым Ю.А., доктором технических наук, профессором Караушем С.А., доктором физико-математических наук, профессором Кузнецовым Г.В. Содержание работы:
В первой главе рассмотрено современное состояние проблемы определения ТФХ материалов методом импульсного лазерного нагрева. Описаны-теоретические и экспериментальные способы определения ТФХ материалов. Показано, что, как правило, в качестве математических моделей теплопереноса используются одномерные упрощенные модели, либо модели, учитывающие незначительную часть факторов, определяющих протекание процессов теплопереноса при лазерном импульсном нагреве материалов. Отмечено, что практически не встречается работ с анализом погрешности методов, которые учитывали бы процессы распространения тепла в образце материала в двух- и трехмерной постановках. Установлена необходимость анализа погрешностей метода с использованием математических моделей, описывающих реальный перенос тепла при лазерном импульсном нагреве, с целью обоснованного выбора параметров лазерного нагрева, размеров исследуемого образца и условий теплообмена при проведении эксперимента по определению ТФХ.
Во второй- главе представлена постановка задач анализа погрешностей импульсного метода определения ТФХ и условий его применения. Описаны методы решения задач. Приведены результаты тестирования метода и алгоритма численного решения.
В третьей главе представлены результаты теоретического анализа основных погрешностей, определения ТФХ материалов методом импульсного лазерного нагрева, который проведен путем численных исследований с учетом изменения основных значимых параметров лазерного нагрева и условий численного эксперимента в диапазонах, соответствующих условиям реального эксперимента. Полученная информация является теоретической основой для определения границ использования метода импульсного лазерного нагрева и выделения реального допустимого диапазона изменения-параметров эксперимента (интенсивности и продолжительности лазерного воздействия, размеров зоны нагрева, геометрических размеров образца, условий теплообмена с окружающей средой). В заключении подведены итоги проведенных исследований. 1 t V
ЗАКЛЮЧЕНИЕ
Результаты выполненной работы и основные выводы заключаются в следующем.
1. Впервые сформулирован подход к анализу погрешностей и условий применения импульсных методов определения ТФХ материалов, в основе которого лежит численное моделирование реальных температурных полей, формирующихся в образцах в процессе экспериментов, из решения нестационарной задачи теплопроводности с нелинейными граничными условиями при высоких градиентах температур в. приповерхностном слое образца, высоких темпах нагрева, локальности теплового воздействия и сравнение вычисленных по результатам * численных экспериментов значений теплофизических характеристик материалов с истинными.
2. Проведен анализ влияния группы факторов на погрешности определения ТФХ материалов на основе одномерной, двухмерной и трехмерной моделей теплопроводности для области решения с неоднородными и нестационарными граничными условиями при высоких градиентах температур в приповерхностном слое образца и темпах нагрева, а также локальности теплового воздействия.
3. Результаты, полученные при решении одномерной задачи теплопроводности в условиях импульсного нагрева, позволяют сделать вывод об увеличении погрешности определения удельной теплоемкости, коэффициентов температуропроводности и теплопроводности с учетом теплообмена на всех открытых границах образца. Для материалов с низкой теплопроводностью при решении одномерной задачи установлены границы диапазонов возможного изменения интенсивности импульсного теплового потока, продолжительности импульса и толщины образца (в зависимости от типа материала), обеспечивающих достоверность определения ТФХ.
4. При решении двухмерной задачи теплопроводности в условиях импульсного нагрева показано влияние на погрешность определения ТФХ локальности воздействия лазерного импульса и неоднородности распределения энергии по поперечной координате. Установлено, что неучет этих факторов приводит к погрешностям определения удельной теплоемкости стали до 25 %, коэффициента температуропроводности до 13%, коэффициента теплопроводности до 38 %. Показано, что неоднородность пространственной энергетической плотности лазерного потока на поверхности пластины приводит к увеличению погрешностей. Так погрешность определения теплоемкости при значении коэффициента сосредоточенности теплового потока к = 0,06 составляет до 30 %, коэффициента температуропроводности до 14 %, коэффициента теплопроводности до 44 %,
5. Результаты, полученные при решении трехмерной задачи теплопроводности в условиях импульсного нагрева, позволяют сделать вывод о влиянии на погрешность определения ТФХ неодномерного I теплопереноса в образце при локальности воздействия лазерного луча: Установлено, что на методику определения ТФХ импульсным методом накладываются жесткие ограничения по размерам образца. Площадь нагреваемой поверхности образца должна соответствовать площади лазерного нагрева.
6. Результаты проведенных исследований показали, что теоретическая оценка погрешностей метода лазерного импульса создает условия для достоверного определения ТФХ любых материалов с установлением методических ошибок, которые можно называть систематическими. Кроме того, знание масштабов этих ошибок создает предпосылки для постановки новых задач по снижению погрешностей определения ТФХ этим методом.
1. Батаев A.A. Композиционные материалы: строение, получение, применение. М.: Логос, 2006. - 398с.
2. Болтон У. Конструкционные материалы: металлы, сплавы, полимеры, керамика, композиты: пер. с англ. 2-е изд. М.: Додека-ХХ1, 2007. - 320 с.
3. Басиев Т. Т., Осико В. В. Новые материалы для ВКР-лазеров // Успехи химии. — 2006. — Т. 75. —№ 10. — С. 939 955.
4. Валиев Р.З. Наноструктурные материалы, полученные интенсивной пластической деформацией. М.: Логос, 2000. - 272 с.
5. Кручинина, И. Ю. Новые материалы для совершенствования характеристик современных электрических машин // Пробл. создания и эксплуатации новых типов электроэнергетического оборудования. 2004 -№6.-С. 110-120.
6. Лякишев, Н. П. Промышленное применение наноматериалов // Наука Москвы и регионов . 2005. - № 1. - С. 44-45.
7. Neue Werkstoffe für zukünftige Gasturbinen // Konstr. + Giessen. 2005. — B. 30,-№ l.-L. 35.
8. Additive in Beschichtungssystemen kleine Dosierung, grosse Wirkung // Zement-Kalk-Gips Int. - 2007. -B. 60. - № 2. - L. 20 - 22.
9. Лыков A.B. Теория теплопроводности. M.: Высшая школа, 1967. - 599с.
10. Филиппов Л.П. Измерение тепловых свойств твердых и жидких металлов при высоких температурах. М.: Изд. Московского университета, 1967. -329с.
11. Шашков А.Г., Волохов Г.М., Абраменко Т.Н., Козлов В.П. Методы определения теплопроводности и температуропроводности. . М.: Энергия, 1973. - 336с.
12. Годовский Ю.К. Тепло физические методы исследования полимеров. — М.: Химия, 1976.-216 с.
13. Чудновский А.Ф. Теплофизика почв.- М.: Наука, 1976. -352с.
14. Платунов Е.С. Тепло физические измерения и приборы. Л.: Машиностроение, 1986. -256 с.
15. Чернышева Т.И., Чернышев В.Н. Методы и средства неразрушающегоVконтроля теплофизических свойств материалов. — М.: Машиностроение, 2001.-240 с.
16. Декуша Л.В. Грищенко Т.Г., Менделеева Т.В., Воробьев Л.И. , Декуша О.Л. Теоретическое обоснование прибора для экспресс-определения коэффициента теплопроводности твердых материалов // Промышленная теплотехника. 2004. - т. 26. - №1. - С. 76-82.
17. Фокин В.М. Неразрушающий контроль теплофизических характеристик строительных материалов. М.: Машиностроение, 2004. - 212с.
18. Cernuschi F., Bison P. G., Figari A., Marinetti S., Grinzato E. Thermal diffusivity measurements by photothermal < and thermographic techniques // Int. J. Thermophys. 2004. - Vol. 25. -No 2. - Pp. 439-457.
19. Пономарев С.В., Мищенко С.В., Дивин А.Г. Теоретические и практические аспекты теплофизических измерений: Монография. В 2 кн. Тамбов: Изд-во Тамб. гос. техн. ун-та, 2006. - Кн.1. - 206с.; Кн.2- 236 с.
20. Сысоев Э.В., Чернышов А.В. Методы и системы* неразрушающе го теплового контроля теплозащитных свойств многослойных строительных конструкций. -М.: Машиностроение 1, 2006. 104с.
21. Ищук И.Н. Идентификация теплофизических свойств материалов на основе использования вариационных методов исчисления // Измерительная техника. 2007 . - №7. - С .34-36.
22. Ищук И.Н. Принципы обнаружения дефектов при помощи термографии // Измерительная техника. 2008. - №4. - С. 49-53.
23. Parker W.J., Jenkins R.J., Butler С.Р. et al. Flash method of determining thermal diffusivity, heat capacity and thermal conductivity // J. of Appl. Physics. 1961. - Vol. 32. - No 9. - Pp. 1675-1684.
24. Regini F., Cezairliyan A. Pulse method of Thermal , diffusivity measurements // Hing. Temperatures. Hing Press. - 1973. - Vol: - No 5. - Pp. 481-501*.
25. Taylor R.E. Critical evaluation of flash method for measuring thermal diffusivity // Rev. int. Htes Temp, et refract. 1975. - V.12.-No 2. - Pp.141 -145'.
26. Варламов Г.Б., Дешко В.И., Карнацкий А .Я. Модифицированный метод Мгновенного источника для определения коэффициента температуропроводности // Промышленная теплотехника. 1987. - Т9. -№3. - С. 80-83.
27. Алабовский А.Р., Варламов Г.Б., Дешко В.И. Влияние радиационно-кондуктивного теплообмена при определении температуропроводности фтористого кальция методом мгновенного источника // Промышленная теплотехника.- 1988.-Т 10.- №4.- С. 83-88.
28. Тарзиманов А.А., Габитов Ф.Р, Поникарова И.Н., Юзмухамедов Ф.Л. Применение метода импульсного нагрева для различных теплофизических исследований // ИФЖ, 1993. -Т. 63. №4. - С. 436-442.
29. Пат. 2184952 РФ. МПК 8 G01N25/18. Способ неразрушающего контроля теплофизических характеристик материалов/ И.Н. Ищук, Т.А. Фесенко, В.В. Обухов. Заявлено 17.07.2000; Опубл.10.072002, Бюлл. №5.- Зс.
30. Фесенко А.И., Ищук И.Н. Метод определения теплофизических характеристик материалов при действии импульсного источника тепла // ИФЖ. 2000. - Т. 73. - №2. - С. 458.
31. Baba.T., Ono A. Improment of the laser flash method.to reduce uncertainty in thermal diffusivity measurements // Meas. Shi. Technol. 2001. - No 12.-Pp. 2046-2057.
32. Taylor R. Construction, of apparatus for heat pulse thermal diffusivity measurements from 300-3000K // J. Phys. E: Sci. Instrum. 1980. - V.13. -P.p 1193-1199.
33. Wisniewski S, Nagorski Z. Badania eksperymentalne wspolczynnika wyrownywania temperatury materialow konstrucyjnych mttoda impulsowq.— Pomiary Automatika kontrola. 1982. - Vol. 26. - №7. - Pp. 258-260.
34. Dai,J. М. . Development of a millisecond pulse-heating apparatus // Int. J. Thermophys. 2002. - Vol. 23. - No 5. - Pp. 1401-1405.
35. Wei Gaosheng, Zhang Xinxin, Yu Fan, Chen Kui. Thermal diffusivity measurements on insulation materials with the laser flash method // Int. J. Thermophys. 2006 - Vol. 27. - No 1. -Pp.235-243.
36. Arkin H. Thermal pulse decay method for simultaneous thermal conductivity and blood perfusion: a theoretical analysis // Trans. ASME. J. Biomech. Eng. -2001.-Vol. 108.-Pp. 208-214.
37. Zhang Haifeng, He Liquin, Cheng Shuxia, Zhai Zaiteng, Gao Dayong. A dualthermistor probe for absolute measurement of thermal diffusivity and thermal conductivity by the heat pulse method // Meas. Sci. and Technol. 2003- Vol. 14 .- No 8.- Pp. 1396-1401.
38. Schriempf J.T. A laser flash technique for determining thermal diffusivity of liquid metals elevated temperatures // Rev. Sci. instrum. 1972. - Vol. 43. -Pp. 781-786.
39. Ohta H. Thermal diffusivity measurements of molten salts using a three-layered cell by the laser flash method // Rev/, sci. Instrum. 1990. - Vol. 61. - No 10. -Pp. 2645-2649.
40. Tada Y., Harada M., Tanigaki M., Equchi W. Laser flash method for measuring thermal conductivity of liquids- application to low Thermal conductivity liquids //Rev. Sci. Instrum. 1978. - Vol. 49. - No 9. - Pp. 1305-1405.
41. Чистяков В.И. Импульсный метод определения коэффициента теплопроводности покрытий // Теплофизика высоких температур. 1973. — Т. 11. -№4. - С. 832-837.
42. Загромов Ю.А., Кац М.Д. Расчет теплофизических характеристик покрытия, нанесенного на металлическое основание, с использованием импульсного метода. ТПУ. Томск, 1987. - 17с. Деп. в ВИНИТИ 09.03.87, № 1686-В87.
43. Кац М.Д., Карауш С.А., Бугаев И.В. Влияние теплопотерь с поверхности двухслойного образца на измерение теплофизических характеристик импульсным методом // ИФЖ. 1991. - Т. 60. - №1. - С. 127-130.
44. Просунцов П.В., Резник C.B. Определение теплофизических свойств полупрозрачных материалов // ИФЖ. 1985. - Т. 49. - № 6. - С.971-976.
45. Kabayabaya Thomas, Yu Fan, Zhang Xinxin. Thermal diffusivity measurement of glass at high temperature by using flash method // J. Therm. Sei. 2004-Vol. 13.- No 1.- Pp. 91-96.
46. Baba T., Cezairliyan A. Thermal diffusivity of POCO AXM-5Q1 graphite in the range 1500 to 2500 К // Int. J. Thermophys. 1994. - Vol. 15. - Pp. 343364.
47. Baba T., Ono A. Development of laser flash apparaturs for thermal diffusivity standard . Proc 4th Asian thermophysical conf. - 1995. - Vol.4. - Pp. 581584.
48. Akoshima, M., Baba Т. Study on a thermal-diffusivity standard for laser flash method measurements // Int. J. Thermophys. 2006 . - Vol. 27. — No 4. -Pp. 1189-1203.
49. Неразрушающий контроль: справочник: в 8 т./ под ред. В.В. Клюева: Т.5. Кн.1: Вавилов В.П. Тепловой контроль. М.: Машиностроение, 2006: - 688 с.
50. Clark III L.M. , Taylor R.E. Radiation loss in the flash method of thermal diffusivity // J. Appl. Phys. 1975. - Vol. 46. - Pp. 714-719.
51. Watt D.A. Theory of thermal diffusivity by pulse technique // Br. J. Appl. Phys. -1966.-Vol. 17.-Pp. 231-240 .
52. Mendelsohn A.R. The effect of heat loss on the flash method* of determining thermal diffusivity // Appl. Phys. Letters. 1963. - Vol. 2. - No 1. - Pp. 19 -21.
53. Josell D., Warren J., Cezairliyan A. Comment on «Analysis for determining thermal diffusivity from thermal pulse experiments» // J. Appl. Phys. 1995. -Vol. 78. - No. 11 . - Pp. 6867-6869.
54. Taylor R.E., Cape A. Finite Pulse-Time effects in the flash' diffusivity // Technique Phys. Letters. 1964. -V. 5. -№ 10. - Pp. 212-213.
55. Larson K.B., Koyama K. Correction-for finite-pulse-time effects in ferry thin samples using the flash method of measuring thermal diffusivity //J. Appl. Phys.- 1967.-Vol. 38. NO 2 . -Pp.465-474.
56. Xye J., Liu X., Taylor R. The effects of a finite pulse time in the flash, thermal diffusivity // Int. J. Thermophys. 1993. - Vol. 14. - Pp. 123-133.
57. Taylor R.E., Clark L.M. III. Finite pulse-time effects in flash diffusivity method // Hing. Temp.-Hing Press. 1974. - Vol: 6. - Pp. 65-72.
58. Heckman R.C. Finite pulse-time and heat-loss effects in thermal diffusivity measurements // J. Appl. Phys. 1973. - Vol. 44. - No 4. - Pp. 1455-1460.
59. Андык B.C., Загромов Ю.А., Кац М.Д. Использование импульсного метода для измерения теплоемкости материалов. ТПУ. Томск, 1987.- 8 с. Деп. в ВИНИТИ 12.11.87, № 17999-В87.
60. Azumi T., Takahaschi Y. Novel finite pulse-width correction in flash method thermal diffusivity measurement // Rev. Sci. Instrum. 1981. - Vol. 52. - No 9.- Pp.1411-1413.
61. Перевозчиков C.M., Загребин Л.Д. Анализ погрешностей измерения теплофизических характеристик неограниченной пластины импульсным методом // Измерительная техника . — 2001. — №12. С. 39-43.
62. Beedham R., Dalmple I.P. Nit measurement of thermal diffusivity bu flash method. An investigation into errors arising from boundary conditions //Revue Int. Hautes Temp. Refract. 1970. - Vol. 7. - Pp. 278-283.
63. Mckay J., Schriempf J.T. Corrections for nonuniform surface-heating errors in flash-method thermal-diffusivity measurements //J. Appl. Phys. 1976 .-V47. — No 4.-Pp. 1668-1671.
64. Клименко М.М., Кржижановский Р.Е., Шерман В.Е. Анализ методических погрешностей измерения температуропроводности импульсным методом с применением лазера // Измерительная техника. — 1980. №6. - С.40-42.
65. Gembarovic J., Gembarovic Jr. Nonlinear Effects in Laser Flash Thermal Diffusivity measurements// Int. J. of Thermophysics. 2004. - Vol. 25. - No 4. -Pp. 1253-1260.
66. Исаев К.Б. Автоматизированный комплекс для исследования теплопроводности материалов // Промышленная теплотехника. 2001. -Т. 23. — №2. - С. 46-55.
67. Cezairliyan A.,Baba Т.А., Taylor R. A Hing-temperature laser-pulse thermal diffusivity apparaturs // Int. J. Thermophys. 1994. - Vol.15 . - Pp. 317-341.
68. Бояринов A.E., Способы и программно-управляемые средства для теплофизических измерений // Материалы 2-й научной конференции Тамбовск. Гос. Тех. Ун-та. Тамбов, 1995. - С. 92-93.
69. Перевозчиков С.М., Загребин Л.Д. Автоматизированная система измерения теплофизических параметров металлов и сплавов // Измерительная техника .- 1998. -№3. С.155-158.
70. А.с. 1465751СССР: МКИ3 G 01 N 25/18. Устройство для определения коэффициента температуропроводности твердых материалов: / В.П. Вавилов, А.И. Иванов, В.В. Ширяев, А.И. Суппес, М.Д. Кац. Заяв. 1.06.87; опубл. 15.03.89, Бюл.№10.- Зс.
71. Cowan R. Proposed method of thermal diffusivity at Hing Temperatures // J. of Appl. Phisics. 1961 .-V. 32. - No 7. - Pp. 189-196.
72. Пелецкий В.Э. Исследование теплофизических свойств веществ в условиях электронного нагрева. — М.: Наука, 1983.- 92с.
73. Григорьянц А.Г. Основы лазерной обработки материалов. М.: Машиностроение, 1989.-304с.
74. Материаловедение и технология металлов / М.Г. Карпман, В.М.Матюнин, B.C. Гаврилюк и др. М.: Высшая школа, 2002. 638с.
75. Загромов Ю. А., Кац М. Д. Усилитель постоянного напряжения для регистрации малых изменений температуры // Материалы межвузовской конференции «Теплофизика и гидродинамика технологических процессов. -Томск, 1989. -С.128-131.
76. Kobayashi М., Baba Т., Ono A. Infrared radiation thermometer for laser flash thermal diffusivity // Thermophys. Prop. 1994. - Vol. 8. - Pp. 143-148.
77. Baba Т., Kobayashi M., Ono A. Nonlinearity effect of radiation thermometry in laser flash thermal diffusivity // Japan. J. Thermophys. 1994. - Prop. 7. -Pp.156-158.
78. Каспаров K.H. Исследование динамики температуры при импульсном нагреве методом фотоэмиссионной пирометрии // Измерительная техника. -2006,-№9.- С. 46-48.
79. Резников Б.И., Менде Н.П„ Попов П.А., Сахаров В.А., Штейнберг A.C. Измерение теплового потока по измерения температуры поверхности в импульсных газодинамических процессах // Письма в ЖТФ, 2008. Т. 34. - вып. 5. - С. 49-54.
80. Карслоу Г., Егер Д. Теплопроводность твердых тел. М.: Наука, 1964. -487с.
81. Троицкий О.Ю. Импульсный неразрушающий контроль теплофизических характеристик материалов // Механика композитных материалов. 1992. -№ 4. -С.534- 538.
82. Рыкалин H.H., Углов A.A., Кокора А.Н. Лазерная обработка материалов. -М.: Машиностроение, 1975. 296 с.
83. Ищук И.Н. Принципы обработки информации на основе численного решения нелинейной задачи теплопроводности. Ч 2. Коэффициентная задача функциональной оптимизации в вариационной постановке // Измерительная Техника, 2008. №2 - С. 48-50.
84. Самарский A.A. Теория разностных схем. М.: Наука, 1983 - 616с.
85. Пасконов В.М., Полежаев В.И., Чудов JI.A. Численное моделирование процессов теплообмена. М.: Наука, 1984. -288с.
86. Вержбицкий В.М. Основы численных методов. — М.: Высшая школа, 2002. 840 с.
87. Кузнецов Г.В., Шеремет М.А. Разностные методы решения задач теплопроводности. — Томск: Изд-во ТПУ, 2007. — 173с.
88. Коздоба JI.A. Методы решения нелинейных задач теплопроводности. М.: Наука, 1975.-227с.
89. Марчук Г.И. Методы расщепления . М.: Наука, 1988. - 264с.
90. Чиркин B.C. Теплофизические свойства материалов ядерной техники. М.: Атомиздат, 1968. - 484 с.
91. Технологические лазеры. Справочник в 2-х т. / Под ред. Г.А. Абильсиинова. М.: Машиностроение, 1991. - Т. 1.-431 с.
92. Кузнецов Г.В., Кац М.Д. Анализ погрешностей определения импульсными методами теплофизических характеристик конструкционных материалов // Известия Томского политехнического университета. — 2008. — Т.313. №4. -С. 5-8.
93. Кириллин В.А., Сычев В.В., Шейдлин А.Е. Техническая термодинамика. -М.: Энергоатомиздат, 1983.-416 с.
94. Кржижановский Р.Е, Штерн З.Ю. Теплофизические свойства неметаллических материалов (карбиды). Справочная книга. JL: Энергия, 1976. - 120с.
95. Шленский О.Ф. Терморазрушение материалов. Полимеры и композиты при интенсивном нагреве. М.: Энергоатомиздат, 1996. - 288с.
96. Сабсай О.Ю., Чалая Н.М. Технологические свойства термопластов (обзор) // Пластические массы. 1992. - № 1. - С. 5-13.
97. Кузнецов Г.В., Кац М.Д. Об условиях применения импульсных методов определения теплофизических характеристик конструкционных материалов // Известия Томского политехнического университета. 2008. -Т. 312.-№4.-С. 10-13.
98. Кузнецов Г.В., Кац М.Д. Анализ погрешности определения теплофизических характеристик конструкционных материалов импульсным методом при конечной длительности импульса лазерного луча // Инженерная физика. 2009. - №8. - С. 45-^18.
99. Кац М.Д. О погрешностях определения теплофизических характеристик конструкционных материалов Импульсным методом // Материалы 14-й Всероссийской научно-технической конференции «Энергетика: Экология, Надежность, безопасность». Томск, 2008. — С. 79-81.
100. Кузнецов Г.В., Кац М.Д. Оценка погрешности определения теплофизических характеристик карбида циркония импульсным методом при пространственной неоднородности лазерного луча. ТПУ. — Томск, 2009.-Юс. Деп. ВИНИТИ 01.07.09, № 419 В2009.
101. Зайдель А.Н. Погрешности измерения физических величин. JL: Наука, 1985. -112 с.
102. Кузнецов Г.В., Кац М.Д. Теоретический анализ методических погрешностей определения теплофизических характеристик конструкционных материалов импульсным методом в образце конечных размеров // Измерительная техника. — 2009. — №4. — С. 34—36.
103. Кузнецов Г.В., Кац М.Д. Погрешности определения теплофизических характеристик карбида циркония импульсным методом, обусловленные локальным воздействием лазерного луча. ТПУ, 2009. Томск. - 12с. Деп. ВИНИТИ 18.06.09, № 382 - В2009.