CC BY
27
DOI: 10.18721/JPM.11406 УДК 51-73: 621.574
МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКИХ ПРОцЕССОВ В ПАРОКОМПРЕССИОННОй СИСТЕМЕ ОХЛАЖДЕНИЯ
Д.Л. Карелин1, А.В. Болдырев2, В.М. Гуреев3, С.В. Болдырев2
1 Казанский (Приволжский) федеральный университет, г. Казань, Российская Федерация;
2 Набережночелнинский институт (филиал) Казанского (Приволжского) федерального университета», г. Набережные Челны, Российская Федерация;
3 Казанский национальный исследовательский технический университет им. А.Н. Туполева - КАИ, г. Казань, Российская Федерация
В работе представлена динамическая модель парокомпрессионной системы охлаждения. Особенностями модели являются учет массы рабочего агента в теплообменниках (испарителе и конденсаторе), изменения во времени паро-содержания этого агента на выходе из расширительного клапана и учет всего спектра режимов двухфазных течений при испарении рабочего агента. В ходе численного моделирования установлено, что для стабилизации температур и массового расхода в парокомпрессионной системе охлаждения требуется большее время, чем для стабилизации частоты вращения вала компрессора. Показана связь отрицательной динамики температуры испарения с начальной температурой окружающей среды, а также связь задержки стабилизации массового расхода рабочего агента с указанной начальной температурой и степенью повышения давления в термодинамическом цикле парокомпрессионной системы охлаждения.
Ключевые слова: динамическая модель, парокомпрессионная система охлаждения, теплоотдача, температурный напор, компрессор
Ссылка при цитировании: Карелин Д.Л., Болдырев А.В., Гуреев В.М., Болдырев С.В. Моделирование динамических процессов в парокомпрессионной системе охлаждения // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2018. Т. 11. № 4. С. 61-76. DOI: 10.18721/ JPM.11406
MODELING OF DYNAMIC PROCESSES IN THE VAPOR COMPRESSION COOLING SYSTEM
D.L. Karelin1, A.V. Boldyrev2, V.M. Gureev3, S.V. Boldyrev2
1 Kazan (Volga region) Federal University, Kazan, Russian Federation;
2 Kazan Federal University - Naberezhnye Chelny Institute, Naberezhnye Chelny, Russian Federation; 3 Kazan National Research Technical University named after A.N. Tupolev - KAI,
Kazan, Russian Federation
In the paper, a dynamic model of a vapor compression cooling system is presented. In addition to the usual one, it takes into account the working agent's masses in the heat exchangers, this agent's vapor content behavior in time at the outlet of the expansion valve, and the whole spectrum of two-phase flow modes during the working agent's evaporation. It was established that it took more time for temperature's and mass flow's (in a vapor compression cooling system) transitions to steady states than
for the rotational speed of the compressor shaft. The connection between the negative dynamics of the evaporation temperature and the initial ambient temperature was shown. Moreover, it was the connection between the delay in stabilization of the mass flow of the working medium and the initial ambient temperature as well as the degree of a pressure increase in the thermodynamic cycle of the vapor compression cooling system.
Keywords: dynamic model, vapor-compression cooling system, heat transfer, temperature difference, compressor
Citation: D.L. Karelin, A.V. Boldyrev, V.M. Gureev, S.V. Boldyrev. Modeling of dynamic processes in the vapor compression cooling system, St. Petersburg Polytechnical State University Journal. Physics and Mathematics. 11 (4) (2018) 61-76. DOI: 10.18721/JPM.11406
Введение
Современная тенденция развития энергетических машин неразрывно связана с проблемой постоянного обеспечения теплового состояния. Эта проблема усложняется ростом плотности тепловой энергии, подлежащей рассеянию системой охлаждения по всему диапазону изменения как нагрузки на энергетическую машину, так и внешних параметров окружающей среды. На данный момент известны различные решения, позволяющие существенно улучшить эффективность охлаждения. Одним из них является подъем температуры теплоносителя выше температуры кипения при нормальных условиях [1, 2] (высокотемпературное охлаждение), но оно имеет существенные недостатки: во-первых, не решает проблемы эффективности охлаждения при температуре окружающей среды, существенно превышающей допустимые значения, во-вторых, применимо не для всех систем [3]. В работе [4] предложен достаточно действенный способ повышения эффективности систем охлаждения и уменьшения габаритных размеров ввиду улучшения теплопередающей способности радиаторов (достигается изменением функциональной схемы) и способа передачи тепла с использованием процессов фазовых переходов теплоносителя. В своих исследованиях автор использовал вместо классических теплоносителей (вода либо водный раствор этилен-гликоля) азеотропную смесь воды (81,8 %) и анилина (18,2 %) с температурой кипения 75 °С. В этой работе автор также отметил, что решение вопроса о высоких давлениях конденсации спиртов, хладонов и фреонов может дать возможность их использования
в качестве теплоносителей.
Эффективность управления тепловой нагрузкой корабельных электрических систем заключает в себе целый ряд проблем, связанных с необходимостью рассеяния большого количества тепла, которое выделяется в ограниченном объеме этими современными энергосистемами, обладающими высокой удельной мощностью [5, 6]. Авторы отмечают, что для стабильной работы таких систем (например, электромагнитная пушка, радар, электромотор и моторный привод) требуется отдельная парокомпрес-сионная система охлаждения. Такие системы нуждаются в особых подходах при проектировании, так как, в отличие от традиционных парокомпрессионных систем, работающих по стационарному циклу сжатия, в них возможны резкие динамические изменения тепловых потоков охлаждаемых объектов. В связи с вышеизложенным, задача проектирования парокомпрессионных систем охлаждения, которые используются для управления их тепловым состоянием, должна решаться с помощью динамических моделей [7].
В статье [8] представлена статическая модель парокомпрессионной системы охлаждения на основе теоретических и эмпирических зависимостей. Массовый расход хладагента, холодопроизводительность испарителя, потребляемая мощность компрессора и коэффициент трансформации цикла определены по значениям температур хладагента в испарителе и конденсаторе, а также степени его перегрева на выходе испарителя.
В то же время в работе [9] статическая математическая модель парокомпрессион-
ной системы охлаждения получена через совместное решение уравнений баланса массы в компрессоре и теплообменниках, а также уравнений теплового баланса в конденсаторе и испарителе. Рассчитаны зависимости количества отводимого тепла от температуры и массового расхода воздуха. Проведено сравнение параметров классической и па-рокомпрессионной систем охлаждения при одинаковых габаритах теплообменников и условиях окружающей среды.
Авторами статьи [10] предложен новый метод расчета температуры конденсации, учитывающий тип выбранного хладагента, изменение температуры и массового расхода охлаждающего воздуха, а также тепловую нагрузку на парокомпрессионную систему охлаждения. На основе указанного метода, в статье [11] разработана методика расчета параметров термодинамического цикла данной системы.
В работе [12] предложена математическая модель и проведено моделирование па-рокомпрессионной системы охлаждения с приводом от газового двигателя. В расчетах авторы использовали: уравнения полуэмпирической модели спирального компрессора; зависимость потребляемой мощности газового двигателя от частоты вращения его вала, температуры окружающего воздуха и мощности компрессора; уравнения модели пластинчатого теплообменника-испарителя.
Следует отметить, что в рассмотренных работах предметом интереса были параметры работы парокомпрессионных систем охлаждения в установившемся режиме. Однако в действительности указанные параметры изменяются во времени как при запуске и отключении, так и при изменении условий окружающей среды и прочих условий. Следовательно, для оптимального управления работой компрессора, для обеспечения требуемых быстродействия и запаса устойчивости парокомпрессионных установок и других параметров актуальным является исследование их динамических характеристик.
Так, в статье [13] представлена математическая модель динамики парокомпрес-сионной системы охлаждения на основе
законов сохранения массы и энергии в испарителе, компрессоре, конденсаторе и расширительном клапане. Исследованы эксплуатационные характеристики системы — холодопроизводительность и холодильный коэффициент цикла. Определены такие параметры, как установившийся массовый расход хладагента, степень его перегрева в испарителе и другие. При анализе временных зависимостей частоты вращения вала компрессора и массового расхода рабочего агента отмечено высокое быстродействие системы.
В работе [14] выполнено экспериментальное и численное исследование динамики теплового насоса, работающего в режиме охлаждения и снабженного приводом спирального компрессора от газового двигателя. Математическая модель представлена уравнениями баланса энергии, критериальными уравнениями теплоотдачи для двухфазных течений в теплообменниках и уравнениями для расчета параметров газового двигателя. Давления всасывания и нагнетания приняты равными давлениям, при которых происходит испарение и конденсация хладагента соответственно. Дифференциальные уравнения решены методом Рунге — Кутта. Исследовано изменение во времени температур испарения и конденсации, холодопроизводительности, расхода топлива газовым двигателем и потребляемой мощности на валу компрессора.
Тем не менее, в представленных работах не рассматривалось изменение во времени паросодержания в испарителе и масс хладагента в теплообменниках, что оказывает значительное влияние на точность результатов моделирования, поскольку указанные характеристики могут существенно изменяться в течение переходных процессов, например во временном интервале между запуском установки и моментом достижения ею стационарного режима работы. В то же время, учет данных факторов позволит повысить точность прогнозирования как времени нагрева и охлаждения энергетических агрегатов, теплообменников и аппаратов, так и других характеристик переходных процессов, информация о которых необходима при проектировании блоков
управления парокомпрессионных систем охлаждения.
В настоящей работе разработана математическая модель динамики парокомпресси-онной системы охлаждения (рис. 1). Модель основана на решении уравнения ротационного пластинчатого компрессора и дифференциальных уравнений теплообменника-конденсатора, теплообменника-испарителя. Кроме того, модель учитывает массы рабочего агента, содержащиеся в теплообменнике-конденсаторе и теплообменнике-испарителе, изменение во времени паросодержания, а также теплоотдачу при кипении рабочего агента в соответствии с диаграммой режимов двухфазных течений. С помощью модели рассчитаны переходные характеристики, в частности изменение во времени температур конденсации и испарения (кипения), а также массового расхода рабочего агента, нагнетаемого компрессором.
Устройство и принцип действия парокомпрессионной системы охлаждения
Моделируемая система работает по термодинамическому циклу, подобному тепловому насосу и кондиционеру. Она включает в себя парокомпрессионный кон-
тур 3^4 (рис. 1, b) циркуляции ра-
бочего агента (фреон-132Ь), содержащий одноступенчатый ротационный пластинчатый компрессор (Com), теплообменник-конденсатор (CHE), расширительный клапан (регулируемый дроссель Thr) и теплообменник-испаритель (EHE). Теплота подводится к теплообменнику-испарителю из контура охлаждения, включающего насос (P), обеспечивающий циркуляцию теплоносителя (тосол-65), и рубашку охлаждения (CJ). Для отвода теплоты теплообменник-конденсатор обдувается воздухом из вентилятора (Air).
Динамическая модель парокомпрессионной системы охлаждения
Принятые допущения. При разработке математической модели динамики паро-компрессионной системы охлаждения приняты следующие допущения:
в дифференциальных уравнениях теплового баланса теплообменников при расчете коэффициента теплопередачи не учитывались термические сопротивления загрязнений и паровой пленки;
температура на границе внутренней стенки трубок теплообменников и паровой пленки считалась равной температуре про-
Рис. 1. Принципиальная схема парокомпрессионной системы охлаждения (a) и ее термодинамический цикл (b):
Com — одноступенчатый ротационный пластинчатый компрессор, CHE — теплообменник-конденсатор, Thr — расширительный клапан, EHE — теплообменник-испаритель, P — насос, CJ — рубашка охлаждения; Air — воздух от вентилятора; T , P, Tv, Pv — температуры и давления процессов конденсации и испарения, соответственно; Qv, Qc — тепловые мощности испарителя и конденсатора; P(i) — термодинамическая диаграмма P-i рабочего агента (зависимость давления от удельной энтальпии для жидкой и паровой фаз); жирным пунктиром показаны изотермы
цесса конденсации Tc и испарения Tev, соответственно (см. рис. 1);
масса рабочего агента рассчитывалась по гомогенной модели, предполагающей равномерное распределение паровой и жидкой фаз и не учитывающей эффект скольжения фаз на границе раздела;
процесс дросселирования рабочего агента в расширительном клапане считался изоэнтальпийным;
температуры и массовые расходы охлаждаемой (тосол-65) и охлаждающей (воздух) сред считались постоянными;
не учитывались тепловые потери в окружающую среду.
Компрессор. Удельная работа ¡а адиабатного сжатия рабочего агента в компрессоре определялась по уравнению:
_ ^ пт
k-Т
грхт
- Т
V" т у
(Т)
где Р1, Р2, Па, — давления в точках 1 и 2термодинамического цикла; у1, м3/кг, — удельный объем; к — показатель политропы.
Потребляемая (индикаторная) мощность N. на валу компрессора:
N _ ^aGcom
IV j —
nad
(2)
где Осот, кг/с, — массовый расход рабочего агента (фреон-132Ь); цаа — адиабатный кпд, который для ротационного пластинчатого компрессора принимался равным 0,82 [15].
Объемный расход пара Уоот (м3/с), всасываемого ротационным пластинчатым компрессором, рассчитывался (без учета толщины пластин) по следующей зависимости:
V _ 2nD I s n X ,
com com com com com com'
(3)
где Х>оот, м, — диаметр рабочей поверхности корпуса (цилиндра); ¡оот, м, — длина ротора (поршня); воот, м, — смещение центров образующих цилиндра и ротора (эксцентриситет); пот, с-1, — число оборотов ротора в минуту; Хоот — коэффициент расхода, зависящий от степени повышения давления Пот = Р2 / Р и размеров компрессора.
Из расчета допустимой нагрузки на пластину ротационных компрессоров, в зависимости от степени повышения давления %оот, приняты соотношения между основными геометрическими размерами его рабочих элементов [15, 16]:
I 2s
_ 1)8; _ 0,140
D
ПРИ ncom _ Р/ Р ^ 2>5;
(4)
L
_ 1,8;
2s„
_ 0,ТТ5
(5)
D D
com com
при пcom _ PJ Р Z 5,0.
В данной работе степень повышения давления принята равной ncom _ 4,52, поэтому с учетом выражений (5) и соотношений для радиуса поршня компрессора Rcom _ Dcoml2 и массового расхода Gcom _ Vcom/vT уравнение массового расхода ротационного пластинчатого компрессора принимает вид:
(/) _ 1,664nCAom
■Л''
(6)
Эмпирическая зависимость Хсот от степени повышения давления поот в уравнении (6) принята следующей [15]:
1 _ 1,0 - a
г Р ^
V рт у
(7)
где а = 0,05 и 0,10 для крупных и для мелких компрессоров, соответственно.
Уравнение динамики температуры испарения. Конструкция теплообменника-испарителя принята кожухотрубной и многоходовой. Рабочий агент (фреон-132Ь) движется внутри пучка труб, который снаружи обтекается тосолом. Баланс тепловой мощности теплообменника-испарителя имеет следующий вид:
0 = в* - О , (8)
^оот.еу'
где , Вт, — изменение энтальпии рабочего агента в теплообменнике-испарителе за единицу времени; 0л>еу, Вт, — мощность теплопередачи в теплообменнике-испарителе; 0 , , Вт, — тепловая мощность, отводимая компрессором от теплообменника-испарителя.
Изменение энтальпии рабочего агента
v
в теплообменнике-испарителе за единицу времени выражается как
dT
Cp' + Cp" Q = -¿lev.. m
dt
(9)
где Ср' ет, Ср' ' ет, Дж/(кгК), - изобарные теплоемкости жидкой и газообразной фаз рабочего агента, соответственно; ткг, — масса рабочего агента в теплообменнике-испарителе; Теу, К, — температура испарения рабочего агента в теплообменнике-испарителе; ^ с, — время.
Тепловая мощность О , отводи-
1 ^сот^
мая компрессором от теплообменника-испарителя выражается как
О = (£ - ¡4)в , (10)
^сот.еу V 1 4' сот'
где ¡1, ¡4, Дж/кг, — удельные энтальпии рабочего агента в точках цикла (см. рис. 1, Ь).
Мощность теплопередачи в
теплообменнике-испарителе, с учетом среднелогарифмического температурного напора, примет следующий вид:
О = Ж 1
^й.еу еу
1
a„
+ ^ + J_
Ь„ a„,
cf .ev
(Tf - Tev ) - (Tf- Tev)
(11)
(
ln
_V
T' - T
1cf 1 ev
'T'tt _ у
Tcf - Tev
\
где F
^ ev
щей
м2, — площадь теплопередаю-поверхности теплообменника-
испарителя; 5p ev, м, — толщина стенки трубки теплообменника-испарителя; Tcf , Tf K, — температуры теплоносителя (тосол-65, англ. cooling fluid) на входе и выходе теплообменника-испарителя соответственно; aev, Вт/(м2-К), — коэффициент теплоотдачи при кипении рабочего агента в круглых трубах; acfev, Вт/(м2-К), — коэффициент теплоотдачи тосола для кожухотрубного многоходового теплообменника-испарителя; bm, Вт/(м-К), — коэффициент теплопроводности материала стенки теплообменника-испарителя.
На основании баланса тепловой мощности (8) и уравнений подводимой и отводимой тепловых мощностей (9) — (11), нами получено следующее уравнение динамики температуры испарения:
mev (t).
= F„.
(Cp'ev + CpV ) dTv
dt
1
a
cf .ev m
Ьm av (t)
(12)
(Tf - Tev) - (Tf - Tev) - (• - iGm (t).
ln
( T' - T ^
Jcf *-ev
rptt _ ГГ1
V Tcf ~ Tev J
Коэффициент теплоотдачи при кипении рабочего агента в круглых трубах следует выражению [22]:
aev (t) = 0,28 . 10-
Ь 'ev (t)
(u'' (t). d . Л0'19 (u' (t). d . ^0'66 (13)
ev\ ' ev in ev\S ev in
V ''ev(t)
V'ev (t)
V 'ev (t) v « 'ev (t) J
ev
( Г V
d .
V ev _ in j
которое справедливо для разделенных режимов (расслоенного с гладкой и волновой границей двухфазного течения);
aev(t) = 0,65 . 10
U ''ev(t) . dv i,
5 Ь 'ev (t)
(u '' (t) . d Л0,73 (u ' (t) . d V0'73
V ''ev(t)
d .
ev _ in
u' (t). d -
ev\ ' ev in
V'ev (t)
(14)
У0,3 ( т V
1,69
V 'еу (0 V в 'еу И) ,
для перемежающихся режимов (снарядного, снарядно-кольцевого двухфазного течения);
aev (t) = 0,018 .
Ь 'ev (t) ( U ''ev (t) . deV_ Л
1,19
(u' (t). d . V0'3 (v ^V0,3
^ ev\ ' ev in
V 'ev (t)
V ''ev (t)
V 'ev (t) « 'ev (t)
(15)
V еу\ / у
для дисперсных режимов (дисперсного, дисперсно-кольцевого и кольцевого двухфазного течения),
где X 'еу, Вт/(м-К), — коэффициент теплопроводности жидкой фазы рабочего агента в испарителе; Пп, м, — внутренний диаметр трубок в испарителе; Ьеу, м, — длина
1
X
трубок в испарителе; и , и, м/с, — приведенные функции скорости движения газообразной и жидкой фаз рабочего агента в испарителе, соответственно; V , V , м2/с, — кинематические коэффициенты вязкости газообразной и жидкой фаз рабочего агента в испарителе; а' еу, м2/с, — коэффициент температуропроводности жидкой фазы рабочего агента в испарителе.
В уравнениях (13) — (15) и и и — приведенные функции скорости движения газообразной и жидкой фаз рабочего агента, рассчитываемые по следующим соотношениям:
() _
G (t) • х (t)
com\ / ev У / . Plv (t) • Sp.ev • Zev
' (t) _ Gom (t) •(Т - Xev (t)) Uev (t)
(Т6)
(Т7)
Р'еу (*) ' Зр.еу ' ¿ег
где — количество параллельных трубок в теплообменнике-испарителе, х^ — массовое паросодержание в рабочей точке 4 (см. рис. 1, Ь);
^) - '4(*)
Xev (t) _
(Т8)
1 (*) - »4(*)
где ¡1'', /4', Дж/кг, — удельные энтальпии рабочего агента в точках цикла (см. рис. 1, Ь).
Коэффициент теплоотдачи тосо-ла для кожухотрубного многоходового теплообменника-испарителя определялся по критериальному уравнению для обтекания потоком жидкости пучка труб с шахматным расположением [15]:
а
cf .ev
_ 0,41
Ucf .ev dev.out
cf .ev
0,6
Л
(19)
cf .ev
V acf .ev у
где X 'о/еу, Вт/(м-К), — коэффициент теплопроводности тосола в испарителе; ¿еуоШ, м, — наружный диаметр трубок в испарителе; , м/с, — максимальная скорость тосола
cf .ev
в пучке трубок в испарителе в самом узком сечении; vcfv, м2/с, — кинематический коэффициент вязкости тосола в испарителе; м2/с, — коэффициент температуро-
a
cf .ev'
проводности тосола в испарителе.
В этом уравнении за определяющий размер принят наружный диаметр (^ел.оШ, м) обтекаемой трубы пучка, а за определяющую скорость (, м/с) — максимальная скорость в пучке в самом узком сечении.
Уравнение для расширительного клапана. После полной конденсации рабочий агент проходит через расширительный клапан, где изоэнтальпийно расширяется (процесс 3 ^ 4, см. рис. 1, Ь) с понижением давления от давления Рс конденсации до давления испарения с переходом части его массы в паровую фазу:
4 =/4, (20)
где /3, /4, Дж/кг, — удельные энтальпии рабочего агента в точках цикла (см. рис. 1, Ь).
Уравнение динамики температуры конденсации. В работе принято классическое конструктивное исполнение конденсатора как теплообменного аппарата, состоящего из медных труб с оребрением, внутри которых течет рабочий агент (фреон-132Ь), а снаружи движется воздух. Тогда баланс его тепловой мощности имеет вид:
Q _Q -Qh ,
(2Т)
где Qc, Вт, — изменение энтальпии рабочего агента в теплообменнике-конденсаторе за единицу времени; Qcomc, Вт, — тепловая мощность, подводимая компрессором в теплообменник-конденсатор; Qthc, Вт, — мощность теплопередачи в теплообменнике-конденсаторе.
Изменение энтальпии рабочего агента в теплообменнике-конденсаторе за единицу времени выражается как
Q _ CP' + CP'-т • T, (22) 2 c dt'
где Cp'c, Cp"c, Дж/(кгК), — изобарные теплоемкости жидкой и газообразной фазы рабочего агента в конденсаторе, соответственно; m, кг, — масса рабочего агента в теплообменнике-конденсаторе; T, К, — температура конденсации рабочего агента в теплообменнике-конденсаторе.
Тепловая мощность, подводимая компрессором в теплообменник-конденсатор, определяется как
Q _ (i2 - i3)G , (23)
s^com.c У 2 3' com'
ev.out
где ¡2, ¡3, Дж/кг, — удельные энтальпии рабочего агента в точках цикла (см. рис. 1, Ь).
Мощность теплопередачи в
теплообменнике-конденсаторе с учетом среднелогарифмического температурного напора имеет следующий вид:
О = Ж •- 1
«¿й.с 1 с
1
а„
X.
1
а„
(Тс - т:) - (Тс - Та")
(24)
1п
Т - Т'
с_а_
гр грп с а J
где ¥с, щей
м2, — площадь теплопередаю-поверхности теплообменника-конденсатора; 8(с, м, — толщина стенки трубки теплообменника-конденсатора; Хт, Вт/(м^К), — коэффициент теплопроводности материала трубки теплообменника-конденсатора; Т, Т', К, — температуры воздуха на входе и выходе теплообменника-конденсатора соответственно; ас, Вт/(м2-К), — коэффициент теплоотдачи при конденсации рабочего агента в круглых трубах; аа, Вт/(м2-К), — коэффициент теплоотдачи для воздуха в конденсаторе-теплообменнике.
На основе баланса тепловой мощности (21) и уравнений подводимой и отводимой тепловых мощностей (22) — (24) получено следующее уравнение динамики температуры конденсации:
тс а)
Ср'с + Ср"с ¿Тс
- ¥,
2 1
= (¡2 - 3)Оо0т а) -
X с/, Вт/(м-К), — коэффициент теплопроводности фреона132Ь; йсЫ, м, — внутренний диаметр трубки теплообменника-конденсатора; Ке, Рг — числа Рейнольдса и Прандтля для потока хладагента внутри трубы, соответственно.
Критерий Рейнольдса для потока хладагента внутри трубы в уравнении (26) вычислялся по формуле:
Ке = . (27)
Среднее значение коэффициента теплоотдачи при конденсации рабочего агента в круглых трубах ас для двухфазной области 2"— 3' (см. рис. 1, Ь), в соответствии с [20], определялось из уравнения:
ас(0 = аг(0 •
5 2,04
9 + Рг0-38
«)
(28)
где а,, Вт/(м2-К), — коэффициент теплоотдачи, используемый для расчета в однофазной области по уравнению (26).
Коэффициент теплоотдачи аа для воздуха в конденсаторе-теплообменнике, представляющем собой пучок труб шахматного расположения с наружным оребрением круглой формы [21]:
а: = 0,23-
X.
и • Ь
\ 0,65
(и ?'54 (и ч0Д4
(29)
V ^с.оШ J
/п
5 5 ,
с
где 5с= (^1с- ^/Я2С - )0'2 — коэф-
/гр Т>\ /Т гГг/\ с 4 с е.ош/ с е.ош'
(Тс - Та) - (Тс - Та). (25) фициент; 5z — коэффициент, учитываю
а„
ас«)
1п
Т - Т'
с а
Т Т"
Для определения коэффициента теплоотдачи а,, Вт/(м2-К), в однофазной области перегрева 2 — 2 ' ' и переохлаждения 3' — 3 (см. рис. 1, Ь) рабочего агента использовалось критериальное уравнение [15, 18]:
а / 0) =
Хс, Ц)
• 0,023 • Ке0'8^) • Ргп (*), (26)
где п — показатель степени для процесса конденсации, п = 0,4 [18, 19];
щий количество вертикальных рядов в теплообменнике-конденсаторе; 51с, £2с, м, — продольный и поперечный шаг трубок, соответственно; ¿соШ, м, — наружный диаметр трубки теплообменника-конденсатора; Ь, м, — шаг оребрения; кГт, м, — высота ребра; Ха с, Вт/(м-К), — коэффициент теплопроводности воздуха в конденсаторе; иа с, м/с, — скорость воздуха в пучке трубок в конденсаторе; \а с, м2/с, — кинематический коэффициент вязкости воздуха в конденсаторе.
Уравнение (29) справедливо в диапазоне чисел Рейнольдса Ке = 300 — 22500. За
Ь
определяющий размер принят шаг оребре-ния Ь, а за определяющую скорость — максимальная скорость в узком сечении:
(
= О,.
рА,
(
к.
1 + 2•-Ь й
1 _ йс. ом 81
с
^Т1
(30)
/п
с.ош / у у
где ра, кг/м3, — плотность воздуха; 8 п, м, — высота и толщина ребра; ¥с, м2, — фронтальная площадь теплообменника-конденсатора; (?а, кг/с, — массовый расход воздуха через конденсатор.
Масса рабочего агента в теплообменнике-испарителе и конденсаторе. Теоретическая масса рабочего агента т (кг) в теплообменнике-испарителе и в теплообменнике-конденсаторе, при условии замкнутой системы, принималась постоянной и определялась для номинального режима работы по формуле
V
| й I
т = т/ + т2 =
14 й1 + 1 Р/^Т й1
(31)
где т/Г, т2, кг,
массы жидкой и газообразной фаз рабочего агента, соответственно; У, м3, — внутренний объем трубок; I, м, — продольная координата; р,, ру, кг/м3, — плотности жидкой и газообразной фаз рабочего агента, соответственно; Б, м2, — площадь проходного сечения трубок; Б,, , м2, — площади жидкой и газообразной фаз рабочего агента, соответственно; Ь, м, — эквивалентная длина трубок (У/Б).
С учетом того, что удельное объемное паросодержание а' для элементарного объема равно
а' =-2—, (32)
2 /
формулу (31) можно записать в следующем виде:
V
т = — Ь
| (Р„ а' + (1 _а') р г) й1
(33)
Удельное массовое паросодержание X
определяется отношением массового расхо да паровой фазы пгу (кг/с) к полному мас совому расходу хладагента ^ + т/ (кг/с):
т„
X =
т„ + т
(34)
7
Если учесть, что т = р -5" • и, то удельное объемное паросодержание (32) можно привести к уравнению вида
а =
1 +
1 _ X Лр
X
Р 7
■у
(35)
где у = и / — коэффициент скольжения фаз на границе раздела (у = 1 для гомогенной модели).
После совместного решения уравнений (33) и (35) получим выражение для расчета массы рабочего агента тп (кг), например, в теплообменнике-испарителе, при номинальном режиме работы с учетом паросо-держания в рабочей точке 4 термодинамического цикла (см. рис. 1, Ь):
Ь„, г _ Р Л
V,.
т=
Р 7 .е
1 + 71
й1, (36)
где
71 = [((1 _ Хет) • ¡¡Ьеу + хет)-1 _ 1](Р_/Р/.„)
— функция, учитывающая изменение удельного массового паросодержания X по длине Ьеу теплообменника-испарителя (при I = 0 паросодержание равно его текущему значению в рабочей точке 4, т. е. X = хеу (см. рис. 1, Ь); при I = Ь^ паросодержание принимает значение X = 1); р^, р/е2, кг/м3, — плотности газообразной и жидкой фаз рабочего агента в испарителе; Уеу, м3, — внутренний объем трубок испарителя.
Соответственно выражение для расчета массы тс (кг) рабочего агента в теплообменнике-конденсаторе, с учетом предположения, что процесс конденсации происходит полностью и удельное массовое паросодержание меняется в диапазоне от 1 до 0, имеет вид
т = Ь1
Ьс 0
Pv.cc
(Р7.С _Рг.с )
Л
1 + 7
й1, (37)
72 = [(1 _ //4 )-1 _ 1] • (РV.с/Р7.с)
и с
Рис. 2. Изменение удельного массового паросодержания X: а — по длине теплообменника-испарителя Ьеу при кипении хладагента; Ь — по длине
теплообменника-конденсатора Ьс при конденсации хладагента; х — массовое паросодержание в рабочей точке 4 цикла, I — продольная координата
— функция, учитывающая изменение паросодержания X по длине Ьс теплообменника-конденсатора (при I = 0 паросодержание равно X = 1, при I = Ьс паросодержание принимает значение X = 0); Ьс, м — длина трубок в конденсаторе; V, м3 — внутренний объем трубок конденсатора; рус, р/с, кг/м3, — плотности газообразной и жидкой фаз рабочего агента в конденсаторе.
На рис. 2 изменение удельного массового паросодержания Х представлено на графиках.
Моделирование, его результаты и их обсуждение
Уравнения динамики (12) и (25) температур в испарителе и конденсаторе (соответственно) решались численным методом Рунге — Кутты 4-го порядка с фиксированным шагом. Коэффициенты теплоотдачи при испарении (аеу, см. формулы (13) — (15)) и конденсации (ас, (28)) в двухфазной области, а также массы рабочего агента в испарителе (теу, (36)) и конденсаторе (тс, (37)) являлись функциями массового расхода (£сот(0, (6)), который, в свою очередь, зависел от частоты вращения вала компрессора. При этом процесс разгона вала компрессора при включении установки описывался экспоненциальной функцией (п, об/мин, — установившееся значение частоты вращения вала; ? , мин, — время разгона вала):
(
п(?) = п 1 - ехр
(
-7?
JJ
(38)
Расчеты выполнены при допущении, что в начальный момент времени ? = 0 массы хладагента в конденсаторе тс.0 (кг) и испарителе теу0 (кг) распределены пропорционально внутренним объемам этих теплообменников V и Уу, при давлении насыщенного пара, соответствующем температуре окружающей среды, равной 25°С для первого случая и 50°С для второго:
mv
V + V
еу с
V,. + V
V;
V
г еу'
(39)
(40)
где т2 = тс + , кг, — суммарная масса рабочего агента в теплообменниках (испарителе и конденсаторе) в номинальном режиме работы парокомпрессионной системы охлаждения.
Значения теплоемкостей жидкой и газообразной фаз (Ср' и Ср'') рабочего агента в конденсаторе и испарителе, удельных энтальпий ¡ в рабочих точках цикла (см. рис. 1, Ь), коэффициентов теплопроводности X, динамической вязкости р, и плотностей жидкой и газообразной фаз (р 'и р'') определялись в каждый момент времени как функции температуры испарения и конденсации, путем интерполяции данных в справочных таблицах [23].
По результатам численного моделирования были построены переходные характеристики системы в моменты запуска из режима ожидания при температурах окружающей среды 25 и 50°С, которые пред-
Рис. 3. Переходные характеристики температур испарения Ту (а, Ь) и конденсации Тс (с, Л) при температурах окружающей среды плюс 25°С (а, с) и 50°С (Ь, с(); установившиеся значения Т и Т показаны пунктиром
Рис. 4. Переходные характеристики массового расхода рабочего агента при температурах окружающей среды плюс 25°С (а) и 50°С (Ь); установившиеся значения массового расхода показаны пунктиром
О 1 2 3 4
Рис. 5. Динамика частоты вращения вала компрессора; установившееся значение частоты показано пунктиром
ставлялись в виде зависимостей температур рабочего агента в испарителе и конденсаторе и соответствующих им удельных объемов рабочего агента (получены интерполяцией данных в таблицах [21]) (рис. 3 — 5).
Отмечено большее время стабилизации температур испарения, конденсации (рис. 3) и массового расхода (рис. 4) относительно установления номинальных оборотов вала компрессора (рис. 5). Так, температуры испарения и конденсации рабочего агента достигают установившегося значения в среднем за 4,5 с, причем независимо от значений температуры окружающей среды в момент запуска системы (и перед ним).
Обнаружено некоторое снижение температуры испарения рабочего агента в течении 0,1 с после пуска системы при температуре окружающей среды плюс 50°С (см. рис. 3, Ь). Наблюдающийся вслед за этим резкий спад температуры испарения вызван сменой режима двухфазного течения с разделенного на перемежающийся (первый излом характеристики). При переходе к дисперсному режиму (второй излом характеристики) зафиксирован плавный рост температуры испарения до установившегося значения. Причиной отрицательной динамики температуры испарения, согласно уравнению (12), является малое значение коэффициента теплопередачи, вызванное сравнительно небольшим логарифмическим температурным напором и невысокими приведенными скоростями и" , и'
(см. формулы (16) и (17)) движения фаз рабочего агента при значительно большей величине тепла, отведенного компрессором.
Для варианта с температурой окружающей среды 25°С, существенно увеличившийся температурный напор в теплообменнике-испарителе компенсирует малое значение коэффициента теплопередачи при малых приведенных скоростях движения рабочего агента. Поэтому температура испарения Теу (см. рис. 3, а) возрастает до того момента времени, когда изменение значений приведенных скоростей приведет к смене режимов двухфазного течения. На перемежающемся режиме происходит уменьшение температуры испарения. После установления в трубках теплообменника-испарителя дисперсного течения, температура испарения Теу плавно возрастает до установившегося значения.
В ходе анализа зависимости изменения переходных характеристик массового расхода хладагента С^оот (см. рис. 4) и числа оборотов в единицу времени поот вала компрессора от времени (см. рис. 5) было установлено, что расход стабилизируется медленнее, чем частота вращения вала, причем задержка составляет 1,2 с при температуре окружающей среды 25°С и 0,8 с - при 50°С.
Очевидно, что на время стабилизации температур испарения Т и конденсации Тс в парокомпрессионной системе охлаждения, согласно уравнениям динамики (12) и (25), оказывает влияние теплосодержание рабочего агента в испарителе и конденсаторе, определяемое его массой, а также величина теплоотдачи, зависящая от режима двухфазного течения в испарителе.
Увеличение времени стабилизации массового расхода компрессора Ооот, согласно уравнению (5), связано с тем, что удельный объем рабочего агента в рабочей точке 1 (см. рис. 1, Ь) зависит от температуры испарения Т^ и утечек Xсот, величина которых возрастает с увеличением степени повышения давления псот (6).
Заключение
В данной работе представлена динамическая модель парокомпрессионной системы охлаждения, дополнительно учитывающая
(посредством критериальных уравнений) массы рабочего агента в теплообменниках (испарителе и конденсаторе), изменение во времени паросодержания рабочего агента в точке 4 термодинамического цикла на выходе из расширительного клапана и весь спектр режимов двухфазных течений при испарении рабочего агента. В ходе численного моделирования установлено, что для стабилизации температур и массового расхода в парокомпрессионной системе охлаждения требуется большее время, чем для частоты вращения вала компрессора. Замедление зависит от значений масс рабочего агента в теплообменниках и коэффициента теплоотдачи в испарителе, существенно изменяющегося в зависимости от
режима двухфазного течения. Кроме того, задержка стабилизации массового расхода связана с утечками рабочего агента, которые увеличиваются пропорционально степени повышения давления псот, а также с влиянием начальных условий работы системы, в частности исходной температуры окружающей среды.
На переходных характеристиках выявлены временные интервалы, когда после пуска системы из состояния ожидания наблюдалось снижение температуры испарения со временем. Установлена связь отрицательной динамики температуры Теу с исходной температурой окружающей среды, а также со сменой режима двухфазного течения.
СПИСОК ЛИТЕРАТУРЫ
1. Патрахальцев Н.Н., Савастенко А.А. Форсирование двигателей внутреннего сгорания наддувом. М.: Легион—Автодата, 2004. 176 с.
2. Кривов В.Г., Синатов С.А., Ким Ф.Г., Устинов Н.А. Теплоотвод в зарубашечное пространство форсированного тепловозного дизеля при его высокотемпературном охлаждении // Двигателестроение. 1986. № 11. С. 5—11.
3. Соколовский Э.И., Улитенко А.И., Пят-кина Т.А. Оптимизация режима теплоотвода от твердотельных СВЧ-приборов в условиях повышенной температуры окружающей среды // Электроника: Межвуз. сб. науч. трудов. Рязань, 1978. С. 128-131.
4. Склифус Я.К. Фазовые переходы теплоносителя в системе охлаждения дизеля тепловоза // Труды Ростовского государственного университета путей сообщения. 2014. № 4 (29). С. 92-95.
5. Dilay Е., Vargas J.V.C., Souza J.A., Ordonez J.C., Yang S., Mariano A.B. A volume element model (VEM) for energy systems engineering // International Journal of Energy Research. 2015. Vol. 39. No. 1. Pp. 46-74.
6. Nunes T.K., Vargas J.V.C., Ordonez J.C., Shah D., Martinho L.C.S. Modeling, simulation and optimization of a vapor compression refrigeration system dynamic and steady state response //Applied Energy. 2015. Vol. 158. No. 15. Pp. 540-555.
7. Catano J., zhang T., Wenc J.T., Jensen M.K., Peles Y. Vapor compression refrigeration cycle for electronics cooling. Part I: Dynamic modeling and experimental validation // International Journal of Heat and Mass Transfer. 2013. Vol. 66. November. Pp. 911-921.
8. Cabello R., Torrella е., Navarro-Esbri J.
Experimental evaluation of the internal heat exchanger influence on a vapor-compression plant performance using R134a, R407C and R22 as working fluids // Applied Thermal Engineering. 2004. Vol. 24. No. 13. Pp. 1905-1917.
9. Карелин Д.Л., Гуреев В.М., Мулюкин
B.Л. Моделирование системы охлаждения с парожидкостной компрессионной установкой // Вестник КГТУ им. Туполева. 2015. № 5.
C. 5-10.
10. Карелин Д.Л., Гуреев В.М. Метод расчета температуры конденсации рабочего агента для парожидкостных компрессионных систем охлаждения // Вестник КГТУ им. Туполева. 2016. № 4. С. 20-24.
11. Карелин Д.Л. Методика расчета параметров термодинамического цикла парокомпрес-сионной системы охлаждения // Труды Академ-энерго. 2017. № 3. С. 23-31.
12. Elgendy Е., Schmidt J., Khalil A., Fatouh M. Modelling and validation of a gas engine heat pump working with R410A for cooling applications // Applied Energy. 2011. Vol. 88. No. 11. Pp. 4980-4988.
13. Lei zhao, zaheeruddin M. Dynamic simulation and analysis of a water chiller refrigeration system //Applied Thermal Engineering. 2005. Vol. 25. No. 14-15. Pp. 2258-2271.
14. Sanaye S., Chahartaghi M., Asgari H. Dynamic modeling of gas engine driven heat pump system in cooling mode // Energy. 2013. Vol. 55. 15 June. Pp. 195-208.
15. Розенфельд Л.М., Ткачев А.Г. Холодильные машины и аппараты. М.: ГИТЛ, 1960. 666 с.
16. Кошкин Н.Н. Тепловые и конструктивные расчеты холодильных машин. Л.: Машиностроение, 1976. 402 с.
17. Klimenko V.V. A generalized correlation for two-phase forced flow heat transfer // International Journal of Heat and Mass Transfer. 1988. Vol. 31. No. 3. Pp. 541-552.
18. 2001 ASHRAE handbook: Fundamentals. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, 2001. 544 р.
19. Jabardo J.M.S., Mammani W.G., Ianella M.R. Modeling and experimental evaluation of an automotive air conditioning system with a variable capacity compressor // International Journal of Refrigeration. 2002. Vol. 25. No. 8.
Pp. 1157-1172.
20. Shah M.M. A general correlation for heat transfer during film condensation in tubes // International Journal of Heat and Mass Transfer. 1974. Vol. 22. No. 4. Pp. 547-556.
21. Авчухов В.В., Паюсте Б.Я. Задачник по процессам тепломассообмена. М.: Энергоато-миздат, 1986. 144 с.
22. Чумак И.Г., Чепурненко В.П., Чуклин С.Г. Холодильные установки. 2-е изд., перераб. и доп. М.: Легкая и пищевая промышленность, 1981. 344 с.
23. Максимов Б.Н., Барабанов В.Г., Серуш-
кин И.Л. Промышленные фторорганические продукты. Ленинград: Химия, 1990. 464 с.
Статья поступила в редакцию 13.05.2018, принята к публикации 02.11.2018.
СВЕДЕНИЯ ОБ АВТОРАХ
КАРЕЛИН Дмитрий Леонидович — кандидат технических наук, доцент кафедры высокоэнергетических процессов и агрегатов Казанского (Приволжского) федерального университета, г. Казань, Российская Федерация.
420008, Российская Федерация, Республика Татарстан, г. Казань, Кремлевская ул.,18 [email protected]
БОЛДыРЕВ Алексей Владимирович — кандидат технических наук, доцент кафедры высокоэнергетических процессов и агрегатов Казанского (Приволжского) федерального университета, г. Казань, Российская Федерация.
423812, Российская Федерация, Республика Татарстан, г. Набережные Челны, пр. Мира, 68/19 (1/18).
ГУРЕЕВ Виктор Михайлович — доктор технических наук, заведующий кафедрой теплоэнергетики и энергетического машиностроения Казанского национального исследовательского технического университета им. А.Н. Туполева — КАИ, г. Казань, Российская Федерация.
420111, Российская Федерация, Республика Татарстан, г. Казань, ул. К. Маркса, 10. [email protected]
БОЛДыРЕВ Сергей Владимирович — кандидат технических наук, доцент кафедры высокоэнергетических процессов и агрегатов Казанского (Приволжского) федерального университета, г. Казань, Российская Федерация.
423812, Российская Федерация, Республика Татарстан, г. Набережные Челны, пр. Мира, 68/19 (1/18).
REFERENCES
[1] N.N. Patrakhaltsev, A.A. Savastenko,
Forsirovaniye dvigateley vnutrennego sgoraniya nadduvom [Boosting internal combustion engines], Legion— Avtodata, Moscow (2004).
[2] V.G. Krivov, S.A. Sinatov, F.G. Kim, N.A. ustinov, Teplootvod v zarubashechnoye prostranstvo forsirovannogo teplovoznogo dizelya pri yego vysokotemperaturnom okhlazhdenii [Heat removal to the volume of the boosted diesel engine with
its high-temperature cooling], Dvigatelestroyeniye. (11) (1986) 5-11.
[3] Eh.I. Sokolovskj, A.I. ulitenko, T.A. Pyatkina, Mode optimization of the heat removal from the solid-state microwave devices under high-temperature ambient conditions // Ehlektronika: Mezhvuz. sb. nauch. trudov. Ryazan' (1978) 128-131.
[4] Ya.K. Sklifus, Fazovyye perekhody
teplonositelya v sisteme okhlazhdeniya dizelya teplovoza [Phase transitions of the coolant in diesel locomotive cooling system], Proceedings of Rostov State Transport University. (4(29)) (2014) 92-95.
[5] E. Dilay, J.V.C. Vargas, J.A. Souza, et al., A volume element model (VEM) for energy systems engineering, International Journal of Energy Research. 39 (1) (2015) 46-74.
[6] T.K. Nunes, J.V.C. Vargas, J.C. Ordonez, et al., Modeling, simulation and optimization of a vapor compression refrigeration system dynamic and steady state response, Applied Energy. 158 (15) (2015) 540-555.
[7] J. Catano, T. zhang, J.T. Wenc, et al., Vapor compression refrigeration cycle for electronics cooling, Part I: Dynamic modeling and experimental validation, International Journal of Heat and Mass Transfer. 66 (Nov.) (2013) 911-921.
[8] R. Cabello, E. Torrella, J. Navarro-esbri, Experimental evaluation of the internal heat exchanger influence on a vapor-compression plant performance using R134a, R407C and R22 as working fluids, Applied Thermal Engineering. 24 (13) (2004) 1905-1917.
[9] D.L. Karelin, V.M. Gureyev, V.L. Mulyukin, Modeling system of cooling with two-phase compression installation, Vestnik KGTU im. Tupoleva. (5) (2015) 5-10.
[10] D.L. Karelin, V.M. Gureyev, One method of calculation of condensing temperature of the working agent for vapour-liquid compression cooling systems, Vestnik KGTU im. Tupoleva. (4) (2016) 20-24.
[11] D.L. Karelin, Metodika rascheta parametrov termodinamicheskogo tsikla parokompressionnoy sistemy okhlazhdeniya [Method for calculating the parameters of the thermodynamic cycle of a vapor compression cooling system]. Trudy Akademenergo. (3) (2017) 23-31.
[12] E. elgendy, J. Schmidt, A. Khalil, M. Fatouh, Modelling and validation of a gas engine heat pump working with R410A for cooling applications, Applied Energy. 88 (11) (2011) 4980-4988.
Received 13.05.2018, accepted 02.11.2018.
[13] Lei zhao, M. zaheeruddin, Dynamic simulation and analysis of a water chiller refrigeration system, Applied Thermal Engineering. 25 (14-15) (2005) 2258-2271.
[14] S. Sanaye, M. Chahartaghi, H. Asgari, Dynamic modeling of gas engine driven heat pump system in cooling mode, Energy. 55 (15 June) (2013) 195-208.
[15] L.M. Rozenfeld, A.G. Tkachev, Kholodilnyye mashiny i apparaty [Refrigerating machines and devices], GITL, Moscow (1960).
[16] N.N. Koshkin, Teplovyye i konstruktivnyye raschety kholodilnykh mashin [Thermal and structural calculations of refrigerating machines], Mashinostroyeniye, Leningrad (1976).
[17] V.V. Klimenko, A generalized correlation for two-phase forced flow heat transfer, International Journal of Heat and Mass Transfer. 31 (3) (1988) 541-552.
[18] 2001 ASHRAE handbook, Fundamentals, American Society of Heating, Refrigerating and Air-Conditioning Engineers, Atlanta (2001).
[19] J.M.S. Jabardo, W.G. Mammani, M.R. Ianella, Modeling and experimental evaluation of an automotive air conditioning system with a variable capacity compressor, International Journal of Refrigeration. 25 (8) (2002) 1157-1172.
[20] M.M. Shah, A general correlation for heat transfer during film condensation in tubes, International Journal of Heat and Mass Transfer. 22 (4) (1974) 547-556.
[21] V.V. Avchukhov, B.Ya. Payuste, Zadachnik po protsessam teplomassoobmena [Book of problems in heat and mass transfer processes], Energoatomizdat, Moscow (1986).
[22] I.G. Chumak, V.P. Chepurnenko, S.G. Chuklin, Kholodilnyye ustanovki [Refrigeration units], 2nd ed., Legkaya i pishchevaya promyshlennost [Light and Food Industry PH], Moscow (1981).
[23] B.N. Maksimov, V.G. Barabanov, I.L. Serushkin, Promyshlennyye ftororganicheskiye produkty [Industrial organofluorine ware], Chemistry PH, Leningrad (1990).
THE AUTHORS
KARELIN Dmitriy L.
Kazan (Volga region) Federal University
18 Kremlyovskaya St., Kazan, Republic of Tatarstan, 420008, Russian Federation [email protected]
BOLDYREV Alexey V.
Kazan Federal University - Naberezhnye Chelny Institute
68/19 (1/18) Mir Ave., Naberezhnye Chelny, Republic of Tatarstan, 423812, Russian Federation [email protected]
GÜREEV Viktor M.
Kazan National Research Technical University named after A.N. Tupolev — KAI 10 K. Marx St., Kazan, Republic of Tatarstan, 420111, Russian Federation [email protected]
BOLDYREV Sergey V.
Kazan Federal University — Naberezhnye Chelny Institute
68/19 (1/18) Mir Ave., Naberezhnye Chelny, Republic of Tatarstan, 423812, Russian Federation [email protected]
© Санкт-Петербургский политехнический университет Петра Великого, 2018
