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

Ill

Tt: ir i iif ii i к ii IE к. hi : 11 [ ii. i ir 11 ii ai. i;f i II. м 111 iimi и 1а»н Ш11 к ио.цллл, ■ л. ни1 иль<1>i:■■ ■■ i. н ni нi:iц ьни i:

УДК 532.5:621.694

Вачагина Е.К. - доктор технических наук E-mail: vachaginae @mail. ru

Казанский Научный Центр Российской Академии Наук

Адрес организации: 420111, Россия, г. Казань, ул. Лобачевского, д. 2/3

Золотоносов Я.Д. - доктор технических наук, профессор

E-mail: [email protected]

Мустакимова С.А. - ведущий программист

E-mail: mustakim@kgasu. ru

Крутова И.А. - аспирант

E-mail: iraida_knyazeva@mail. ru

Казанский государственный архитектурно-строительный университет

Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1

Определение эквивалентного диаметра пружинно-витых труб Аннотация

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

Ключевые слова: теплообмен, гидродинамика, пружинно-витые каналы, ковариантные и контравариантные компоненты скорости, линии тока, вторичные течения.

В литературе для определения основных гидравлических характеристик каналов используется параметр, известный в литературе как эквивалентный диаметр[1]. Так для круглых труб эквивалентный диаметр равен обычному диаметру, для кольцевого канала -разности диаметров. Представляет научный и практический интерес определение эквивалентного диаметра пружинно-винтового канала в целях их широкого внедрения в инженерную практику [2-7].

Границы пружинно-витых каналов образованы плотно примыкающими друг к другу свернутыми по спирали витками проволоки [6, 7].

Параметрическое уравнение одного витка центральной винтовой линии - оси проволоки имеет вид [8]:

р(и) = (D/2)(e(u) + wuk),0 < и < 2p ,

или в координатном виде:

(d / 2)cosи, y = (d /2)sin и, ,

z = (d /2)wu, 0 < и < 2p где и - параметр; D/2 - расстояние от центральной винтовой линии - оси проволоки до центра трубы; w=S/(nD) - геометрический параметр; S - расстояние между двумя витками проволоки.

Центральная винтовая линия оси проволоки характеризуется следующими векторами сопровождающего репера:

х = i

т {и ) =

g {и)+wk 1

w- +1

b {и )

- wg {и)+k yjw2 +1

v (u) = -e (u) ,

где e (u) = cosui + sinuj , g(u) = -sinui + cosuj - круговые векторные функции скалярного аргумента u; i, j, k - орты декартовой системы координат.

Параметрическое уравнение поверхности свернутой по спирали проволоки диаметром d может быть представлено в виде:

r (u, v) = (D / 2)(e (u) + wuk ) + (d / 2)(v (u) cos v + b (u) sin v) j

или в координатном виде:

D

x = — 2

(

D

cos u + K\ - cos vcos u + Sin v

WSin u

4

w + 1

У = -

sin u + к

V

-w cos u

- cos v sin u + sin v

V

sin v

w2 + 1

J 0

Б

2 = —I wu + к- I-

21 4мГ+1

где и, V - параметры поверхности; к = ё / Б .

Полагая 2=0 в третьем соотношении (1) выразим из него и:

sin v

u = -к-

(1)

Х^+7- (2)

w^/w +1

Подставляя соотношение (2) в первые два уравнения системы уравнений (1) получим параметрическое уравнение границы канала в плоскости 2=0:

D

x = — 2

(l - к cos v)cos

^ к sin v ^

y =

D 2

V wVw" + 10

^ к sin v ^

2

к w

^ к sin v ^

(l - к cos v)sin

4

w2 + 1

-sin v sin

V

V wyfw" + 10

к w

4

sin v cos

w + 1

V wVw2 + 1 00 ^ к sin v ^

V wVw2 + 1 0y

-p < v < P

(3)

Условие плотного прилегания витков проволоки друг к другу можно представить в виде системы уравнений:

У

(v0) = D (1 - кcosvo)sin

^ к sin v0 ^

V wyfw' +10 f

к w

^ к sin v„ ^

w + 1

V wyfw* + 1 00 \ Л

= o,

dy ( ) D к ( + 1 - к cos v0) dv 2 w (w2 +1)

4

w +1 cos v„ cos

к sin vn

V w^fw" + 10

2

+w sin v0 sin

к sin v„

V w\Zw" +10

= 0,

(4)

где ±у - значение параметра, при котором витки плотно прилегают друг к другу. Смысл

первого уравнения состоит в том, что витки проволоки плотно примыкают друг к другу в точке с координатой у=0. Смысл второго уравнения состоит в том, что в точке у=0 функция у(V) принимает экстремальное минимальное значение, что соответствует условию ёу/^=0.

Система (4) с помощью некоторых преобразований может быть приведена к виду:

Л

2

к sin v„

к cos v„ = sin

" V wVW + 1 0 2 / 2 \ 2 к cos vQ -(w + 1)cos V0 + KW = 0,

Второе уравнение системы (5) является квадратным уравнением относительно cosv0. Допустимое решение этого уравнения запишется как:

(w2 +1) -yj(w2 +1)

(5)

.22 - 4к w

cos v0 =

(6)

2к

при этом для sinv0 имеем следующее выражение:

1 -

(w2 + l) + ^(w2 +1)2 - 4k

2 2 w

Тогда первое уравнение системы (5):

. 2

k cos v„ = sin

2k

k sin v„

V wVw +10

совместно с полученными выражениями для со8У0 и 81пу0 можно рассматривать как нелинейное уравнение относительно w. Решая это уравнение с помощью численных методов [9], например, с помощью метода деления отрезка пополам, получим зависимость w0=w0(к), которая определяет конкретное значение геометрического параметра w0 (значение шага винтовой спирали) в зависимости от параметра к при плотном примыкании витков проволоки друг к другу. При известном значении w0 значение v0 получим из (6) с помощью формулы:

v0 = Arccos

(w02 + l)-i(w02 + l)2 - 4k

2 2 w

0

2k

так как из формул очевидно следует, что это значение удовлетворяет неравенству

0<у0<к/2.

На рис. представлены поперечные сечения пружинно-витого канала.

-»■ ..... чч. 4 Ч/ч

0.5 \\

* ' У \ % \

// / \ * *

!* / \ \ \

ilTv i» \ 0.5 0I5 1 !

/ '

/ ¡ !

'Л / / ;

\ s 0.5 —/ /

...... ...... .....

а б

Рис. Поперечное сечение пружинно-витых труб: к=1/2 - сплошная линия; к=1/5 - штриховая линия; к=1/10 штрих-пунктирная линия; а - поперечное сечение проволоки; б - поперечное сечение проточной части канала

Для определения площади поперечного сечения используется стандартная формула в полярной системе координат [10]:

п П<Р) 1 п

F = / dj / Ыг = — / Г \<р)dj .

-п о 2 -п

Учитывая, что уравнение границы задано в параметрическом виде (3) преобразуем предыдущую формулу к виду:

v

F = 1 j0 г 2 (v) j dv 2 -v 4 ' dv

¿ 0

(7)

Для получения dq/dv воспользуемся формулами перехода от декартовой к полярной системе координат:

ÍX = г cos j y = г sin j '

Дифференцируя эти соотношения по параметру V выразим из полученных соотношений:

ё? ^ = (_ ¡о (8)

¿¡V Зч г V ск ёу 0 Производные —, — определим дифференцируя по V параметрическое уравнение

¡¡V ¡V

границы (3):

ёх Б к (к соз V-1-w2)

¡V 2

w (1+м>2)

(

к 81П V

- w соз 1

V ч w\w2 +1 0

81П V +

+ \ 1 + w з1п

Л

к з1п V

V +1 0

ёу Б к (к соз V-1-w2)

Л

0

¡V 2 w (1+w2)

г \ \

I

1+w2 соз

( . \ к з1п V

Ww +1

соз V+

к з1п V

1

V\1 + w 0

0

Подставляя (9) в (8) и в (7) окончательно получим:

1 ( Б 2 V) ( ёу ёх Л

— — х— - у—

2 V 2 0 - V ч ёv

ёV

где подынтегральное выражение имеет вид:

(

^ ёу ёх^ х--у —

V ёV ёV 0

4+4w04 + 3к 2 + 4w02'

(9)

(2+к2)) соз V+к -4(1+w02) -4(1+w02) cos2v+к cos3v

V 00

4 w0

1+w0

получено подстановкой выражений для х и у из (3), выражений — и ёу - из формул (9).

ёv ёv

Точное значение площади поперечного сечения канала может быть только численно, так как только численно могут быть определены пределы интегрирования. Однако, как показали расчеты значение Vo мало отличается от к/2. Заменяя пределы интегрирования на к/2, принимая Б/ё»1 и w0»к/к, можно получить приближенное значение площади Б:

Б :

Б

5 3 3 4 -3р к - 6р к + 6рк - 3кк

, 2 2Б 2\ 4/, . 2\ 6р к 12 + к 1 + р 16 + 4к i

3

( 2 2\2 |р + к I

Смоченный периметр Р может быть вычислен как длина параметрически заданной кривой по формуле [1]:

Р =

Б \ V

2 )-\

| ^/(Х/Л) + (у / dv) dv .

Подставляя выражения из (9) в формулу для периметра, получим:

Р =

( Б ^

V 2 )

к (1 + -М 2 - к соз V )

^—^—

-V

('+w:)

,Jl+2w^+cos(2v)ёv '

<

3

и

Этот интеграл также определяется только численно.

Для определения приближенного значения смоченного периметра при тех же допущениях, что и при вычислении площади живого сечения можно использовать формулу:

f i^ ^ 1 +1 — I - к cos v

к к

P"(7)!-

4 У 2

+(к)')

1+2 (7

+ cos

(2v)dv.

Значение интеграла в этой формуле также нельзя выразить в элементарных функциях даже для конкретных значений пределов интегрирования, так как он сводится к эллиптическим интегралам и может быть получен только численно. Эквивалентный диаметр рассчитывается по формуле:

D =

экв

4 F

P

= 2

J

-v

dy dx

x--y —

V dv dv

dv

2

vn к (1 + w - к cos v

w

w - к cos v) I-

T-+ 2 w02 + cos (2v )dv

(1 + w0 )

Аналогично тому, как это делалось при вычислении площади живого сечения и смоченного периметра можно получить приближенную формулу для вычисления эквивалентного диаметра:

D,

экв

к-

-3к5к - 6к3к3 + 6кк4 - 3кк5 + 6к2 к2 (2 + к2) + к4 (6 + 4к2)

V 2 У

1/2 2

WT

i з/ 2 2Y 2 Ък \л + к J p

1 +

2

к 2 V

- к cos v

V к

1 + 2

к к

+ cos

(2v)dv

Были рассчитаны конкретные значения площади живого сечения, смоченного периметра и эквивалентного диаметра, которые представлены в табл.

Таблица

Характеристики поперечного сечения пружинно-витого канала

к=1/2 к=1/5 к=1/10

Vo 1,55805 1,56999 1,57070

wo 0,16121 0,06379 0,03185

S/d 1,01294 1,00203 1,00051

F/(D/2)2 1,2064 2,2428 2,6704

F/(D/2)2 1,2015 2,2428 2,6704

P/(D/2) 4,1029 5,35507 5,80667

P/(D/2) » 4,11562 5,35539 5,80669

DK/(D/2) 1,17614 1,67525 1,83954

Dk/7/2) » 1,17227 1,67515 1,83953

Заключение

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

Список библиографических ссылок

1. Касаткин А.Г. Основные процессы и аппараты химической технологии. - М.: Химия, 1971. - 782 с.

2. Золотоносов Я.Д., Князева И.А. Изготовление теплообменных аппаратов нового поколения на базе пружинно-витых каналов // Тезисы докладов XV Международного симпозиума «Энергоэффективность и энергосбережение в республике Татарстан». - Казань, 2015.

3. Евсеев Е.С., Золотоносов А.Я., Золотоносов Я. Д. Высокоэффективные теплообменные аппараты на базе теплообменных элементов в виде пружинно-витых труб // Труды АкадемЭнерго. - М., 2008, № 4. - С. 18-33.

4. Басова О.А., Золотоносов Я.Д. Змеевиковый теплообменник на базе пружинно-витых каналов // Сборник научных трудов. - Казань. 2009. - С. 33-35.

5. Князева И.А., Золотоносов Я.Д. Змеевиковый теплообменник и способ его изготовления // Тезисы доклада 66 Всероссийской научной конференции по проблемам архитектуры и строительства. - Казань: КГАСУ, 2014. - 173 с.

6. Патент РФ № 62694 на полезную модель. Теплообменный элемент / Золотоносов А.Я., Золотоносов Я.Д., Конахина И.А. № 2006143517/22, заявл. 07.12.06; опубл. 27.07.07; Бюл. № 12.

7. Патент РФ № 64750 на полезную модель. Теплообменный элемент / Золотоносов А.Я., Золотоносов Я.Д. № 2007107173/22, заявл. 26.02.07; опубл. 10.07.07; Бюл. № 19.

8. Норден А.П. Краткий курс дифференциальной геометрии. - М.: Физматгиз, 1958. -244 с.

9. Бахвалов Н.С. Численные методы. - М.: Наука, 1975. - 631 с.

10. Кудрявцев И.А. Курс математического анализа. - М.: Высшая школа, 1981. - 687 с.

Vachagina E.K. - doctor of technical sciences E-mail: vachaginae @mail. ru

Kazakh Scientific Center of the Russian Academy of Sciences

The organization address: 420111, Russia, Kazan, Lobachevsky st., 2/3 E-mail: vachaginae @mail. ru

Zolotonosov Ya.D. - doctor of technical sciences, professor

E-mail: [email protected]

Mustakimova S.A. - lead programmer

E-mail: mustakim@kgasu. ru

Krutova I.A. - post-graduate student

E-mail: iraida_knyazeva@mail. ru

Kazan State University of Architecture and Civil Engineering

The organization address: 420043, Russia, Kazan, Zelenaya st., 1

The definition of equivalent diameter of spring-twisted pipes

Resume

The method of calculating the area of the living section, wetted perimeter and equivalent diameter of the cylindrical spring twisted channel was presented. The inner boundary of such channels is formed closely adjacent to each other coils coiled in a spiral with some constant pitch of the wire. For this the equation of the surface of wire twisted in a helical spiral was written in parametric form. Was formulated in terms of tight junction turns of wire to each other. These conditions are reduced to the dependence of one geometric parameter from the other geometric parameter. It makes sense the relationship of the step helix to the diameter of the wire, and has a sense of the relationship of the wire diameter to the diameter of the axial bending of the Central helix of wire. The required dependence is obtained from the search of the root of the nonlinear equation in the where parameter can be considered as a parameter. To search root used numerical methods to find the root of nonlinear equation, including the method of halving. The equation of the boundary of the cross section was obtained as the equation of a curve in cross-section of his plane. Thus, the equation of this curve is the border was also presented in a parametric form. To compute the area of the living section was used standard formulas to find the area of a curvilinear sector in polar coordinates. The formula for calculating the cross sectional area has been transformed considering the fact that the curve bounding the cross section is given in

parametric form. To calculate the wetted perimeter was used the standard formula of mathematical analysis to determine the length of a curve given parametrically. Calculation of equivalent diameter using the proposed method allows producing engineering calculations series of heat exchangers in which heat exchange elements are spring-twisted channels.

Keywords: heat transfer, fluid flow, and spring-twisted channels, covariant and contravariant components of velocity, streamlines, secondary flows.

Reference list

1. Kasatkin A.G. Basic processes and apparatuses of chemical technology. - M.: Chemistry, 1971. - 782 p.

2. Zolotonosov Ya.D., Knyazeva I.A. Manufacturing of heat exchangers of a new generation on the basis of spring-twisted channels // Proceedings of XV International Symposium «Energy Efficiency and conservation in the Republic of Tatarstan». - Kazan, 2015.

3. Evseev E.S., Zolotonosov A.Ja., Zolotonosov Ya.D. Highly efficient heat transfer devices based on heat-exchange elements in the form of spring-twisted pipe // transactions of AcademEnergo. - M., 2008, № 4. - P. 18-33.

4. Basova O.A., Zolotonosov Ya.D. Coil a heat exchanger on the basis of spring-twisted channels // Collection of scientific works. - Kazan, 2009. - P. 33-35.

5. Knyazeva I.A., Zolotonosov Ya.D. Coil a heat exchanger and method of its making // Abstracts of the 66 all - Russian scientific conference on problems of architecture and construction. - Kazan: KSUAE, 2014. - 173 p.

6. RF patent № 62694 utility model. The heat exchange element / Zolotonosov A.Ja., Zolotonosov Ya.D., Konahina I.A. № 2006143517/22, Appl. 07.12.06; publ. 27.07.07; bull. № 12.

7. RF patent № 64750 utility model. The heat exchange element / Zolotonosov A.Ja., Zolotonosov Ya.D. № 2007107173/22, Appl. 26.02.07; publ. 10.07.07; bull. № 19.

8. Norden A.P. A Short course in differential geometry. - M.: Fizmatgiz. 1958. - 244 p.

9. Bakhvalov N.S. The Numerical methods. - M.: Science, 1975. - 631 p.

10. Kudryavtsev I.A. A course of mathematical analysis. - M.: Higher school. 1981. - 687 p.