УДК 697.922.564
Давыдов А.П. - кандидат технических наук, доцент E-mail: [email protected]
Валиуллин М.А. - кандидат технических наук, доцент E-mail: posohin@kgasu. ru
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1
Методика расчета всасывающих воздухопроводов с присоединением расхода по пути (Метод коэффициентов расхода)
Аннотация
Постановка задачи. Целью работы является разработка методики расчета воздухоприемных устройств различной геометрии при стационарной и не стационарной режимах работы.
Результаты. При удалении вредных выделений от технологического оборудования используются воздухоприемники различной конструкции. Высокая эффективность работа таких устройств зависит от точности расчета указанных устройств, обеспечивающих заданный объем удаляемого воздуха. Регулирование расхода воздуха по площади воздухоприемника осуществляется изменением коэффициента расхода, который определяется предлагаемой методикой расчета.
Выводы. Значимость для строительной отрасли заключается в возможности рассчитать воздухоприемное устройство любой геометрии при стационарном и не стационарном режимах работы.
Ключевые слова: методика расчета всасывающих воздухопроводов, равномерно-всасывающий воздухопровод, отверстие, внутренний экран, коэффициент расхода.
В вентиляционной практике для расчета всасывающих воздухопроводов наиболее распространенными являются методики В.Н. Посохина [1, 2]. В известной мере это объясняется достаточной упрощенностью этих методов расчета. Упрощения приводят на практике к существенным отклонениям от поставленных задач, а также завышенным энергетическим затратам (большое сопротивление всасывающих воздухопроводов, завышенные объемы удаляемого воздуха). Специальные исследования по оптимизации конструкций всасывающих воздухопроводов отсутствуют. Цель работы - разработка методики расчета воздухоприемных устройств различной геометрии при стационарной и нестационарных условиях работы.
Рассмотрим собирающий воздухопровод произвольной геометрии (рис. 1). Пусть направление движение воздуха в воздухоприемнике произвольное вдоль осей ху2 со скоростью W. Вход воздуха проходит непрерывно через стенки воздухоприемника в плоскостях ХОУ (4 -1-3-3 ), ХО2 (1-4-4 -1 ), УО2 (1-2-3-4) соответственно со скоростью Уху, Ух2, Ууг. Стенки воздухоприемника бесконечно тонкие. Каждая грань воздухоприемника имеет собственную проницательную способность характеризуемую коэффициентом расхода на всасывание Ци,^,, ^и.хг, ^и.уг и на истечение ^в.хг, ^в.уг. Скорость истечения соответственно Жх, Жу, Ж2. Течение воздуха рассматривается как течение не сжижаемой жидкости. Течение изотермическое. Скорость всасывания принимается V; скорость в воздухопроводе - Жху2.
Основной задачей является - рассчитать воздухопровод таким образом, чтобы удельные расходы по площадям всасывания соответствовали поставленным требованиям. В общем случае удельный расход воздуха через грани воздухоприемника при нестационарной работе можно записать в виде:
Я,=/(х,у,г,0. (1)
Или в соответствии плоскостям всасывания:
дхо2=/(х;2;0; дуо^Ду^^); Цхоу=/(х;у;Г). (2)
Рис. 1. Расчетная схема воздухопровода (схема течения)
Стационарное всасывание воздуха
При стационарной работе воздухоприемного устройства вход воздуха проходит через соответствующие грани. По оси х грань 1-1'-4'-4 со скоростью входа Ухг и переменным коэффициентом расхода ¡в.х.2. По оси у грань 1-2-3-4 со скоростью входа Ууг и переменным коэффициентом расхода ¡ву.2. По оси 2 грань 4-3-3'-4' со скоростью входа Уху и переменным коэффициентом расхода ¡вх ху.
Истечение воздуха происходит через соответствующие грани. По оси х через грань 1'-2'-3'-4' скорость истечения Жх коэффициент расхода истечения х. По оси у через грань 2-2'-3'-3 скорость истечения Жу коэффициент расхода истечения ¡и.у. По оси 2 через грань 1-1'-2'-2 скорость истечения Ж2 коэффициент расхода истечения ¡¡и.2.
При неравномерном заборе воздуха по плоскостям всасывания скорости всасывания будут равны:
V. =
V =
у2
V =
ху '
дчх. ^ 2 + +
дх 0 к д.2
2 + (дЧу,
ду ) [ &
дЧу Л 2 + (дЧху
дх 1 ду
(за)
(зб)
(зв)
Скорость воздуха во всасывающих плоскостях можно определить по выражению:
(4а)
V = %х №.
х2 ((ж:))р' V = % \Щ.
У2
((су )) р' V = %. р/7
((Ж.)) р'
(4б) (4в)
где фг - коэффициент скорости; ег - коэффициент сжатия струи; Сг - коэффициент местного сопротивления.
Для нашего случая, можно принять, ф = 1, е = 1 учитывая, что стенка имеет
бесконечно малую толщину. Обозначая величину соотношения
локальный коэффициент расхода можем записать:
^ =
Эх 0 I
2СР'
тв.уг =
=
Эу
Эг
А
\2
Эх
Эу
Чу Р .
2СРу '
ар
2СР.'
((с" )
=
как
(5а) (5б) (5в)
Нестационарное всасывание воздуха
Для нестационарного потока скорость входа можно выразить в зависимости от изменения удельного расхода по площади всасывания:
V =
V =
V =
(ба)
(бб) (6в)
В случае нестационарной (квазистационарной) работы всасывающих устройств коэффициент расхода можно определить из выражений:
Двх: =
т у: =
т ху =
СхР
2СР
-с.х
эv а
Эх
\ 2
Сх
аур
Эх
2 /
ЭVydX Эх
2СР
сд: р
2СР
(7а) (7б) (7в)
Выражения (7а^7в) позволяют рассчитать величину коэффициента расхода для любой всасывающей поверхности воздухопровода.
Для определения необходимой величины коэффициента расхода, регулирующего расчетный расход воздуха, необходимо знать распределение удельного расхода по плоскостям всасывания, а также закономерность изменения давления по пути движения воздуха внутри воздухопровода.
Наиболее простой случай работы всасывающего воздухоприемного устройства -воздухопровод постоянного поперечного через одну грань которого равномерно входит поток воздуха.
Рис. 2. Расчетная схема воздухопровода постоянного поперечного сечения, входящий в поток не вносит составляющей количества движения
1
Рассмотрим всасывающий воздухопровод постоянного поперечного сечения /=1у*1г и длиной 1к(х) (рис. 2). Примем, что воздух поступает в воздухопровод непрерывно и равномерно через грань 'х*'2, скорость входящего потока обозначим через Ув.х. Для рассматриваемого случая коэффициент расхода будет равным:
( гг.—^—т"—^
=
| 11 дд*
dx ) у dz
fr (8)
2dPx ■ (8)
При условии, что воздух поступает в воздухопровод непрерывно и равномерно через грань lxxlz, а также ^dqxz/dz = 0 (по оси z отсутствует изменение давления), Увх = const и dqxz/dx = Увх скорость всасывания по оси х в расчетном сечении не меняется и остается постоянной, уравнение (8) примет вид:
mx = (9)
При условии, что у = и 2Р0 и dPx=plx*dfx коэффициент расхода в точке на ю Ц р
расстоянии 1Х находим из выражения:
= (10)
\рх„
Для определения величины ¡лХ необходимо знать величину давления в начале воздухопровода -рю и в расчетном сечении - рь.
Поскольку ось х расположена горизонтально, представляется возможным исключить из учета силы тяжести.
Выделим два сечения 1 и 2, соответственно на расстоянии 1Х и lx+dlx. Допустим одновременно, что давление в плоскости, перпендикулярной оси х постоянно, а распределение давления подчиняется гидростатическому закону.
В общем виде уравнение изменения проекции количества движения системы (11), равное сумме проекций импульсов внешних сил и при отсутствии инерционных сил, можно записать:
I р(К + dWx)] - | ± | рУв2х соз^ = Рх,^ - (Рх + dРx)[/,] ~[{гх)и ]■ (11)
Решим уравнение (11) относительно величины давления учитывая, что:
¡рК ^ = РРКЛ, (12)
где вХ - коэффициент Буссинеско.
После преобразования уравнения (12) и пренебрегая составляющими второго порядка малости, учитывая, что cos 0 = 0 получим:
2рвхКСр.х = dpx - dт■ (13)
Записывая относительно арх, получим:
- dpx = - 2рвхКсР.х- dт. (14)
Здесь представляется целесообразней выразить среднюю скорость в расчетном сечении - Жср.х через конечную среднюю скорость - Жср.к в конце воздухопровода. В нашем случае равномерного присоединения потока по длине:
КР.х = . (15)
1к
Подставляя выражение (15) в уравнение (14) опуская индексы ср. получим:
-йРх = -2рЬК„2 '^Г - dt. (16)
к
Принимая
d я wK23i2 dix
dt = (17)
d 2 l2
экв k
В этом случае уравнение (14) примет вид:
d 2 nR W 2 ixdlx Г я W23lx2dlx
d P 2 lt
экв. k
(18)
В этом уравнении не раскрыта зависимость изменения коэффициента сопротивления трения - X по длине воздухопровода для движения воздуха с переменным по пути расходом. Обычно здесь принимается - X = const, что для приближенных расчетов допустимо, но, конечно, не отражает истинной картины. Основная сложность заключается в том, что для различных рассматриваемых сечений воздухопровода число Re, по которому принято определять величину коэффициента трения, неодинаково.
Определения коэффициента трения в настоящее время становится возможным лишь по эмпирическим формулам. Вполне приемлемо дающая хорошую сходимость с опытными данными для области работы воздухопроводов является формула А. Д. Альтшуля [3]:
K 684 Re
l = 0,111 —^ + —
d эк
Подставив это значение (19) в уравнение (18) имеем:
(19)
dp, =-2pbW
ldl
(
- 0,11
68
+ Re
Wp
3/2 dlx
2d„,
(20)
где k - коэффициент абсолютной шероховатости, м; dKKe.- эквивалентный диаметр, м;
d wj—
Re = -
Ч -
критерий Рейнольдса; 0,11 и 68 - эмпирические коэффициенты.
Записывая уравнение (20) в интегральной форме получим:
J dplx = -J 2pbW2 ^ -J 0,11
к 68
dm Re
/
, wpK 3i2x dix
' 2dK.„ l,2
(21)
Интегрируя уравнение (21) в установленных пределах при условии, что cos О = 0 т.е. входящий поток присоединяется к основному под углом в = 900, а = в = const.
p = Pi0 + PW2 ^+0,11
/ \ 0,25
yd3Ke + Re
Wl l3
ср.к ^x
эке k
2d,
(22)
где
Pl0 =
V2
ех
2
0
2m
p.
Решая совместно уравнение (10) и уравнение (22) получим:
mlx = m,,
vx
v2 + 2ml W
(i
V h 0
W2l 3 (
+0,11m/
0 2dK„ll2
\ 0,25
к 68 ^
— + —
эе; V dK Re 0
(23)
Выражение (23) позволяет определить коэффициент расхода по длине линейного воздухопровода постоянного поперечного сечения для случая, когда входящий поток воздуха не вносит составляющей количества движения.
Используя рассматриваемую методику, представляется возможным рассчитать воздухоприемное устройство любого назначения и любой геометрии [4, 5].
Список библиографических ссылок
1. Посохин В. Н., Сафиуллин Р. Г., Бройда В. А. Вентиляция. М. : АСВ, 2015. 624 с.
2. Посохин В. Н. Аэродинамика вентиляции. М. : АВОК-ПРЕСС, 2008. 209 с.
3. Альштуль А. Д. Гидравлические сопротивления. М. : Недра, 1970. 248 с.
4. Давыдов А. П., Валиуллин М. А. Уравнение движения воздуха в кольцевом равномерно-всасывающем воздухопроводе переменного сечения // Известия КГАСУ. 2011. № 2 (16). С. 142-146.
5. Давыдов А. П., Ланцов А. Е., Ахмерова Г. М. Работа равномерно-всасывающей панели с отверстиями круглого сечения // ООО «Скрипта». 2013. С. 461-465.
6. Давыдов А. П., Валиуллин М. А. Коэффициент расхода всасывающей щели с внутренним экраном // Известия КГАСУ. 2011. № 3 (17). С. 85-89.
7. Давыдов А. П., Ланцов А. Е., Ахмерова Г. М. Зависимость коэффициента расхода иссасывающего отверстия круглого сечения от числа Яе: сб. докладов международной научно-технической конф. «Теоретические основы расчета систем ТГВ» / МГСУ. М., 2013. С. 201-204.
0
8. Усманов Р. А., Давыдов А. П., Валиуллин М. А. Расчет воздухопроводов равномерной раздачи и всасывания с внутренними экранами, установленными под углом // Вестник КТУ. 2014. № 12. С. 138-140.
9. Посохин В. Н., Зиганшин А. М., Варсегова Е. В. К расчету потерь давления в местных сопротивлениях. Сообщение 3 // Известия высших учебных заведений. Строительство. 2016. № 6. С. 58-65.
10. Kumada T., Hirota T., Tamura N., Ishiguro R. Heat and mass transfer with liquid evaporation into a turbulent air stream // Letters in Heat and Mass Transfer. 2012., № l. P. 1-9.
11. Volchkov E. P. Concerning the heat and mass transfer features on permeable surfaces // Inter. J. Heat Mass Transfer. 2013. P. 755-760.
Davydov A.P. - candidate of technical sciences, associate professor E-mail: [email protected]
Valiullin M.A. - candidate of technical sciences, associate professor E-mail: posohin@kgasu. ru
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya st., 1
Method of calculation of suction air pipes with flow connection along the way
(Method of flow coefficients)
Abstract
Problem statememt. The aim of the work is to develop a methodology for calculating air-intake devices of various geometries for stationary and non-stationary operation modes.
Results. When removing harmful emissions from the process equipment used air inlets of various designs. The high efficiency of such devices depends on the accuracy of the calculation of these devices, providing a given volume of air being removed. The air flow rate is regulated by the air intake area by changing the flow coefficient, which is determined by the proposed calculation procedure.
Conclusions. Significance for the construction industry is the ability to calculate the air receiver device of any geometry for stationary and non-stationary modes.
Keywords: calculation technique of suction air ducts, uniformly suction air line, aperture, internal screen, flow coefficient.
References
1. Posohin V. N., Safiullin R. G., Broida V. A. Ventilation. M. : ASV, 2015. 624 p.
2. Posohin V. N. Aerodynamics of ventilation. M. : AVOK-PRESS, 2008. 209 p.
3. Alshtul A. D. Hydraulic resistance. M. : Nedra, 1970. 248 p.
4. Davydov A. P., Valiullin M. A. The equation of motion of air in the ring evenly-suction air duct of variable cross section // Izvestiya KGASU. 2011. № 2 (16). P. 142-146.
5. Davydov A. P., Lantsov A. E., Akhmerova G. M. Work evenly - the suction plate with the holes of circular cross section // OOO «Scripta»2013. P. 461-465.
6. Davydov A. P., Valiullin M. A. Flow Coefficient suction slot with inner screen // Izvestiya KGASU. 2011. № 3 (17). P. 85-89.
7. Davydov A. P. Lantsov A. E., Akhmerova G. M. Dependence of the flow issasyvayuschego holes of circular cross section of the number Re: Proceedings of the International / MGSU. M., 2013. P. 201-204.
8. Usmanov R. A., Davydov A. P., Valiullin M. A. Calculation of the air ducts still -dimensional distribution and absorption of internal screens installed at an angle // Herald KTU. 2014. № 12. P. 138-140.
9. Posokhin V. N., Ziganshin A. M., Varsegova E. V. The calculation of the pressure losses in local resistances. Message 3 // Izvestija vysshih uchebnyh zavedenij. Stroitel'stvo. 2016. № 6. P. 58-65.
10. Kumada T., Hirota T., Tamura N., Ishiguro R. Heat and mass transfer with liquid evaporation into a turbulent air stream // Letters in Heat and Mass Transfer. 2012., № l. P. 1-9.
11. Volchkov E. P. Concerning the heat and mass transfer features on permeable surfaces // Inter. J. Heat Mass Transfer. 2013. P. 755-760.