Обоснование применения грунтовых теплообменников для перегрева паровой фазы сжиженного углеводородного газа Текст научной статьи по специальности «Энергетика и рациональное природопользование»

ТЕПЛОСНАБЖЕНИЕ, ВЕНТИЛЯЦИЯ, КОНДИЦИОНИРОВАНИЕ ВОЗДУХА, ГАЗОСНАБЖЕНИЕ И ОСВЕЩЕНИЕ

УДК 696.2:662.767.3

ОСИПОВА НАТАЛИЯ НИКОЛАЕВНА, докт. техн. наук, доцент, osnat75@mail. ги

БЫЧКОВА ИРИНА МИХАЙЛОВНА, аспирант, [email protected]

СВИРИДОВА КСЕНИЯ СЕРГЕЕВНА, магистрант, [email protected]

Саратовский государственный технический университет имени Ю.А. Гагарина,

410054, г. Саратов, ул. Политехническая, 77

ОБОСНОВАНИЕ ПРИМЕНЕНИЯ ГРУНТОВЫХ ТЕПЛООБМЕННИКОВ ДЛЯ ПЕРЕГРЕВА ПАРОВОЙ ФАЗЫ СЖИЖЕННОГО УГЛЕВОДОРОДНОГО ГАЗА

Предложена схема газоснабжения на базе сжиженного углеводородного газа с естественной регазификацией, оснащенной грунтовым теплообменником, позволяющая обеспечить безгидратное редуцирование паровой фазы перед подачей к газоиспользую-щим приборам. Исследован процесс теплообмена грунтового теплообменника с грунтовым массивом. Установлено, что наличие грунтового теплообменника позволит обеспечить перегрев паровой фазы сжиженного углеводородного газа и исключить образование гидратов в регуляторах давления при дросселировании паров.

Ключевые слова: сжиженный углеводородный газ; грунтовый теплообменник; естественная регазификация; паровая фаза; перегрев; гидратообразование.

NATALIA N. OSIPOVA, DSc, A/Professor, osnat75@mail. ru

IRINA M. BYCHKOVA, Research Assistant, [email protected]

KSENIYA S. SVIRIDOVA, Undergraduate Student, [email protected]

Yuri Gagarin State Technical University of Saratov, 77, Politekhnicheskaya Str., 410054, Saratov, Russia

© Осипова Н.Н., Бычкова И.М., Свиридова К.С., 2017

SOIL RATIONALE FOR HEAT EXCHANGERS IN VAPOR PHASE OVERHEATING IN LIQUEFIED PETROLEUM GAS

The paper proposes a gas supply scheme based on liquefied petroleum gas with natural re-gasification equipped with a soil heat exchanger. This scheme provides a free-hydrate reduction of the vapor phase before its supply to gas-based units. The heat exchange process of the soil heat exchanger is studied in this paper. The proposed soil heat exchanger will provide overheating of the vapor phase of liquefied petroleum gas and avoid the formation of hydrates in pressure-sensitive detectors during vapor throttling.

Keywords: liquefied petroleum gas; soil heat exchanger; natural regasification; vapor phase; overheating; hydrate formation.

Сжиженные углеводородные газы всегда содержат в себе влагу в растворенном состоянии [1]. При дросселировании влажных паров сжиженного газа в регуляторах давления часто образуются гидраты. Образование гидратов негативно сказывается на работе системы газоснабжения, вызывает перебои в подаче газа, требует проведения демонтажа регуляторов давления и их очистки от гидратов. Как показали исследования, на процесс образования гидратов оказывают влияние параметры дросселируемой паровой смеси пропан-бутана: температура, давление, степень сухости, компонентное содержание [2, 3]. На практике вопрос подготовки паровой фазы к дросселированию для систем с искусственной регазификацией решается достаточно просто - производится перегрев паровой фазы в испарителях сжиженного газа до температуры безгидратного редуцирования. В системах с естественной регазификацией сжиженного углеводородного газа предусматриваются различные технические решения для предупреждения гидратообразования [4, 5], в том числе и с перегревом паровой фазы [6, 7]. Одним из способов обеспечения безгидратной эксплуатации систем снабжения сжиженным газом является использование грунтовых теплообменников для перегрева паровой фазы перед ее редуцированием [8].

Использование грунтовых теплообменников формирует систему газоснабжения с естественной регазификацией сжиженного углеводородного газа и исключает использование дополнительных энергоресурсов (электричество, пар и т. д) для испарения газа. Процессы теплообмена резервуаров сжиженного углеводородного газа (сферической изотермической полости) с грунтовым массивом представлены в работе [9].

Рассмотрение процесса теплообмена шахтных грунтовых теплообменников, представляющих собой ряд вертикально размещаемых в грунте труб с грунтовым массивом, подробно изложено в научной работе [8]. В то же время процесс теплообмена горизонтально проложенных грунтовых теплообменников с грунтовым массивом, представляющих собой протяженные цилиндрические трубопроводы, размещаемые вокруг резервуара по дну котлована, изучен недостаточно полно. Для определения механизма перегрева паровой фазы сжиженного углеводородного газа в грунтовом теплообменнике необходимо провести исследование процесса теплообмена.

Схема системы газоснабжения с естественной регазификацией сжиженного углеводородного газа с перегревом паров в грунтовом теплообменнике

с последующим редуцированием паровой фазы в регуляторе давления представлена на рисунке.

Схема системы газоснабжения с перегревом паров в грунтовом теплообменнике:

1 - подземный резервуар; 2 - трубопровод паровой фазы СУГ; 3 - трубчатый грунтовый теплообменник; 4 - изолированное укрытие редуцирующего узла; 5 -регулятор давления; 6 - газопровод низкого давления; 7 - тепловая изоляция; 8 -кожух резервуарной головки

Схема газоснабжения с естественной регазификацией работает следующим образом. Сжиженный углеводородный газ за счет теплообмена подземного резервуара 1 с грунтом из жидкого состояния переходит в газообразное и по трубопроводу паровой фазы 2 подается в трубчатый грунтовый теплообменник 3. За счет естественного притока тепла из грунта пары сжиженного углеводородного газа, продвигаясь по трубопроводу грунтового теплообменника, повышают свою температуру. Перегретая паровая фаза поступает в регулятор давления 5, размещаемый в изолированном укрытии 4. Для сохранения перегрева паровой фазы, полученной за счет теплообмена в горловине резервуара и грунтовом теплообменнике, трубопроводы системы газоснабжения, сообщающиеся с наружным воздухом и промерзающим грунтом, покрываются тепловой изоляцией 7. Перегрев паровой фазы позволяет провести редуцирование в регуляторе давления без образования гидратов. После редуцирования паровая фаза низкого давления подается к потребителю в газоиспользующие приборы по газопроводу 6.

При постановке задачи исследований учитывалось влияние температурного поля подземного резервуара сжиженного газа на температурные режимы эксплуатации грунтового теплообменника. В то же время влияние температурного поля грунтового теплообменника на температурное поле самого резервуара не учитывалось, поскольку диаметры резервуара и трубопровода теплообменника определяются отношением Д,/а?то > 20, что формирует второй порядок малости влияния.

Для решения поставленной задачи воспользуемся методом тепловых источников и принципом наложения температурных полей (суперпозиции) [10]. Представим создаваемое температурное поле в грунтовом массиве указанной системы газоснабжения как сочетание собственного температурного поля грунта, формируемого по глубине массива и создающего распределение температур в районе заложения грунтового теплообменника Цх) и температурного поля, возбуждаемого в массиве с нулевой температурой поверхности линейным источником теплоты (грунтовым теплообменником) переменной интенсивности д(х).

Представим трубопровод грунтового теплообменника длиной I в виде последовательно чередующихся отрезков i = 1, 2, 3,... п протяженностью Ах. В пределах каждого отрезка интенсивность линейного источника примем постоянной и равной дг-(х). В точке А с координатами хА и уА грунтового массива воздействие элемента dx линейного источника теплоты дг-(х) формирует температуру dí, определяемую по выражению

=

д1 (х)йх

4 %к

I

I

^¡(х^Ха^ТУА у1(2х - хА)2 + у^2

(1)

где х, у - координаты элемента линейного источника теплоты; хА, уА - координаты точки А в грунтовом массиве; X - коэффициент теплопроводности грунтового массива Вт/м-к; д, - интенсивность линейного источника теплоты, Вт/м.

Линейный источник теплоты в пределах каждого отрезка эквивалентен точечному источнику интенсивностью Qi = qiАх. Тогда температура массива в точке А, формируемая под воздействием всех точечных элементов грунтового теплообменника, определится по выражению

1 п t / =—Те,

4%Х

;=1

I

I

(2)

Естественная температура массива на глубине заложения трубопровода грунтового теплообменника

^ (ХА ). (3)

Формируемое температурное поле вокруг трубопровода грунтового теплообменника определится выражением

^ = -

I

4%\

I

I

л[(х—ха)Г+УА у1(2х - хА )

2 2 2 + уА

+ (е ( ха ) .

(4)

1=1

Граничные условия формируемой задачи:

- на поверхности грунтового массива при у ^ 0; = ¿н;

- в грунтовом массиве при у ^ да; х ^ да: ?(х) = ¿е(х).

Для решения поставленной задачи с учетом ограничивающих условий воспользуемся приближенным численным методом (методом последовательных приращений).

При разбиении общей длины грунтового теплообменника I на п отрезков протяженностью Ах температура грунта ^ на поверхности расчетного отрезка трубопровода Ах в т. Ы, обусловленная воздействием температурных полей, определится

1 п гр 4%ХУ е

I I

л/(х _ )2 + у2м 7(2_ хN)2 + х

+ (х„ ) , (5)

где Qi - интенсивность 7-го источника (стока) тепла, Вт; х7 - координата 7-го источника (стока) тепла, м; хЫ, уЫ - координаты точки Ы, лежащей на середине п-го отрезка трубопровода грунтового теплообменника, м; ¿е(хп) - естественная температура грунта на глубине заложения трубопровода грунтового теплообменника, м.

Уравнение теплового баланса любого отрезка п грунтового теплообменника имеет вид

сG ^ _1 ) = а Ах, (6)

п \ кон нач/ 1п ' V /

где сп - массовая теплоемкость паровой фазы Вт-ч/(кг-К); О - расход сжиженного углеводородного газа, кг/ч; *нач; *кон - температура паровой фазы СУГ в начале и в конце п-го отрезка грунтового теплообменника, °С; Ах - протяженность п-го отрезка грунтового теплообменника, м; дп - удельный приток тепла к паровой фазе СУГ на п-м отрезке грунтового теплообменника, Вт/м.

В свою очередь, удельный приток тепла от грунтового массива к паровой фазе СУГ на любом п-м участке грунтового теплообменника определяется уравнением теплопередачи цилиндрической стенки с учетом ламинарного течения среды при естественной регазификации сжиженного углеводородного газа

Ч„„ + („

2%

Чп =

* _| " нач кон

гр,п

I , йиз 2

-1п—^ +-

X йнар а йвн

из тр п тр

(7)

где Хиз - коэффициент теплопроводности гидроизоляции грунтового теплообменника, Вт/ (м-К); а?из - наружный диаметр грунтового теплообменника с учетом гидроизоляции, м; й™р - наружный диаметр трубопровода грунтового теплообменника, м; й™ - внутренний диаметр трубопровода грунтового

теплообменника, м; ап - коэффициент теплоотдачи от внутренней поверхности трубопровода к паровой фазе СУГ, Вт/(м2-К).

Поскольку в пределах расчетного отрезка Ах изменение температуры паровой фазы СУГ весьма незначительно, положим в качестве допущения

г =г . (8)

нач кон V /

Тогда уравнение (7) примет следующий вид:

2лГг - г 1

_ |_ Гр,и нач J („)

= ± + -1_. ()

X с1нар а с1вн

из тр п тр

Определение коэффициента теплоотдачи ащ Вт/(м2К) проводилось в соответствии со следующими исходными данными:

- грунтовый теплообменник выполнен из стальной бесшовной трубы по ГОСТ 8732-78 диаметром 28x3 мм, внутренний диаметр указанной трубы составляет ^вн = 0,022 м;

- расчетный расход газа О = 1,2 кг/ч (0 = 0,6 м3/ч) (газоснабжение индивидуального жилого здания коттеджного типа при наличии газовой четы-рехконфорочной плиты и двухконтурного газового котла для отопления и горячего водоснабжения с тепловой нагрузкой 10 кВт).

При заданном диаметре трубопровода и расходе паровой фазы СУГ критерий Рейнольдса составляет Яе = 2412, что соответствует турбулентному режиму течения.

Для вычисления коэффициента теплоотдачи в условиях вынужденного стабилизированного турбулентного потока в трубах и каналах используется формула [11]

Ыи = 0,023Рг0'4 Яе0'8, (10)

где Рг - критерий Прандтля.

В результате преобразования выражения (10) имеем

( ^ У'8

9,96 • 10-5 Рг04 X

Рв

а =-ч ^ 'V , (11)

п , ( )

где X - теплопроводность паровой фазы СУГ, Вт/(м-К); р - плотность паровой фазы пропан-бутановой смеси, кг/м3; ^ - динамическая вязкость паровой фазы сжиженного углеводородного газа, Па-с.

Система уравнений (1) - (11) представляет собой математическую модель теплообмена грунтового теплообменника сжиженного газа с грунтовым массивом.

Реализация модели производится методом последовательных приращений с учетом разбиения грунтового теплообменника на п участков при условии п > 10. Начальная температура для первого участка грунтового теплообменника г1нач - температура паровой фазы СУГ на выходе из резервуара Конечная температура для первого участка г1кон и т. д., конечная температура на выходе из грунтового теплообменника для последнего участка - ^

,кон 1ш •

В качестве практического применения рассматриваемой математической модели определим степень перегрева паров в пароперегревателе при следующих исходных данных:

- климатический район эксплуатации резервуарной установки г. Екатеринбург;

- расчетная температура наружного воздуха tB = -31 °С;

- геометрический объем подземного резервуара сжиженного газа V = 3,84 м3;

- расход газа G = 1,2 кг/ч;

- массовая теплоемкость паровой фазы сп = 0,47 Вт-ч/(кг-К);

- конструкция грунтового теплообменника - стальная труба диаметром 28x3 мм с полимерной гидроизоляцией толщиной 2 мм и теплопроводностью

= 0,174 Вт/(мК);

- трубопровод грунтового теплообменника уложен по дну котлована ре-зервуарной установки вокруг резервуара (глубина заложения Нп = 2,9 м) на расстоянии от него 1 м. С учетом того, что равноудаленное расстояние от трубопровода паровой фазы любого точечного источника грунтового теплообменника от резервуара определяется расстоянием, которое может быть представлено эквивалентным радиусом, можно записать

lDp

R =-р— +1,0 м. (12)

Vb ^ +1 , ( )

Уравнение теплового баланса для участка n грунтового теплообменника длиной dln имеет следующий вид:

qdl (t -t) = cGdt , (13)

in n^ гр n' n '

где dtn - приращение температуры паровой фазы на элементе грунтового теплообменника длиной dln; tIV - температура на поверхности гидроизоляции трубопровода грунтового теплообменника, принимая равной температуре окружающего грунтового массива, °С; qn - удельный теплоприток к участку грунтового теплообменника, Вт/м2, определяемый по формуле (9).

Значение температуры на поверхности гидроизоляции = const определяется в соответствии с уравнением (5). При этом координаты k-й точки имеют следующий вид:

ID

Хк=Т^Б~+1,0; Ук =ндоп + К, (14)

где hH - глубина заложения грунтового теплообменника, м; Ндоп - высота дополнительного слоя грунта, принимаемого равным высоте снежного покрова.

Разделим в уравнении (13) переменные

cGdt

dl=-^— (15)

n q (t -1 ) ( )

гр,- -,Пs

и проинтегрируем полученное выражение:

idln = - ir^". C6)

0 qn i„ *гр,п tn,u

где 1п - длина грунтового теплообменника, м; 1п - температура паровой фазы СУГ на выходе из резервуара, °С; - температура паровой фазы СУГ в конце грунтового теплообменника, °С. В результате имеем:

I = 1п ~?п , (17)

п а ? - ?кон

~п гр,п п,п

откуда температура паровой фазы СУГ в конце грунтового теплообменника

Г>н = - + г . (18)

Чп г а \ т 'гр,и у '

ехр1

Результаты численной реализации математической модели (1) - (11) с учетом выражений (14) - (18) представлены в таблице.

Формирование температур в системе газоснабжения с естественной регазификацией

Температура (°С) грунта на глубине заложения грунтового теплообменника Температура (°С) паровой фазы СУГ в элементах резервуарной установки

в естественном состоянии с учетом теплового воздействия резервуара подземный резервуар грунтовый теплообменник

на входе на выходе

+2,03 +0,62 -8,63 -8,63 +0,52

Как видно из таблицы, прокладка грунтового теплообменника вокруг резервуара по дну котлована обеспечивает перегрев паровой фазы СУГ с -8,63 до +0,52 °С или на 9,15 °С.

Согласно исследованиям [2, 3], установлено, что область образования гидратов пропан-бутановой смеси при давлении от 0,5 МПа и ниже находится в температурном интервале от -20 до -5 °С.

Таким образом, проведенными исследованиями установлено, что применение грунтового теплообменника позволяет обеспечить подачу в редуцирующий узел перегретой паровой фазы и провести процесс снижения давления в регуляторе в режиме, исключающем кристаллизацию влаги в дросселирующем органе.

Библиографический список

1. Системный анализ возникновения источников свободной воды и ее накопления в подземных резервуарных установках сжиженного углеводородного газа / А.П. Усачев, А.Л. Шурайц, Л.Л. Феоктистов, М.С. Недлин, А.В. Рулев, Т.А. Усачева // Нефтегазовое дело: транспорт и хранение. - 2009. - Т. 7. - № 1. - С. 98-101.

2. Осипова, Н.Н. Исследование процесса гидратообразования при редуцировании влажного газа / Н.Н. Осипова // Приволжский научный журнал. - 2012. - № 3 (23). - С. 112-117.

3. Максимов, С.А. Предупреждение гидратообразования в резервуарных установках с естественной регазификацией сжиженного газа / С.А. Максимов // Научная дискуссия: инновации в современном мире : материалы IX Международной заочной конференции. -М., 2013. - С. 177-180.

4. Шурайц, А.Л. Расчет экономии энергетических ресурсов на нужды регазификации сжиженного углеводородного газа за счет использования природной теплоты грунта / А.Л. Шурайц, А.П. Усачев, А.В. Рулев // Проблемы и методы обеспечения надежности и безопасности систем транспорта нефти, нефтепродуктов и газа : материалы Международной научно-практической конференции в рамках нефтегазового форума и XXII Международной специализированной выставки «Газ. Нефть. Технологии - 2014». - Уфа, 2014. - С. 123-131.

5. Осипова, Н.Н. Объективный выбор толщины тепловой изоляции участков трубопроводной обвязки узла редуцирования с целью предупреждения гидратообразования / Н.Н. Осипова, Б.Н. Курицын, С.А. Максимов // Вестник МГСУ. - 2011. - № 7. - С. 525-530.

6. Курицын, Б.Н. Использование природного тепла грунта для обогрева регуляторов давления сжиженного газа / Б.Н. Курицын, С.А. Максимов // Разработка современных технологий и материалов для обеспечения энергосбережения, надежности и безопасности объектов архитектурно-строительного и дорожного комплекса : сб. статей Международного научно-практического симпозиума «Социально-экономические проблемы жилищного строительства и пути их решения в период выхода из кризиса». - Саратов, 2010. - С. 96-101.

7. Курицын, Б.Н. Разработка и обоснование технических решений по предупреждению гидратообразования в системах резервуарного снабжения сжиженным газом / Б.Н. Курицын, Н.Н. Осипова, С.А. Максимов // Приволжский научный журнал. - 2013. -№ 1 (25). - С. 73-80.

8. Курицын, Б.Н. Применение грунтовых теплообменников в системах регазификации сжиженных углеводородных газов/ Б.Н. Курицын, Г.П. Чирчинская, М.В. Павлутин // Газ России. - 2005. - № 1. - С. 31-33.

9. Курицын, Б.Н. Моделирование теплообмена при хранении сжиженного газа в подземных резервуарных установках под воздействием естественных температур грунта и наружного воздуха / Б.Н. Курицын, Н.Н. Осипова, С.А. Максимов // Научный вестник Воронежского государственного архитектурно-строительного университета. Строительство и архитектура. - 2012. - Вып. 2 (26). - С. 35-46.

10. Пентегов, И.В. К теории метода тепловых источников, используемого при анализе тепловых процессов в электротехнических системах / И.В. Пентегов // Электротехнические и информационные комплексы и системы. - 2014. - Т. 10. - № 3. - С. 5-15.

11. Михеев, М.А. Основы теплопередачи / М.А. Михеев, И.М. Михеева. - М. : Энергия, 1977. - 344 с.

References

1. Shuraits A.L., Usachev A.P., Feoktistov A.A., Nedlin M.S., Rulev A.V., Usacheva T.A. Sistemnyi analiz vozniknoveniya istochnikov svobodnoi vody i ee nakopleniya v podzemnykh rezervuar-nykh ustanovkakh szhizhennogo uglevodorodnogo gaza [Systems analysis of the emergence of free water sources and its accumulation in underground tank installations for liquefied petroleum gas]. Neftegazovoe delo [Oil and Gas Business]. 2009. V. 7. No. 1. Pp. 98-101.(rus)

2. Osipova N.N. Issledovanie processa gidratoobrazovaniya pri reducirovanii vlazhnogo gaza [A study of hydrate formation process at wet gas reduction]. Privolzhsky Scientific Journal. 2012. No. 3 (23). Pp. 112-117. (rus)

3. Maksimov S.A. Preduprezhdenie gidratoobrazovaniya v rezervuarnykh ustanovkakh s estestvennoi regazifikatsiei szhizhennogo gaza [Prevention of hydrate formation in tanks with natural liquid gas regasification]. Proc. 9th Int. Conf. 'Innovation in the Modern World'. Moscow, 2013. Pp. 177-180. (rus)

4. Shuraits A.L., Usachev A.P., Rulev A.V. Raschet ekonomii energeticheskikh resursov na nu-zhdy regazifikatsii szhizhennogo uglevodorodnogo gaza za schet ispol'zovaniya prirodnoi tep-loty grunta [Energy savings for regasification of liquefied petroleum gas due to natural soil heating]. Proc. Int. Sci. Conf. 'Problems and Methods of Reliable and Safe oil Transportation Systems and Products'. Ufa, 2014. Pp. 123-131. (rus)

5. Osipova N.N., Kuritsyn B.N., Maksimov S.A. Ob"ektivnyi vybor tolshchiny teplovoi izolyatsii uchastkov truboprovodnoi obvyazki uzla redutsirovaniya s tsel'yu preduprezhdeniya gidra-

toobrazovaniya [Selection of thermal insulation thickness of pipeline portions with reduction assembly to prevent hydrate formation]. VestnikMGSU [Scientific and Technical Journal on Construction and Architecture]. 2011. No. 7. Pp. 525-530. (rus)

6. Kuritsyn B.N., Maksimov S.A. /spol'zovanie prirodnogo tepla grunta dlya obogreva regulya-torov davleniya szhizhennogo gaza [Natural soil heating for pressure-sensitive detector heating of liquefied gas]. Coll. Papers Int. Sci. Symp. 'Socioeconomic Housing Problems and Their Solution During Crisis Recovery'. Saratov, 2010. Pp. 96-101. (rus)

7. Kuritsyn B.N., Osipova N.N., Maksimov S.A. Razrabotka i obosnovanie tekhnicheskikh reshenii po preduprezhdeniyu gidratoobrazovaniya v sistemakh rezervuarnogo snabzheniya szhizhennym gazom [Development and support of technical solutions on prevention of hydrate formation in liquefied gas systems]. Privolzhsky Scientific Journal. 2013. No. 1 (25). Pp. 73-80. (rus)

8. Kuritsyn B.N., Chirchinskaya G.P., Pavlutin M.V. Primenenie gruntovykh teploobmennikov v sistemakh regazifikatsii szhizhennykh uglevodorodnykh gazov [The use of soil heat exchangers in regasification systems of liquefied petroleum gas]. Gaz Rossii. 2005. No. 1. Pp. 31-33. (rus)

9. Kuritsyn B.N., Osipova N.N., Maksimov S.A. Modelirovanie teploobmena pri hranenii szhizhennogo gaza v podzemnyh rezervuarnyh ustanovkah pod vozdejstviem estestvennyh temperatur grunta i naruzhnogo vozduha [Modeling of heat transfer in the storage of liquefied natural gas in underground tank installations under the influence of the natural rate of soil tires and the outdoor]. Scientific Bulletin of Voronezh State Architecture and Construction University. Construction and architecture. 2012. V. 2 (26). Pp. 35-46. (rus)

10. Pentegov I.V. K teorii metoda teplovykh istochnikov, ispol'zuemogo pri analize teplovykh protsessov v elektrotekhnicheskikh sistemakh [The analysis Heat sources method in analyzing thermal processes in electrical systems]. Elektrotekhnicheskie i informatsionnye kompleksy i sistemy. 2014. V. 10. No. 3. Pp. 5-15. (rus)

11. Mikheev M.A., Mikheeva I.M. Osnovy teploperedachi [Fundamentals of heat transfer]. Moscow: Energiya Publ., 1977. 344 p. (rus)