Mathematical model of the condensation process in a cylinder of a piston engine
- Authors: Kolunin A.V.1, Lazarev E.S.2, Kaminskiy V.N.1, Korytov M.S.3, Ruzimov A.O.4
-
Affiliations:
- Moscow Polytechnic University
- Omsk State Transport University
- Siberian State Automobile and Highway University
- Ministry of Defense of the Russian Federation
- Issue: Vol 90, No 5 (2023)
- Pages: 395-404
- Section: Environmentally friendly technologies and equipment
- URL: https://journals.eco-vector.com/0321-4443/article/view/607362
- DOI: https://doi.org/10.17816/0321-4443-607362
- ID: 607362
Cite item
Abstract
BACKGROUND: In recent years, there has been a trend of increasing activity towards the development of polar territories. A characteristic feature of the North is negative ambient temperatures that have a negative impact on the condition of piston engines of ground transport, mobile and stationary power plants and labor saving tools. An engine is the least adapted unit for use in such conditions. There is a chain of negative factors that consistently links negative ambient temperatures, in which the equipment is operated, and the condition of the mechanisms and engine systems. The primary link of this chain is condensation processes. The existence of condensation processes during low-temperature operation of the engine has been experimentally proved. The latter takes place when warming up in conditions of negative ambient temperatures. The question «How much water changes the state during the warm-up period?» arises.
AIMS: Development of a mathematical model that makes possible to obtain unbiased information about the activity of condensation processes and to estimate the amount of water that changes the state during the warm-up period.
METHODS: Solving the given tasks is based on classical theories describing operational processes of boilers. The high labor intensity and significant financial costs in organizing such experiments require the search for new research methods. Mathematical models help to solve the task of defining the mass amount of water condensing in a cylinder of a piston engine computationally.
RESULTS: The mathematical model that is characterized by its adaptation to piston engines and is capable of determining the mass amount of water changing the state during the warm-up period iteratively, using the differences in partial pressures and the density of the mass flow of water condensate, has been developed.
CONCLUSIONS: The existence of water has a negative impact on conditions of a piston engine. The information about the amount of water condensing in a cylinder during the warm-up period stimulates to continue studies in the field of motor oils watering, active acids formation and corrosive wear of surfaces of details.
Full Text
ВВЕДЕНИЕ
На рубеже XX–XXI веков правительства ряда государств обратили свой взор на заполярные территории нашей планеты. Интерес обосновывается комплексом аргументов. Несметные, не тронутые до настоящего времени запасы полезных ископаемых. Большое геополитическое значение, заключающееся в объединении воедино окраины материков Евразии и Северной Америки, Северный Ледовитый океан с островами, а также прилегающие части Атлантического и Тихого океанов. Кроме того, Северный морской путь — кратчайшая коммуникация между Европейской частью России и Дальним Востоком; Северо-Западный проход — морской путь между Атлантическим и Тихим океанами и воздушный мост между Северной Америкой и Юго-Восточной Азией. Развитие и освоение Арктики может осуществляться на основе применения передовых технологий и средств наземного транспорта приспособленных к суровым условиям заполярных территорий. Именно наземные средства принимают на себя большую часть задач по обеспечению перевозок и выполнению работ технологического характера. Особенностью этих мест являются отрицательные температуры, сокращающие жизненный цикл техники.
Анализ работ и статистических данных выявил один из наименее приспособленных агрегатов к таким условиям, коим является поршневой двигатель. Низкая приспособленность объясняется работой в широких температурных и нагрузочных диапазонах, с химически активными соединениями. При низкотемпературных режимах повышается хрупкость конструкционных и вязкость эксплуатационных материалов, ухудшаются условия трения сопряжённых поверхностей деталей, снижается полнота сгорания топлива, повышая активность образования химически активных соединений. Снижается надёжность работы, сокращается жизненный цикл.
Одной из причин, оказывающих негативное влияние на состояние двигателя, является наличие процессов конденсации водяных паров на поверхностях деталей и масла [1].
Вода образуется по причине окисления водорода топлива атомами кислорода. По разным источникам выхлопные газы поршневого двигателя содержат водяные пары порядка 8–12% по объёму. Водяные пары, содержащиеся в рабочем теле двигателя, в выхлопных газах, в картерных газах, в условиях оптимального температурного режима без особого вреда выводятся через выпускные и вентиляционные системы. При этом следует учитывать то, что выходу двигателя на такой режим может предшествовать низкотемпературное состояние. Изначально, при окислении углеводородов топлива в камере сгорания выделяется большое количество тепла, что приводит к разогреву продуктов сгорания до 1500–2000 ºС и подъёму давления в цилиндре до 6–10 МПа. Такие параметры обеспечивают высокое давление насыщения, и газообразное состояние воды в смеси газов. Однако, холодные детали и масло снижают давление насыщения вблизи своих поверхностей. При контакте с поверхностью стенки цилиндра, имеющей температуру ниже температуры насыщения при текущем парциальном давлении водяных паров, возникает конденсация. Активность конденсационного процесса, характеризуемого потоком массы, не является стабильной. Значительные перепады температур и давлений определяют высокие дестабилизационные условия. При прогреве двигателя температура стенки цилиндра возрастает и с момента достижения точки росы конденсация прекращается.
Актуальность работы заключается в негативном влиянии воды на состояние двигателя. Вода неизбежно поступает в картерное пространство через сопряжение деталей цилиндропоршневой группы, растворяется в моторном масле, нарушая коллоидную стабильность, инициирует усиление межмолекулярных взаимодействий продуктов низкой агрегатной устойчивости [2]. Незначительное обводнение моторного масла может привести к деактивации пакета присадок, введённого в базовое масло при его производстве. При многократном пуске двигателя в условиях отрицательных температур такое воздействие имеет накопительный характер.
Не следует сбрасывать со счетов влияние воды на конструкционные материалы. Являясь коррозионно-активным соединением, она инициирует процессы коррозионного изнашивания поверхностей деталей. В значительной степени такому изнашиванию подвержены детали цилиндро-поршневой группы.
Постановка задачи
Эти процессы могут иметь место в механизмах и системах двигателя. Необходимо систематизировать области конденсационных процессов. К первичным областям следует отнести пространства, в которых происходит непосредственно образование воды по причине окисления атомов водорода топлива. А именно, камера сгорания и внутренняя поверхность гильзы цилиндра. Ко вторичным – внутренние повехности деталей и сборочных единиц обеспечивающих отвод газов. И наконец, поверхности в объёме картерного пространства представляют третью область таких процессов [3]. Определение расчётного массообразования воды представляет непростую задачу, требующую системного подхода с учётом множества факторов влияющих на активность процессов. В соответствии с представленной системой областей, необходимо разработать математическую модель, описывающую конденсационные процессы в цилиндре двигателя. Математическая модель, перспективно обеспечивающая решение задачи по определению количества воды, конденсирующейся в цилиндре за период прогрева. Создание самой модели обуславливается такими факторами как: значительные перепады температур и давлений, турбулизация газов, изменение направлений и скоростей движения поршня, активность теплообменных процессов.
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ
При разработке математической модели приняты допущения, не оказывающие существенного влияния на погрешность расчётов.
Допущение 1. Толщина металлической стенки цилиндра и теплоотдача от охлаждающей жидкости не оказывают термического сопротивления. Масляная плёнка и водный конденсат имеют ничтожно малую толщину слоёв. Для упрощения алгоритма математической модели можно пренебречь термическим сопротивлением охлаждающей жидкости, стенки цилиндра, масляной плёнки и водного конденсата. Температуры стенки цилиндра и камеры сгорания принять равными температуре охлаждающей жидкости.
Допущение 2. При движении газа вблизи стенки цилиндра за счёт действия сил вязкости образуется слой замедленного движения вещества [4]. Этот же слой оказывает сопротивление свободному перемещению молекул газов, в том числе и водяных паров. В силу высокой турбулизации, вызванной движением поршня, концентрацию водяных паров на границе пограничного слоя считать равной концентрации во всём объёме цилиндра.
Допущение 3. Парциальное давление на стенке цилиндра приравнивается к давлению насыщения при температуре насыщения водяных паров, равной температуре охлаждающей жидкости.
Допущение 4. Возвратно-поступательное движение поршня обеспечивает высокую турбулизацию газов в цилиндре, интенсивное перемешивание и высокую равномерность концентрации компонентов в газовой смеси. Поэтому концентрацию водяных паров считать равномерно распределёнными в объёме цилиндра, кроме пристеночного слоя.
Допущение 5. Поскольку продолжительность полного рабочего цикла ничтожно мала относительно времени прогрева охлаждающей жидкости и количество передаваемого от рабочего тела тепла к стенке цилиндра за один цикл не будет производить значительного нагрева охлаждающей жидкости и масла, то температуру стенки цилиндра в рамках одного цикла считать постоянной.
Допущение 6. Влагосодержание топлива и воздуха не принимать в расчёт по причине низких абсолютных значений.
В качестве исходных данных принимаются следующие показатели:
- Индикаторные характеристики двигателя;
- Элементный состав топлива;
- Ход поршня, S, м;
- Диаметр цилиндра, d, м;
- Радиус кривошипа, R, м;
- Площадь камеры сгорания, м2;
- Угол поворота коленчатого вала, φ, °п.к.в.;
- Коэффициент избытка воздуха, α;
- Уравнение регрессии, описывающее изменение температуры охлаждающей жидкости;
- Частота вращения коленчатого вала, n, об/мин.
Массовое количество конденсата, образующегося на внутренней поверхности гильзы цилиндра, зависит от активности конденсационных процессов и продолжительности периода определяемого динамикой приращения температуры.
Каждый рабочий цикл характеризуется образованием определённого количества водного конденсата. Зная количество прошедших за расчётное время циклов при конкретной температуре охлаждающей жидкости и закон изменения температуры охлаждающей жидкости, интегрированием по времени можно получить общую массу конденсата:
, (1)
где Jцикл — поток массы за цикл, кг/цикл; τпрогр — время нагрева охлаждающей жидкости до рабочего уровня, с.
Поток массы компонента смеси [5]:
, (2)
где J — поток массы, кг/с; — масса диффундирующего компонента смеси, кг; τ — время, с; j — плотность потока массы, кг/(с·м2); F — площадь конденсации, м2.
Масса конденсата, выпадающего на стенках цилиндра, будет суммироваться по продолжительности каждого такта по углу поворота коленчатого вала. За нулевое значение угла принимается положение ВМТ перед тактом впуск. Разделив весь цикл на 4 отдельных такта, преобразуем выражение (1) к виду:
(3)
Каждое из четырёх слагаемых в правой части уравнения (3) соответствует отдельному такту.
Так как влагосодержание воздуха при отрицательных температурах ничтожно мало, то первые два слагаемых уравнения (3) можно отбросить как несущественные. Тогда среднее за цикл значение потока массы определится следующим выражением:
, (4)
где Jцикл — поток массы за цикл, кг/цикл.
Для решения задачи расчёта массы выпадающего конденсата на стенке цилиндра рассмотрим кольцевой участок стенки с высотой определённой перемещением поршня S. Площадь этого участка определится из произведения трёх сомножителей:
, (5)
где F — площадь элементарного цилиндрического участка, м2; d — диаметр цилиндра, м; S — перемещение поршня, м
Движение поршня вдоль оси цилиндра подчиняется косинусоидальному закону, таким образом, можно выразить высоту рассматриваемого участка следующей зависимостью [6]:
, (6)
где R — радиус кривошипа, м; λ — отношение радиуса кривошипа к длине шатуна, ; φ — угол поворота коленчатого вала в градусах, °п.к.в.
Тогда площадь внутренних поверхностей цилиндра, на которых осуществляется процесс конденсации, определяется суммой рабочей поверхности и поверхности камеры сгорания:
, (7)
где Fкс — площадь поверхностей камеры сгорания, м2.
Плотность потока массы может быть определена с помощью выражения [5]:
, (8)
где ρ — плотность смеси газов, кг/м3; β — коэффициент массоотдачи, м/с; mп0 — массовая доля концентрации пара в основном объёме цилиндра; mп.гр — массовая доля концентрации пара на поверхности конденсации.
Плотность газов будет меняться в зависимости от температуры, давления в цилиндре и состава.
Зависимость давления и температуры в цилиндре от угла поворота коленчатого вала выражена в индикаторной диаграмме конкретной модели двигателя.
Плотность продуктов сгорания (ρпс, кг/м3) согласно источника [7] и уравнения состояния идеального газа представляется следующей зависимостью:
, (9)
где — теоретический объём азота в продуктах сгорания, м3/кгтопл; — объём трёхатомных газов в продуктах сгорания, м3/кгтопл; — объём водяных паров в продуктах сгорания, м3/кгтопл; — количество избыточного воздуха для сгорания 1 кг топлива, м3/кгтопл; Т — температура, К; p — давление, Па
Значения объёмов компонентов могут быть определены согласно формул [7]:
- теоретический объем азота в продуктах сгорания, м3/кг:
(10)
здесь V0 — теоретический объем сухого воздуха, м3/кг, необходимого для полного сгорания топлива (при коэффициенте избытка воздуха α = 1) описывается как:
(11)
где C — содержание углерода в топливе, %; S — содержание серы в топливе, %; H — содержание водорода в топливе, %; O — содержание кислорода в топливе, %;
- избыточный объем воздуха, м3/кг:
(12)
- теоретический объем трехатомных газов в продуктах сгорания, м3/кг:
(13)
- теоретический объем водяных паров в продуктах сгорания, м3/кг:
(14)
где W — содержание воды в топливе, %; dвл — влагосодержание воздуха, г/кг сухих газов.
В качестве примера, значения объёмов продуктов сгорания дизельного топлива для коэффициента избытка воздуха α = 1,7 сведены в таблицу 1.
Таблица 1. Значения объёмов продуктов сгорания на 1 кг топлива
Table 1. Values of combustion products volume per 1 kg of fuel
, |
|
|
|
8,67 | 1,61 | 1,40 | 7,68 |
, (15)
где — объём компонента смеси, м3; – массовая доля компонента смеси:
, (16)
где — молярная масса компонента смеси, кг/моль; — объёмная доля компонента смеси газов:
, (17)
где — давление парогазовой смеси, Па:
. (18)
По данным таблицы 1 произведён расчёт масс компонентов и сведён в таблицу 2.
Таблица 2. Значения масс продуктов сгорания на 1 кг топлива
Table 2. Values of combustion products mass per 1 kg of fuel
|
|
| |
10,87 | 3,17 | 1,13 | 9,97 |
Для определения коэффициента массоотдачи применены критериальные уравнения, полученные на основе тройной аналогии подобия процессов теплоотдачи, массоотдачи [5]:
, (19)
где — диффузионный критерий Нуссельта; d — характеристический размер (в рассматриваемой задаче диаметр цилиндра), м; D — коэффициент диффузии, м2/с.
Коэффициент диффузии (D) [8]:
, (20)
где — коэффициент диффузии при нормальных условиях, м2/с; Т — температура, К; = 273К; p — давление, Па; = 101325 Па.
Ввиду схожести свойств воздуха и продуктов сгорания принято: и u = 0,8 [4]. Тогда:
, (21)
где — диффузионный критерий Нуссельта; — диффузионный критерий Нуссельта, полученный на основании тройной аналогии; — объёмная доля сухих газов; — показатель характеризующий разность парциальных давлений.
Для случая течения газа в трубе, что наиболее близко к случаю течения продуктов сгорания в цилиндре поршневого двигателя диффузионный критерий Нуссельта, полученный на основании тройной аналогии, определится выражением [5]:
, (22)
где Re — критерий Рейнольдса; — диффузионный критерий Прандтля; — объёмная доля водяных паров.
При этом критерий Рейнольдса:
, (23)
где — скорость движения среды, м/с; — вязкость смеси, м2/с.
Скорость движения среды может быть приравнена к скорости движения поршня [5]:
, (24)
где: n — частота вращения коленчатого вала, об/мин; ω — угловая скорость вращения коленчатого вала, рад/с:
. (25)
При коэффициенте избытка воздуха α = 1,7 состав продуктов сгорания близок к среднему [7]. Вязкость продуктов сгорания аппроксимируются выражением:
, (26)
где — температура продуктов сгорания, °С.
Диффузионный критерий Прандтля:
, (27)
, (28)
где — парциальное давление сухих газов, Па.
, (29)
где — парциальное давление пара вдали от стенки, Па; — парциальное давление пара на границе раздела фаз, Па.
. (30)
Учитывая закон Амага и Авогадро, можно получить:
, (31)
. (32)
Согласно принятому допущению № 3, парциальное давление на стенке цилиндра приравнивается к давлению насыщения при температуре насыщения водяных паров, равной температуре охлаждающей жидкости. Уравнение давления насыщения для диапазона температур от -60 до 0°С представляется экспоненциальной формой [9]:
, кПа, (33)
где — температура охлаждающей жидкости, °С.
Для продуктов сгорания дизельного топлива , где — отношение парциального давления сухих газов к общему давлению смеси, тогда согласно [5] коэффициенты в уравнении (21): А = 0,71; k = – 0,9; l = – 0,1.
Описанные зависимости, в определённой последовательности, сведённые в схему, изображённую на рис. 1, в совокупности представляют собой алгоритм математической модели.
Рис. 1. Схема алгоритма математической модели конденсационных процессов в цилиндре поршневого двигателя.
Fig. 1. A diagram of the algorithm of the mathematical model of condensation processes in a cylinder of a piston engine.
Дискуссия
Относительная массовая концентрация паров вдали от поверхности конденсации является величиной переменной, т.к. масса газа в цилиндре ограничена и отсутствует подпитка новым веществом. В результате в каждый следующий момент времени парциальное давление водяных паров уменьшается в зависимости от уже сконденсированной массы воды.
Существенными факторами являются перепады давлений и температур. От температуры и давления зависит коэффициент диффузии. Также давление в цилиндре в свою очередь оказывает влияние на парциальное давление водяных паров. Т.е. парциальное давление водяных паров меняется не только вследствие конденсации, а в большей степени за счёт процессов сгорания и расширения газов.
Также количество уравнений в настоящей модели велико и одна только запись результирующего уравнения является весьма громоздкой.
С применением математической модели построена графическая зависимость массы воды по смене агрегатного состояния, в зависимости от начальной температуры прогреваемого двигателя. Такая зависимость представлена на рис. 2.
Рис. 2. Зависимость массы воды по смене агрегатного состояния, от начальной температуры, за период прогрева двигателя.
Fig. 2. Dependence of the mass of water on the change of the aggregate state, from the initial temperature, during the engine warm-up period of the engine.
В силу перечисленных выше факторов практическая реализация модели выполняется с помощью численного интегрирования с шагом угла Δφ методом трапеций с последовательным расчётом концентраций и парциальных давлений.
Реализация математической модели применительно к двигателю КамАЗ-740.30 рабочим объёмом 10,9 литра, при температуре минус 32 °С, позволило определить расчётное количество воды по смене агрегатного состояния. Оно составило 2,9 грамм.
Активность, наличие / отсутствие конденсационных процессов зависит от комплекса факторов. Перепады температур и давлений определяют значительные дестабилизационные условия. При предметном и глубоком рассмотрении таких процессов возникает ряд безответных на настоящее время вопросов. Резонно было бы задаться вопросом о том, что происходит в дальнейшем с водой жидкого состояния в условиях камеры сгорания? Неизвестно, какая часть воды направляется в систему выпуска отработанных газов, а какая — в картерное пространство через сопряжения деталей цилиндро-поршневой группы.
Значительные дестабилизационные условия могут привести к частично- или полнообратному процессу смены состояний. Возможны многократные смены состояний одной группы молекул в зависимости от режимов и мгновенных положений кривошипно-шатунного механизма, а также места положения этой группы в объёме цилиндра, в рамках рабочего цикла. В настоящее время существуют работы, описывающие конденсационные процессы в поршневых двигателях [2, 4]. Известна математическая модель конденсационного процесса в картерном пространстве поршневого двигателя, прогреваемого при отрицательных температурах холодного климата [7].
ЗАКЛЮЧЕНИЕ
Математическая модель конденсационного процесса в цилиндре разработана впервые, представляет новые научные знания, позволяет решать задачи по определению массы воды, конденсирующейся на внутренней поверхности цилиндра за период прогрева.
Математическая модель применима при проектировании двигателей специфично ориентированных на условия отрицательных температур. Техника северного исполнения должна быть приспособлена к применению в таких условиях. Решение проблемы видится в применении средств тепловой подготовки, исключающих такие явления. По результатам расчётов можно судить об активности процессов, инициированных отрицательными температурами, определять потребность в средствах тепловой подготовки двигателей, определять потребную мощность этих средств, определять приемлемые теплоносители и энергоносители для их функционирования.
Дополнительно
Вклад авторов. А.В. Колунин ― поиск публикаций по теме статьи, написание текста рукописи; Е.С. Лазарев ― экспертная оценка, утверждение финальной версии; В.Н. Каминский ― редактирование текста рукописи, создание изображений; М.С. Корытов ― редактирование текста рукописи; А.О. Рузимов ― создание изображений. Авторы подтверждают соответствие своего авторства международным критериям ICMJE (все авторы внесли существенный вклад в разработку концепции, проведение исследования и подготовку статьи, прочли и одобрили финальную версию перед публикацией).
Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи.
Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования.
Additional information
Authors’ contribution. A.V. Kolunin ― search for publications, writing the text of the manuscript; E.S. Lazarev ― expert opinion, approval of the final version; V.N. Kaminsky ― editing the text of the manuscript, creating images; M.S. Korytov ― editing the text of the manuscript; A.O. Ruzimov ― creating images. Authors confirm the compliance of their authorship with the ICMJE international criteria. All authors made a substantial contribution to the conception of the work, acquisition, analysis, interpretation of data for the work, drafting and revising the work, final approval of the version to be published and agree to be accountable for all aspects of the work.
Competing interests. The authors declare that they have no competing interests.
Funding source. This study was not supported by any external sources of funding.
About the authors
Alexander V. Kolunin
Moscow Polytechnic University
Author for correspondence.
Email: kolunin2003@mail.ru
ORCID iD: 0000-0001-7151-8489
SPIN-code: 7483-9619
Associate Professor of the Power Plants for Transport and Small Energy Department
Russian Federation, 38 Bolshaya Semenovskaya street , 107023 MoscowEvgeniy S. Lazarev
Omsk State Transport University
Email: Incoe@yandex.ru
ORCID iD: 0000-0002-0226-3678
SPIN-code: 7069-0551
Senior Lecturer of the Heat Power and Thermal Engineering Department
Russian Federation, OmskValery N. Kaminskiy
Moscow Polytechnic University
Email: kamr@mail.ru
ORCID iD: 0000-0002-5430-4304
SPIN-code: 8509-5210
Dr. Sci. (Tech.), Professor of the Power Plants for Transport and Small Energy Department
Russian Federation, 38 Bolshaya Semenovskaya street , 107023 MoscowMichael S. Korytov
Siberian State Automobile and Highway University
Email: kms142@mail.ru
ORCID iD: 0000-0002-5104-7568
SPIN-code: 2921-4760
Associate Professor, Dr. Sci. (Tech.), Professor of the Road Transportation Department
Russian Federation, OmskAndrey O. Ruzimov
Ministry of Defense of the Russian Federation
Email: ruzim2009@mail.ru
ORCID iD: 0000-0002-2960-767X
Employee
Russian Federation, MoscowReferences
- Kolunin AV, Gelver SA, Buryan IA. The influence of the Arctic climate on the watering of motor oils in the operating conditions of motor transport. J. Phys.: Conf. Ser. 2019;1260:062012. doi: 10.1088/1742-6596/1260/6/062012
- Kolunin AV, Dudkin VM, Korneev SV. Water pollution and colloidal stability of motor oils. Chem Technol Fuels Oils. 2006;42:273–275. doi: 10.1007/s10553-006-0071-6
- Kolunin A.V., Kaminsky V.N., Kostyukov A.V., et al. Mathematical model of the condensation process in the crankcase space of a piston engine heated at negative temperatures of a cold climate. J. Phys.: Conf. Ser. 2021;1260:112011. doi: 10.1088/1742-6596/1260/11/112011
- Schlichting G. Boundary layer theory. Moscow: Nauka; 1974. (In Russ).
- Isachenko VP. Heat exchange during condensation. Moscow: Energiya; 1977.(In Russ).
- Kolchin AI, Demidov VP. Calculation of automobile and tractor engines. Moscow: Vysshaya shkola; 1980.(In Russ).
- Thermal calculation of boilers (normative method) 3rd ed. Saint Petersburg: NPO CCTI; 1998.(In Russ).
- Grigoriev VA., Zorin VM. Theoretical foundations of heat engineering. Thermal engineering experiment. Guide. 2nd ed. Moscow: Energoatomizdat; 1988. (In Russ).
- AVOK Reference manual. Moscow: LLC IIP AVOK-PRESS; 2014. (In Russ).