CC BY
25
УДК 532.5:621.694
Багоутдинова А.Г. - кандидат технических наук E-mail: bagoutdinova@rambler.ru
Казанский (Приволжский) федеральный университет
Адрес организации: 420008, Россия, г. Казань, ул. Кремлевская, д. 18 Вачагина Е.К. - доктор технических наук E-mail: vachaginae @mail. ru
Казанский Научный Центр Российской Академии Наук
Адрес организации: 420111, Россия, г. Казань, ул. Лобачевского, д. 2/3 Золотоносов Я.Д. - доктор технических наук, профессор E-mail: zolotonosov@mail.ru
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1
Математическое моделирование труб с винтовой поверхностью теплообмена
Аннотация
Постановка задачи. Для повышения теплогидродинамической эффективности трубчатых аппаратов необходимо разрабатывать новые конструкции теплообменных элементов, отличающихся высокой тепловой эффективностью, технологичностью и невысокой стоимостью производства.
В связи с этим, авторами предложена серия винтовых труб с высокой степенью развития поверхностей теплообмена.
Результаты. В работе построены математические модели труб с винтовой поверхностью теплообмена различных конфигураций. Параметры, входящие в математические модели, позволяют исследовать и в дальнейшем оптимизировать внутреннюю геометрию труб.
Выводы. Значимость полученных результатов для строительной области состоит в том, что предлагаемые теплообменные элементы устанавливаются в инновационных теплообменных аппаратах, монтируемых в индивидуальных тепловых пунктах. Кроме того, строгое математическое описание труб с винтовой поверхностью теплообмена позволит при производстве предлагаемых теплообменных элементов использовать современные компьютерные технологии.
Ключевые слова: моделирование, винтовая поверхность, винтовая линия, гидромеханика, теплообменный аппарат.
Важным направлением решения задач ресурсо- и энергосбережения является разработка и внедрение новых видов конвективных поверхностей теплообмена, отличающихся высокой тепловой эффективностью, технологичностью и невысокой стоимостью производства.
В различных видах промышленности и ЖКХ одним из важнейших видов технологического оборудования являются теплообменные аппараты конвективного типа, общая масса которых достигает до 40 % массы всего оборудования. Модернизация, ремонт такого оборудования и создание новых производственных мощностей требуют огромных капитальных вложений, поэтому внедрение более эффективных и дешевых теплообменных устройств с учетом их доли в общей массе оборудования приводит к существенной экономии материальных и энергетических ресурсов.
В связи с этим, авторами [1-6] предложена серия инновационных теплообменных элементов с винтовой поверхностью теплообмена. На рис. 1 представлен один из вариантов такого элемента.
Рис. 1. Труба с винтовой поверхностью теплообмена
В работе [7] приведена оценка степени развития поверхности у предлагаемых теплообменных элементов по сравнению с поверхностью гладкотрубного элемента. Показано, что этот показатель для труб с винтовой поверхностью теплообмена при определенных геометрических соотношениях может изменяться от 1,4 до 2.
Применение теплообменных поверхностей при у>1,15 считается целесообразным, поскольку положительно влияет на увеличение теплового потока, передаваемого через стенку и обеспечивает снижение металлоемкости теплообменного аппарата [8]. Следовательно, рассматриваемые винтовые трубы могут быть использованы в качестве эффективных теплообменных элементов трубчатых теплообменников.
Для организации серийного производства предлагаемых теплообменных элементов с использованием современных компьютерных технологий необходимо точное математическое описание геометрии их поверхностей.
В работе [9] рассмотрен универсальный метод построения поверхности, образованной движением непрерывной кривой р вдоль криволинейной направляющей
g: r = r (t).
Согласно этому методу радиус -вектор точки поверхности записывается в виде суммы:
r (t, 5) = r (t) + p (t, 5),
где 5 - произвольный параметр в некоторой плоскости кривой у. Радиус-вектор точки границы поперечного сечения в выбранной плоскости зап исывается в виде:
Р (t, 5) = Pl (t, 5) e1 (t) + Р2 (t, 5) e2 (t), где e1 (t) и e2 (t) - единичные векторы, лежащие в плоскости, пересекающей направляющую кривую в каждой точке, соответствующей параметру t.
Следовательно, _ _
r (t, 5) = r (t) + pi (t, 5) e, (t) + P2 (t, j) e2 (t). (1)
В качестве плоскости, которая пересекает направляющую кривую, наиболее удобно использовать плоскость, нормальную к этой кривой в каждой точке, определяемой параметром t. В качестве ортов e1 (t) и e2 (t) наиболее удобно выбирать вектор главной нормали v (t) и вектор бинормали ff (t) направляющей кривой в точке, соответствующей значению параметра t. В этом случае радиус-вектор точки поверхности представляется в виде суммы r (t, s ) = r (t) + p (t, s), где 5 - произвольный параметр в нормальной плоскости кривой у, а p (t, 5) - радиус-вектор точки границы поперечного сечения в нормальной плоскости, описывается равенством:
Р (t, 5) = Pv (t, 5)V (t) + Pf (t, 5) ff (t).
Тогда радиус-вектор точки поверхности представляется в виде:
r (t, 5) = r (t) + Pv (t, 5) V (t) + pp (t,j) f (t). (2)
Наиболее часто в качестве параметра 5 выбирается полярный угол ф в нормальной плоскости кривой у, отсчитываемый от главной нормали по направлению к бинормали; p (t,j) - соответствующий «полярный радиус».
Тогда: _
r (t, j) = r (t) + p (t, j) cos jV (t) + p (t, j) sin jf (t), (3)
где p(t, ф) - переменный радиус границы поперечного сечения канала.
Единичные векторы касательной t, нормали v и бинормали b вычисляются по формулам:
dr / dt • dr / dt
t =i
b = [t xv ] (4)
|dr / d^ \dt /dí| или по формулам:
• dr / dt - \dr / dt x d2r / dt2 ] r r" r-|
тТ-л, f = ,U-i, v =[f xt]■ (5)
|dr / dt| |dr / dt x d2 r / dt2 J v 7
t = i
Для описания поверхности винтовой трубы, образованной движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на эллиптическом цилиндре выберем в качестве ортов следующие векторы:
е, (t ) =
Г 0 Л 0
V10
■ (t ) =
Г costЛ sin t
V 0 0
а в качестве направляющей кривой - винтовую линию, расположенную на эллиптическом цилиндре:
у:х=Т1С081, z=bt, 0<<пп,
где п - количество витков винтовой линии.
В качестве образующей кривой выберем треугольник, определяемый заданной высотой к, проведенной к основанию длиной ¡=1а+1к и проекциями двух сторон и ¡к на третью сторону, к которой проведена высота в плоскости с ортами е1 и е2 :
где
р: р ^, 5 ) = ^ ^) + / (5 ) ¿2 ^), - ¡л < ^ < ¡к,
/■ ) = |(к( + ¡к)/¡а)^ + к' "¡а <^< (5) {"(к (( + ¡к )/¡к ) + к, 0 < 5 < ¡к.
Тогда уравнение боковой поверхности канала запишется в виде:
5 ) =
bt
+ s
Г 0 Л
V10
Л
+ f (s)
cos t sin t
0
Л
- Id £ 5 £
(6)
(7)
где Ь^+^/^п).
На рис. 2 представлена труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на эллиптическом цилиндре, поверхность которой построена по уравнению (7).
Рис. 2. Труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на эллиптическом цилиндре
Для описания поверхности винтовой трубы, образованной движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на круговом цилиндре выберем в качестве ортов следующие векторы:
Г0А ГcostА
ei(t ) =
0 1
,(t ) =
sin t
0
а в качестве направляющей кривои - винтовую линию, расположенную на круговом цилиндре:
g: x = r0 cos t, y = r0 sin t, z = bt, 0 £ t £ pn. Тогда уравнение боковой поверхности запишется как:
5) =
r0 cos tЛ
r0 sin t
bt
+ s
Г 0 Л 0
V10
Л
+ f (s)
t Л
cos sin t 0
- h £ s£ 4,
(8)
где /(ъ) определяется с помощью (6), а b=(¡d+¡k)/(2ж).
На рис. 3 представлена труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на круговом цилиндре, поверхность которой построена по уравнению (8).
Рис. 3. Труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на круговом цилиндре
Запишем уравнение поверхности винтовой трубы, образованной движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на эллиптическом конусе. Пусть у - угол конусности. Выберем в качестве ортов следующие векторы:
м 'tgy cos tЛ M ' cost N
ei (t )= tgy sin t и e2 (t) = sin t
, i 0 ,-tgy 0
а в качестве направляющей кривои винтовую линию, лежащую на эллиптическом конусе: g: х = (r + bttgy)cos t,y = (r2 + bttgy)sin t,z = bt,0 < t < pn. Тогда уравнение боковой поверхности канала запишется как: ((ri + bttgy) cos Л
s
tgy sin t 1
s ) =
(r2 + bttgy )sin t bt
f (s)
cost
л/1+tg
У
sin t
- tgy.
где /(з) определяется с помощью (6), а Ь = (( + ¡к) / ().
Запишем уравнения поверхности трубы в координатном виде:
(у + / (з))
- ld < s < lk,
х (t, s ) = У (t, s ) =
( + bttgy)
(r2 + bttgy) z (t, s ) = bt +
Ф + tg2y + (stgy + f (s ))
Vi+tg2y
s - f (s )tgy
л/1+tg1
cos t;
sin t;
y
Для удобства введем новый параметр О:
Следовательно,
в = bt +
t = в / b -
s - f (s)tgy . Vi+tg2y
s - f (s) tgy bj 1 + tg:
y
Тогда
х (в, s ) = У (в, s) =
( ( r + b
>/b -
s - f (s )tgy
W i+tg:
y
r2 + b
4b -
s - f (s )tgy by¡1 + tg2y
Л
tgy
л Л
tgy
(stgy + f (s))
>/i + tg2y (stgy + f (s))
>/í+ig
y
>/b -
1/b -
s - f (s )tgy
W i+tg:
y
s - f (s)tgy by¡ i + tg2y
z (t, s ) = t
(9)
Здесь 0<з<4+4, 0<в<+1к.
На рис. 4 представлена труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на эллиптическом конусе.
Л11ЦЦЩ|
Рис. 4. Труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на эллиптическом конусе
Для описания поверхности винтовой трубы, образованной движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на круговом конусе выберем в качестве ортов векторы:
(tgy cos t •• ' costл
tgy sin t и e2(t) = sin t
v 1 0 v-tgy 0
e1(t ) =
а в качестве направляющей кривои винтовую линию, лежащую на круговом конусе: g: х = (r0 + bttgy) cos t, y = (r0 + bttgy) sin t, z = bt,0 < t < p n. Здесь у - угол конусности.
Тогда уравнение боковой поверхности трубы запишется как: (( r0 + bttgy) cos tЛ
-(t, * ) =
(r0 + bttgy )sin t bt
' tgy cos t Л , f(*) 'cost Л
* tgy sin t , 1 0 sin t
л/1+tg2y + tg2y
- h < * < h
где /(я) определяется с помощью (6), а Ь = (( + 1к)/ (+ ).
Запишем уравнения поверхности трубы в координатном виде:
( яЩу + / (я ))
х(t, *) = У ( * ) =
ф + tg2y * - f (*)tgy Vi + tg2y
Для удобства введем новый параметр О:
(Г, + bttgy ) +
(Г, + bttgy )-z (t, * ) = bt + -
Vi + tg2y ,
(*tgy + f (* ))
cos t;
sin t;
в = bt +
Следовательно,
t = в / b -
* - f (*)tgy.
л/i+tgy
* = *
* - f (*)tgy .
W1 + tg2y
Тогда
х (в, * ) =
У (в, * ) =
( (
+ b
v v ( (
Г, + b
v v
в / b - *- f(*)tgy 41+tg2y
в / b -
* - f (*)tgy
W1 + tg2y
Л Л tgy 0 0 + (*tgy+f (*))
л/1 + tg2y
Л Л tgy 0 0 + (*tgy + f (*))
Ф + tg2y
в / b-
b + tg2y .
в / b-
* -f(*)tgy bi + tg2y
?(t, * ) = в.
Здесь 0<з<Шк, 0<в<Шк.
На рис. 5 представлена труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на круговом конусе.
Рис. 5. Труба с винтовой поверхностью теплообмена, образованная движением треугольника в плоскости осевого сечения вдоль винтовой линии, расположенной на круговом конусе
Одной из наиболее современных технологий изготовления винтовых труб является метод ротационной ковки (или ротационного обжатия), реализуемый на специальных ротационно-ковочных машинах, где уравнения поверхностей используются при работе аппаратно-программного комплекса, обеспечивающего непрерывное производство труб [10].
На рис. 6 изображен теплообменный элемент, полученный методом ротационной ковки.
Рис. 6. Теплообменный элемент, полученный методом ротационного обжатия
Предлагаемый класс теплообменных элементов разработан для модернизации и реконструкции существующего морально устаревшего и физически изношенного парка трубчатого теплообменного оборудования с помощью замены гладкотрубных теплообменных элементов винтовыми трубами, а также для производства на их базе широкой номенклатуры инновационного энерго- и ресурсосберегающего теплообменного оборудования для различных отраслей промышленности.
Список библиографических ссылок
1. Теплообменный элемент : пат. 119452 Рос. Федерация. № 2012109355/06 ; заявл.
12.03.12 ; опубл. 20.08.12. Бюл. № 23.
2. Теплообменный элемент : пат. 139566 Рос. Федерация. № 2013125072/06 ; заявл.
25.07.13 ; опубл. 20.04.14. Бюл. № 11.
3. Теплообменный элемент : пат. 146241 Рос. Федерация. № 201446289 ; заявл.
20.04.14 ; опубл. 10.10. 14. Бюл. № 28.
4. Теплообменный элемент : пат. 140683 Рос. Федерация. № 2013154131 ; заявл. 05.12.13 ; опубл. 20.05. 14. Бюл. № 14.
5. Теплообменный элемент : пат. 137749 Рос. Федерация. № 2013144123/02 ; заявл. 01.10.13 ; опубл. 27.02.14. Бюл. № 6.
6. Багоутдинова А. Г., Золотоносов Я. Д., Посохин В. Н. Конструкции эффективных теплообменных элементов для скоростных теплообменников // Сантехника. Отопление. Кондиционирование. 2014. № 7. С. 72-75.
7. Багоутдинова А. Г., Золотоносов Я. Д., Посохин В. Н. Математическое моделирование винтовых теплообменных элементов // Известия высших учебных заведений. Строительство. 2014. № 8. С. 41-46.
8. Письменный Е. Н., Баранюк М. М., Вознюк М. М. Равноразвитые поверхности теплообмена и методика численных исследований их теплогидравлических характеристик // Промышленная теплотехника. 2012. Т. 34. № 1. С. 45-54.
9. Багоутдинова А. Г., Вачагина Е. К., Золотоносов Я. Д. Математическое моделирование теплообменных поверхностей пружинно-витых каналов // Известия КГАСУ. 2017. № 3 (41). С. 146-157.
10. Багоутдинова А. Г., Золотоносов А. Я., Золотоносов Я. Д., Яхнев М. Н. Современные теплообменные элементы теплообменного оборудования и технологии их изготовления // Известия КГАСУ. 2013. № 2 (24). С. 150-156.
Bagoutdinova A.G. - candidate of technical sciences E-mail: bagoutdinova@rambler.ru Kazan (Volga region) Federal University
The organization address: 420008, Russia, Kazan, Kremlevskaya st., 18 Vachagina E.K. - doctor of technical sciences E-mail: vachaginae @mail. ru
Kazan Scientific Center of Russian Academy of Sciences
The organization address: 420111, Russia, Kazan, Lobachevsky st., 2/3 Zolotonosov Ya.D. - doctor of technical sciences, professor E-mail: zolotonosov@mail.ru
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya st., 1
Mathematical modeling of pipes with a helical heat exchange surface Abstract
Problem statement. To improve the thermohydrodynamic efficiency of tubular devices, it is necessary to develop new designs for heat exchange elements that are characterized by high thermal efficiency, manufacturability and low production costs.
In this regard, the authors proposed a series of helical pipes with a high degree of development of heat exchange surfaces
Results. Mathematical models of pipes with a helical heat exchange surface of various configurations are constructed in this paper. The parameters included in the mathematical models allow us to investigate and further optimize the internal geometry of the pipes.
Conclusions. The significance of the obtained results for the construction area is that the proposed heat exchange elements are installed in innovative heat exchangers installed in individual heat points. In addition, that a strict mathematical description of pipes with a helical heat exchange surface will allow using modern computer technologies when manufacturing the proposed heat exchange elements.
Keywords: modeling, helical surface, helix, hydromechanics, heat exchanger.
References
1. Heatexchange element : patent 119452 of the Rus. Federation. № 2012109355/06 ; decl.
12.03.12 ; publ. 20.08.12. Bull. in № 23.
2. Heatexchange element : patent 139566 of the Rus. Federation. № 2013125072/06 ; decl.
25.07.13 ; publ. 20.04.14. Bull. in № 11.
3. Heatexchange element : patent 146241 of the Rus. Federation. № 201446289 ; decl.
20.04.14 ; publ. 10.10.14. Bull. in № 28.
4. Heatexchange element : patent 146241 of the Rus. Federation. № 2013154131 ; decl. 05.12.13; publ. 20.05.14. Bull. in № 14.
5. Heatexchange element : patent 137749 of the Rus. Federation. № 2013144123/02 ; decl. 01.10.13 ; publ. 27.02.14. Bull. in № 6.
6. Bagoutdinova A. G., Zolotonosov Ya. D., Posokhin V. N. Designs of effective heatexchange elements for high-speed heat exchangers // Santehnika. Otoplenie. Conditionirovanie. 2014. № 7. P. 72-75.
7. Bagoutdinova A. G., Zolotonosov Ya. D., Posokhin V. N. Mathematical modeling of helical heat exchange elements // Stroitelstvo. 2014. № 8. P. 41-46.
8. Pismennyi E. N., Baranyuk M. M., Voznyuk M. M. Equalized heat transfer surfaces and the technique of numerical studies of their thermo hydraulic characteristics // Promischlennaya teplotechnika. 2012. V. 34. № 1. P. 45-54.
9. Bagoutdinova A. G., Vachagina E. K., Zolotonosov Ya. D. Mathematical modeling of heat-exchange surfaces of spring-twisted channels // Izvestiya KGASU. 2017. № 3 (41). P.146-157.
10. Bagoutdinova A. G., Zolotonosov A. Ya., Zolotonosov Ya. D., Yakhnev M. N. Modern heat-exchange elements of heat-exchange equipment and technologies for their production // Izvestiya KGASU. 2013. № 2 (24). P. 150-156.
