Оценка термических характеристик геоматериалов по данным компьютерного моделирования

В статті представлені результати досліджень теплових властовостей геоматеріалів в широкому діапазоні температур. За даними термодинамічних розрахунків досліджено вплив факторів на ефективну та істинну теплоємність. За допомогою метода Сімпсона обчислено теплові ефекти піролізу....

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2011
Автор: Клюев, Э.С.
Формат: Стаття
Мова:Russian
Опубліковано: Інститут геотехнічної механіки імені М.С. Полякова НАН України 2011
Назва видання:Геотехническая механика
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/33585
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Оценка термических характеристик геоматериалов по данным компьютерного моделирования / Э.С. Клюев // Геотехническая механика: Межвед. сб. науч. тр. — Днепропетровск: ИГТМ НАНУ, 2011. — Вип. 95. — Бібліогр.: 15 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-33585
record_format dspace
spelling irk-123456789-335852012-05-29T12:55:51Z Оценка термических характеристик геоматериалов по данным компьютерного моделирования Клюев, Э.С. В статті представлені результати досліджень теплових властовостей геоматеріалів в широкому діапазоні температур. За даними термодинамічних розрахунків досліджено вплив факторів на ефективну та істинну теплоємність. За допомогою метода Сімпсона обчислено теплові ефекти піролізу. The results of investigations of geomaterials’ thermal properties in wide temperature range were presented in the article. Due to data of thermodynamic calculations the influence of factors on effective and real heat capacities was investigated. By means of Simpson method the thermal effects of pyrolysis were determined. 2011 Article Оценка термических характеристик геоматериалов по данным компьютерного моделирования / Э.С. Клюев // Геотехническая механика: Межвед. сб. науч. тр. — Днепропетровск: ИГТМ НАНУ, 2011. — Вип. 95. — Бібліогр.: 15 назв. — рос. http://dspace.nbuv.gov.ua/handle/123456789/33585 662.62:53 ru Геотехническая механика Інститут геотехнічної механіки імені М.С. Полякова НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
description В статті представлені результати досліджень теплових властовостей геоматеріалів в широкому діапазоні температур. За даними термодинамічних розрахунків досліджено вплив факторів на ефективну та істинну теплоємність. За допомогою метода Сімпсона обчислено теплові ефекти піролізу.
format Article
author Клюев, Э.С.
spellingShingle Клюев, Э.С.
Оценка термических характеристик геоматериалов по данным компьютерного моделирования
Геотехническая механика
author_facet Клюев, Э.С.
author_sort Клюев, Э.С.
title Оценка термических характеристик геоматериалов по данным компьютерного моделирования
title_short Оценка термических характеристик геоматериалов по данным компьютерного моделирования
title_full Оценка термических характеристик геоматериалов по данным компьютерного моделирования
title_fullStr Оценка термических характеристик геоматериалов по данным компьютерного моделирования
title_full_unstemmed Оценка термических характеристик геоматериалов по данным компьютерного моделирования
title_sort оценка термических характеристик геоматериалов по данным компьютерного моделирования
publisher Інститут геотехнічної механіки імені М.С. Полякова НАН України
publishDate 2011
url http://dspace.nbuv.gov.ua/handle/123456789/33585
citation_txt Оценка термических характеристик геоматериалов по данным компьютерного моделирования / Э.С. Клюев // Геотехническая механика: Межвед. сб. науч. тр. — Днепропетровск: ИГТМ НАНУ, 2011. — Вип. 95. — Бібліогр.: 15 назв. — рос.
series Геотехническая механика
work_keys_str_mv AT klûevés ocenkatermičeskihharakteristikgeomaterialovpodannymkompʹûternogomodelirovaniâ
first_indexed 2025-07-03T14:16:20Z
last_indexed 2025-07-03T14:16:20Z
_version_ 1836635579229405184
fulltext УДК 662.62:53 Э.С. Клюев, асп. (ИГТМ) ОЦЕНКА ТЕРМИЧЕСКИХ ХАРАКТЕРИСТИК ГЕОМАТЕРИАЛОВ ПО ДАННЫМ КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ В статті представлені результати досліджень теплових властовостей геоматеріалів в широкому діапазоні температур. За даними термодинамічних розрахунків досліджено вплив факторів на ефективну та істинну теплоємність. За допомогою метода Сімпсона об- числено теплові ефекти піролізу. THE ESTIMATION OF THERMAL CHARACTERISTICS OF GEOMATERIALS USING DATA OF COMPUTER SIMULATION The results of investigations of geomaterials’ thermal properties in wide temperature range were presented in the article. Due to data of thermodynamic calculations the influence of factors on effective and real heat capacities was investigated. By means of Simpson method the thermal effects of pyrolysis were determined. Актуальность работы. Одним из интенсивно развиваемых направлений использования угля является термохимическая деструкция, характер которой зависит, в основном, от свойств перерабатываемого материала и конструкции тепловых установок [1-3]. Основная проблема, которая должна при этом ре- шаться, заключается в поиске рациональных и эффективных путей такой пе- реработки, что требует выявления оптимальных условий процесса и разра- ботки экспресс-методов оценки термических превращений. С точки зрения технологии большой интерес представляют данные о теплофизических свой- ствах сырья и их изменении в ходе термической переработки, так как от этих свойств в значительной мере зависит продолжительность и энергетическую эффективность процесса. Поэтому их оценка имеет большое научно- практическое значение, так как расход тепла на высокотемпературную пере- работку определяется не только количеством перерабатываемого материала, но и его теплоемкостью. Степень исследования проблемы. Разработкой методов теоретического и экспериментального определения теплоемкости твердого топлива, а также построением обобщенных зависимостей, учитывающих влияние ряда факто- ров, занимались в свое время как отечественные, так и зарубежные ученые. Первые работы в этом направлении были начаты еще в середине прошлого века. За это время накоплен огромный экспериментальный материал, иссле- довано различные типы топлив: торф, бурый уголь, горючие сланцы, камен- ный уголь разной степени метаморфизма и др. [4,5] Первые шаги к выявлению взаимосвязи между строением углей и их теп- лоемкостью были сделаны В. Фрицем и Г. Мозером [6]. Опытным путем в ка- лориметрической установке с помощью метода смешения было установлено среднюю удельную теплоемкость 40 проб каменных углей и сделан вывод о закономерном уменьшении этого параметра с увеличением степени метамор- физма. На основе полученных данных предложено эмпирическое уравнение  dVс 008,011012  , где с – удельная теплоемкость сухой массы угля, Дж/(кгК); V d – выход лету- чих веществ на сухую массу топлива, %. Идеи [6] были развиты в исследованиях Д.В. ван Кревелена [7]. Его работа заключалась в применении к органической массе углей эмпирического пра- вила Неймана-Коппа, согласно которому молярная теплоемкость твердых со- единений равна сумме атомных теплоемкостей входящих элементов. В последствии в трудах А.А. Агроскина [8,9] было учтено влияние внеш- ней влаги и минеральных примесей на среднюю удельную теплоемкость и предложено формулу для приблизительного ее расчета при обычных темпе- ратурах   rrrd WAВVс 4187794008,0112,10  , где В r – содержание органической массы в рабочем топливе, %; A r – зольность рабочего топлива, %; W r – содержание влаги в рабочем топливе, %. В работе [10] установлено, что теплоемкость угольных смесей является аддитивной величиной. Расчетным путем доказано, что с повышением степе- ни метаморфизма угля максимум теплоемкости сдвигается в сторону более высоких температур. Так, если у бурого угля максимум значения теплоемко- сти получается при 350 о С, то у углей марок Г и Ж – при 400 о С, а у К и ОС – при 450-470 о С. Как известно [11], в термохимии и термодинамике твердого топлива раз- личают эффективную и истинную теплоемкости, а также тепловой эффект ре- акции. Немаловажное значение имеет также разработка методики определе- ния этих характеристик в широком диапазоне температур. Этими проблемами занимались А.А. Агроскин, В.Б. Глейбман, Е.И. Гончаров, Н.С. Грязнов и др. Базируясь на методе диатермической оболочки, в работе [12] определены эф- фективная и истинная теплоемкость, а также теплота реакций пиролиза углей Донбасса различных марок. Исследования проводились с сухими пробами, измельченными до 0,25 мм, при скорости нагрева 10 град/мин в непрерывном токе сухого азота, очищенного от примеси кислорода. Установлено, что эф- фективная теплоемкость исследуемых проб линейно возрастает до темпера- туры начала разложения. При этом истинная теплоемкость углей с повыше- нием температуры равномерно растет, а абсолютные ее значения уменьшают- ся с ростом степени метаморфизма. Комплексный подход к определению удельной теплоемкости каменных углей был предпринят Л.И. Гладковым и А.Н. Лебедевым [13], в результате которого выведены эмпирические зависимости, учитывающие влияние золь- ности, влажности, выхода летучих веществ и температуры процесса. Для рас- чета теплоемкости минеральных включений углей См, кДж/(кгК) предложено линейное уравнение вида ТСм 5105071,0  , (1) где Т – температура процесса, К. А.А. Агроскин [14] для расчета теплоемкости влажного угля при условии, что влага находится в несвязанном виде, предложил использовать выражение   r c r вл WСWС  1 . (2) Расчет органической массы углей следует проводить по формуле  TVС daf орг   303)13(10368,0 6 . (3) Зависимости (1) – (3) являются, по-видимому, одними из наиболее надеж- ных из ныне предложенных, что дает возможность рассчитывать теплоем- кость реальных углей с учетом их зольности и влажности. Существенным недостатком выполненных работ является то, что приве- денные обобщенные закономерности изменения термических характеристик твердых топлив получены лишь в узком диапазоне температур. Так, выраже- ние (3) применимо лишь до температуры 500 о С. Кроме того, в литературе от- сутствуют какая-либо информация об измерении теплоемкости при высоких температурах (1500 о С и выше) ввиду сложности комплектации эксперимен- тального оборудования и трудностей, связанных со специфической природой углей. Поэтому данные, полученные с помощью компьютерного моделирования термохимических высокотемпературных процессов, помогут провести более точный анализ теплофизических свойств исследуемых материалов, в т.ч. теп- лопроводности. В связи с этим целью работы является изучение влияния температуры до 2800 К на изменение теплоемкости геоматериалов для их дальнейшего энер- гетического использования с минимальными энергозатратами. При химических превращениях любых веществ наблюдается выделение или поглощение тепла, что существенно изменяет температурную зависи- мость теплоемкости: процессы, сопровождающиеся поглощением тепла (эн- дотермическим эффектом) повышают эффективную теплоемкость термоди- намической системы, а экзотермические процессы, протекающие с выделени- ем тепла, снижают ее. В этой связи имеет смысл говорить об эндотермиче- ском максимуме и экзотермическом минимуме процесса. Если привести сис- тему к некоторой температуре Т * и подвергнуть ее достаточно продолжитель- ному выдерживанию при этой температуре, то по мере того как система будет переходить в новое равновесное состояние теплоемкость ее будет изменятся, стремясь к некоторому значению, которое можно было бы назвать равновес- ной (истинной) теплоемкостью термодинамической системы при температуре Т * . В целях определения теплоемкости геоматериалов методом компьютерно- го моделирования были проведены исследования углей Львовско-Волынского бассейна и отходов углеобогащения (шлам) ЦОФ «Червоноградская». Резуль- таты элементного и технического анализов исследованных проб приведены в таблице 1. Таблица 1 – Результаты элементного и технического анализа исследован- ных проб Исследуемый геоматериал С р ,% H р , % N р , % O р , % р общS ,% W р , % A d ,% V daf , % ММ * , % Сапропелевый уголь 34,8 2,5 0,7 3,2 0,3 1,8 56,7 23,4 62,8 Гумусовый уголь 78,4 5,9 1,7 9,9 0,5 2,2 1,4 33,1 1,8 Шлам, проба № 2 26,9 3,1 2,8 19,5 1,7 2,1 43,9 17,9 49,6 Шлам, проба № 3 26,2 3,1 2,8 19,8 1,7 1,6 44,8 20,2 50,6 Шлам, проба № 8 25,6 3,0 2,9 19,9 1,8 1,4 45,4 18,5 51,4 Шлам усредненный 26,8 1,8 2,9 20,4 2,8 1,2 44,1 21,8 50,7 * ММ – содержание минеральных примесей в пробе. На рисунках 1,2 представлены графики зависимости эффективной тепло- емкости для различных проб исследуемых материалов от температуры. 0 1 2 3 4 5 6 7 8 9 10 0 500 1000 1500 2000 2500 3000 Т, К С, кДж/(кг К) 1 2 1 – гумусовый уголь; 2 – сапропелевый уголь Рис. 1 – Температурная зависимость эффективной теплоемкости углей Львовско- Волынского бассейна 0 1 2 3 4 5 6 7 8 9 10 0 500 1000 1500 2000 2500 3000 Т, К С, кДж/(кг К) 0 1 2 3 4 5 6 7 8 0 500 1000 1500 2000 2500 3000 Т, К С, кДж/(кг К) а) б) 0 2 4 6 8 10 12 0 500 1000 1500 2000 2500 3000 Т, К С, кДж/(кг К) в) г) а) усредненная проба; б) проба № 2; в) проба № 3; г) проба № 8 Рис. 2 – Температурная зависимость эффективной теплоемкости отходов углеобогащения ЦОФ «Червоноградская» Анализ кривых показывает, что первый максимум эффективной теплоем- кости, обусловленный эндотермическими реакциями разложения глинистого вещества и пирита, имеет место при температуре около 1000 К. Для углей (рис. 1) этот максимум смещен на 100-200 град в сторону более низких тем- ператур по причине их невысокой степени метаморфизма. При этом даль- нейшее снижение эффективной теплоемкости в интервале 1000-1400 К свиде- тельствуют об экзотермических реакциях поликонденсации органической массы угля. Второй, значительно меньший по абсолютной величине, макси- мум вызван эндотермическими реакциями разложения карбонатов и процес- сом выделения двуокиси углерода. Дальнейший нагрев приводит к росту ве- личины эффективной теплоемкости до 9-9,1 кДж/(кгК), что непосредственно связано с увеличением энергозатрат на процесс термической переработки. При этом для гумусового угля (рис. 1, кривая 1) абсолютные значения тепло- емкости существенно выше, чем у сапропелевого угля (рис. 1, кривая 2), вследствие низкого содержания минеральных примесей (табл. 1), теплоем- 0 1 2 3 4 5 6 7 8 9 10 0 500 1000 1500 2000 2500 3000 Т, К С, кДж/(кг К) кость которых ниже теплоемкости органической массы угля. Кроме того, этот факт предопределяет наличие двух дополнительных эндоэффектов в пробе сапропелевого угля при температуре 1900 К и 2100 К. Что касается характера изменения кривых для отходов обогащения угля (рис. 2), то в данном случае наблюдается иная картина. Для всех проб после 1000 К отмечается от 4 до 6 эндотермических эффектов, свидетельствующих, согласно табл. 1, о наличии значительных количеств минеральных примесей. Относительно низкие значения эффективной теплоемкости в интервале 1500- 2000 К объясняются, как отмечалось, экзотермическими реакциями поликон- денсации присутствующего в шламах угля. Кроме того, при температуре 2000 К для всех проб наблюдается эффект, достигающий величины 3 кДж/(кгК). В ходе дальнейшего нагрева, способствующего протеканию эндо- и экзотерми- ческих реакций, на кривых заметны переходные тепловые процессы, в итоге которых при температуре 2700 К достигается очередной эндотермический максимум с последующими экзоэффектом. Следует подчеркнуть, что это яв- ление наблюдается только у исследуемых отходов, а в углях – отсутствует. Результаты определения истинной теплоемкости исследованных геомате- риалов приведены в табл. 2. Установлено, что истинная теплоемкость всех проб монотонно возрастает с повышением температуры. Более высокие значения истинной теплоемкости гумусового угля по сравнению с теплоемкостью других материалов объясня- ются его низкой зольностью. Компьютерная обработка этих эксперименталь- ных данных позволила установить корреляционную связь между истинной теплоемкостью и температурой нагрева отходов в режиме до 2800 К, которая описывается уравнением 91,0;1021060886,1 2274   RTTCp . (4) Воспользовавшись первым законом термодинамики, можно получить формулу Кирхгофа, выражающую зависимость теплового эффекта от темпе- ратуры химического превращения      dTССТEТE Т Т ххо о   , (5) где Е(Т)и Е(То)– тепловые эффекты при температуре Т и То соответственно, Дж/кг; Cх и Сх – теплоемкости системы в исходном и конечном состояниях, Дж/(кгК). Применительно к процессам термической переработки твердого топлива есть смысл принять То достаточно низкой (например, 20 о С), чтобы пренеб- речь химическими превращениями при этой температуре. Тогда формула (5) запишется в следующем виде     Т Т эф Т Т истТТ оо о dTС dTС E , (6) где  ТТо E  – суммарный тепловой эффект реакции пиролиза в интервале тем- ператур от То до Т, Дж/кг; Сист – истинная теплоемкость, Дж/(кгК); Сэф – эф- фективная теплоемкость, Дж/(кгК). Таблица 2 – Истинная теплоемкость геоматериалов Темпера- тура про- цесса, К Сапропе- левый уголь Гумусо- вый уголь Шлам усред- ненный Шлам, проба № 2 Шлам , про- ба № 3 Шлам, проба № 8 500 1,254 1,597 1,203 1,312 1,306 1,296 600 1,380 1,804 1,286 1,409 1,402 1,390 700 1,498 2,002 1,353 1,501 1,493 1,478 800 1,632 2,190 1,434 1,607 1,598 1,580 900 1,631 2,344 1,467 1,647 1,636 1,616 1000 1,681 2,447 1,489 1,679 1,668 1,647 1100 1,715 2,514 1,499 1,697 1,685 1,663 1200 1,742 2,564 1,514 1,716 1,703 1,682 1300 1,765 2,602 1,531 1,735 1,723 1,700 1400 1,785 2,634 1,547 1,752 1,739 1,717 1500 1,802 2,662 1,549 1,768 1,751 1,727 1600 1,818 2,688 1,559 1,778 1,77 1,741 1700 1,819 2,713 1,562 1,782 1,768 1,745 1800 1,752 2,737 1,525 1,752 1,741 1,721 1900 1,759 2,761 1,543 1,772 1,762 1,741 2000 1,765 2,785 1,548 1,780 1,770 1,749 2100 1,766 2,808 1,557 1,793 1,784 1,764 2200 1,688 2,831 1,546 1,782 1,766 1,751 2300 1,694 2,854 1,537 1,770 1,770 1,754 2400 1,698 2,876 1,545 1,776 1,780 1,768 2500 1,701 2,899 1,531 1,767 1,772 1,761 2600 1,701 2,921 1,528 1,765 1,770 1,759 2700 1,696 2,940 1,522 1,759 1,763 1,746 2800 1,253 2,955 1,508 1,745 1,743 1,723 Проанализировав полученные экспериментальные данные, можно сделать вывод о том, что в данном случае левая часть выражения (6) представлена не- явно, и произвести расчет теплового эффекта в аналитическом виде не пред- ставляется возможным. Поэтому следует воспользоваться численными мето- дами интегрирования. В данной работе применялся метод Симпсона, осно- ванный на приближении подынтегральной функции на отрезке интерполяци- онным полиномом второй степени [15]. Полученные результаты приведены в табл. 3, которые свидетельствуют о значительном поглощении подводимого тепла в процессе пиролиза до температуры 2800 К, поскольку все величины по знаку оказались положительны. Суммарный тепловой эффект реакций раз- ложения исследованных проб составляет от +2801,9 (гумусовый уголь) до +5955,2 кДж/кг (проба шлама № 3). Таблица 3 – Тепловой эффект разложения геоматериалов Исследуемый геома- териал Суммарный те- пловой эффект пиролиза, кДж/кг Сапропелевый уголь +4261,2 Гумусовый уголь +2801,9 Шлам усредненный +3673,2 Шлам, проба № 2 +5521,5 Шлам, проба № 3 +5955,2 Шлам, проба № 8 +4539,1 Характерной особенностью изучаемых закономерностей является то, что экзотермический эффект реакций пиролиза оказался равным нулю для всех проб геоматериалов, что свойственно лишь породе. Поэтому выделяемого в процессе деструкции тепла не достаточно для самостоятельного поддержания процесса термической деструкции, что свидетельствует о необходимости до- полнительно подводить тепло извне. Выводы. 1. Определены закономерности изменения эффективной и истинной теп- лоемкостей, а также значения суммарного теплового эффекта термических процессов при переработке сапропелевого и гумусового углей Львовско- Волынского бассейна, а также отходов углеобогащения Червоноградской ЦОФ в интервале температур до 2800 К. 2. Установлено, что с постепенным ростом температуры кривые эффек- тивной теплоемкости характеризуются наличием ряда экстремальных точек, обусловленных одновременным протеканием реакций деструкции органиче- ской и минеральной частей исследуемого материала. 3. Исследованиями показано, что истинная теплоемкость монотонно рас- тет с повышением температуры, а их абсолютные значения для разных мате- риалов незначительно отличаются друг от друга. Исключение составляет лишь гумусовый уголь, для которого эта величина в 1,3-1,5 раз больше, вследствие его низкой зольности. 4. Путем математической обработки данных компьютерного моделирова- ния установлена корреляционная связь между истинной теплоемкостью и температурой нагрева, что может быть использовано для расчета удельных энергетических затрат при термической переработки углесодержащих мате- риалов. 5. С применением численного метода Симпсона вычислены суммарные тепловые эффекты пиролиза исследованных геоматериалов, которые свиде- тельствуют о значительном поглощении подводимого тепла в процессе пиро- лиза, поскольку все величины по знаку оказались положительны. СПИСОК ЛИТЕРАТУРЫ 1. Mountouris A. Plasma gasification of sewage sludge: Process development and energy optimization / A. Mountouris, E. Voutsas, D. Tassios // Energy Conversion and Management. – 2008. – № 49. – рр. 2264 – 2271. 2. Пат. 56306 Україна, МПК 7 С 10 J 3/46, F 23 C 3/00. Спосіб термічної переробки шламів та неконди- ційного вугілля / Булат А.Ф., Возіянов В.С., Курносов С.А., Клюєв Е.С., Приходченко В.Л., Слащова О.А., Осінній В.Я., Давидов С.Л.; заявник та патентовласник Інститут геотехнічної механіки ім. М.С. Полякова НАН України. – u 2010 07414; заявл. 14.06.2010; опубл. 10.01.2011, Бюл. №1. 3. Осенний В.Я. Плазменные процессы при переработке отходов углеобогатительных фабрик / В.Я. Осенний, Э.С. Клюев, С.Л. Давыдов // Сучасні ресурсоенергозберігаючі технології гірничого виробниц- тва: Науково-виробничий збірник / Кременчуцький державний університет ім. М. Остроградського. – Кре- менчук, 2010. – Вип. 1/2010 (5). – С. 112 – 120. 4. Глущенко И.М. Термический анализ твердых топлив / И.М. Глущенко. – М.: Металлургия, 1968. – 192 с. 5. Влияние условий нагрева низкосортных углей и отходов углеобогащения на продукты термодеструк- ции / В.Л. Приходченко, Е.А. Слащева, В.Я. Осенний [и др.] // Геотехническая механика: Межвед. сб. науч. трудов / Институт геотехнической механики им. Н.С. Полякова НАН Украины. – Днепропетровск, 2010. – Вып. 89. – С. 65 – 72. 6. Fritz W. Spezifische Wärme, Wärmeleitfähigkeit und Temperaturleitfähigkeit von Steinkhole und Koks / W. Fritz, H. Moser // Feuerungstechnik. – 1940. – s. 97 – 107. 7. Van Krevelen D.W. Physikalische Eigenschaften und Chemische Struktur der Steinkohle-Brennst / D.W. Van Krevelen // Chemie. – 1953. – Bd. 34. – N 11/12. – s. 167 – 182. 8. Агроскин А.А. Теплоемкость углей / А.А. Агроскин, Е.И. Гончаров // Кокс и химия. – 1965. – № 7. – С. 8 – 13. 9. Агроскин А.А. Теплоемкость и особенности молекулярной структуры веществ углей / А.А. Агроскин, В.Б. Глейбман // Химия твердого топлива. – 1975. – № 3. – С. 61 – 67. 10. Николаев И.Н. Зависимость теплоемкости углей и угольных смесей от температуры их нагрева / И.Н. Николаев, Н.И. Козлова, В.А. Калинкина, А.А. Степанчиков // Кокс и химия. – 1977. – № 5. – С. 6 – 10. 11. Агроскин А.А. Теплофизика твердого тела / А.А. Агроскин, В.Б. Глейбман. - М.: Недра, 1980. – 256 с. 12. Агроскин А.А. Теплоемкость и теплота реакций пиролиза углей Донбасса / А.А. Агроскин, Е.И. Гон- чаров, Л.А. Макеев [и др.] // Кокс и химия. – 1970. – № 5. – С. 8 – 13. 13. Гладков Л.И. Теплоемкость твердого топлива и угольной пыли / Л.И. Гладков, А.П. Лебедев // Изв. ВТИ. – 1948. – № 8. – С. 18 – 20. 14. Агроскин А.А. Физика угля / А.А. Агроскин. – М.: Недра, 1965. – 351 с. 15. Мудров А.Е. Численные методы для ПЭВМ на языках Бейсик, Фортран и Паскаль / А.Е. Мудров. – Томск: МП «РАСКО», 1991. – 272 с.