student
VAC 05.17.00 Химическая технология
VAC 05.23.00 Строительство и архитектура
UDK 69 Строительство. Строительные материалы. Строительно-монтажные работы
Based on fixed discrete vortex modeling of detached flow at the entry in the round suction duct in the range of action is realized. Fine circular screen, drop supplied is found in the action range of suction duct. Calculated coefficients are compared with received experimental dates. Bond formula between localized resistance’s coefficient of round inlet duct and coefficient of flow‘s compression is proposed. Results of experimental research of detached flows at the entry of slit-like suction duct, different mechanical screen s supplied are applied. Screens axes which are promoting maximal air cost reduction are determinate.
detached flow, discrete vortex, mechanical screen, slit thinness
Введение. С отрывными течениями приходится сталкиваться во всех технических отраслях, имеющих дело с течением жидкости или газа [1–3]. В промышленной вентиляции отрыв потока изучают на входе во всасывающие отверстия местных вентиляционных установок
[4–14]. Профилирование по границам отрывных областей на входе во всасывающие каналы открытых местных вентиляционных отсосов позволяет снизить их аэродинамическое сопротивление, что приводит снижению энергоемкости систем вентиляции [4, 5]. Для снижения расхода воздуха, поступающего через неплотности аспирационных укрытий, имеет смыл напротив повысить аэродинамическое сопротивление входа во всасывающие каналы. Данный подход продемонстрирован в работах [9–13]. Целью данной статьи является выявление путем вычислительного и натурного экспериментов закономерностей отрыва потока на входе во всасывающие каналы круглой и щелевидной формы, а также определение параметров механических экранов, способствующих наибольшему снижению расхода воздуха, поступающего через неплотности аспирационных укрытий.
Численное и экспериментальное исследование отрыва потока на входе во всасывающие каналы с механическими экранами
Рассматривается область течения на входе в круглый всасывающий канал, в спектре действия которого находится тонкий кольцевой экран, снабженный выступом (рис. 1). Моделирование отрыва потока осуществляется на основе стационарных дискретных вихрей.
Рис. 1. Область отрывного течения: а) граница области; б) дискретизация границы области
С острых кромок А, В, С происходит срыв потока и образуются свободные поверхности тока. Необходимо определить их положение, скорость потока в любой заданной точке и коэффициент местного сопротивления при входе во всасывающее отверстие. Математическая постановка задачи состоит в решении двумерного уравнения Лапласа для потенциальной функции j: , при заданных значениях граничной нормальной составляющей скорости
, где x – точка границы S. Функция
выражает влияние свободных вихрей, находящихся на свободных поверхностях тока, положение которых заранее неизвестно.
Данное уравнение сводится к граничному сингулярному интегральному уравнению:
,
где – плотность циркуляции присоединенного вихревого слоя;
− плотность циркуляции свободного вихревого слоя, непрерывно размещенного на поверхности
, образующейся при срыве потока с острой кромки;
– произвольная точка границы S. Функция
численно равна величине скорости в точке
вдоль направления единичного вектора n
, вызываемой вихрем единичной циркуляции, расположенной в точке
.
Для осесимметричной задачи в цилиндрической системе координат:
,
,
,
,
,
.
Вихревой слой, моделирующий поверхность трубы и экрана, заменяется бесконечно тонкими вихревыми кольцами постоянной интенсивности ,
(рис. 1, б). Вихри лежащие на острых кромках А, В, С считаются свободными [15]. Между присоединенными вихрями располагались контрольные точки. Введем обозначения:
– точка расположения k-го присоединенного вихря;
– p-я контрольная точка. Тогда скорость в точке
вдоль единичного направления n, индуцированная вихрем
, расположенным в точке
, определится из выражения
.
Полагалось, что на свободной поверхности тока, стекающей с кромки С, интенсивность вихрей постоянна и равна . Расстояние между свободными вихрями есть величина постоянная и равная h. Первое приближение для свободной линии тока выбиралось следующим образом: первые 3 вихря располагались параллельно
и ниже точки срыва, начиная с острой кромки, остальные – параллельно -
и левее этой точки.
Для свободных поверхностей тока, сходящих с кромок А и В, начального приближения не задавалось. Циркуляции на этих линиях тока заранее неизвестны, но они равны по абсолютной величине и противоположны по знаку, поскольку вращение частиц жидкости, сорвавшихся с этих кромок, осуществляется в противоположных направлениях.
Обозначим N – количество присоединенных вихрей, включая два свободных на кромках А и В; номер вихря в точке А обозначим N1; номер вихря в точке В − N; – точка расположения свободного вихря;
– интенсивность вихрей на свободных поверхностях тока, срывающихся с острых кромок А и В;
– эти же интенсивности найденные на предыдущей итерации;
– количество свободных вихрей, расположенных на свободной поверхности тока, срывающейся с острой кромки C;
– тоже, но с кромки А;
– количество свободных вихрей на всех линиях тока, на первой итерации
.
Система линейных алгебраических уравнений для определения неизвестных интенсивностей присоединенных вихрей имеет вид:
(1)
где , на первой итерации
.
После определяются неизвестные циркуляции , где
присоединенных вихрей и величина интенсивности свободных вихрей
. Скорость в любой точке
области вдоль любого заданного направления вычисляется по формуле:
На первой итерации, после определения неизвестных циркуляций вихрей, строятся все три поверхности тока, стекающие с острых кромок. После того как они удалятся от всасывающего проема трубы на расстояние больше 10 калибров, их построение прекращается. При построении поверхностей тока шаг во времени выбирается достаточно малым. Построение поверхности тока начинается с острой кромки.
После определения второго приближения для свободных поверхностей тока необходимо снова решить систему уравнений (1) и определить циркуляции присоединенных вихрей. Затем строится третье приближение свободных поверхностей тока и т.д.
Если расстояние от точки до вихря меньше радиуса дискретности h/2, то скорость, вызываемая этим вихрем, определяется из следующего выражения:
.
С использованием разработанной компьютерной программы решена задача о моделировании отрыва потока на входе в круглый всасывающий канал с кольцевым механическим экраном, снабженным выступом с целью определения оптимальных геометрических размеров и удаления этого экрана, способствующих наибольшему аэродинамическому сопротивлению входа среды в этот канал.
Из расчетов линий тока (рис. 2) видно, что снабжение кольцевого экрана выступом значительно изменяет структуру течения. Отрывная зона на входе в трубу расширяется, образуется дополнительная область отрыва потока под выступом.
Рис. 2. Линии тока при обтекании тонкого экрана с центральным круглым отверстием при ;
;
;
; а) без выступа; б) с выступом
Предложена следующая формула связи между к.м.с. всасывающего круглого патрубка и коэффициентом сжатия струи
:
, (2)
расчеты по которой дают приемлемые для практики результаты (рис. 3).
Как видно из данных, представленных на рис. 3, снабжение кольцевого экрана выступом заметно повышает аэродинамическое сопротивление всасывающего канала. Максимальное значение к.м.с. достигается при коэффициенте сжатия струи
, что следует из анализа функции (2) с помощью дифференциального исчисления.
Вычислительный эксперимент позволил найти оптимальные геометрические размеры экрана и его удаление
, способствующие наибольшему росту к.м.с. Для круглых всасывающих отверстий, выступающих из плоской стенки на расстояние полурадиуса, предлагается устанавливать на расстоянии 0.1 – 0.2 радиуса диск с центральным отверстием такого же диаметра, что и всасывающий канал, и шириной кольца равной 0.7 – 1 радиуса, что позволяет увеличить к.м.с. на 20% относительно случая отсутствия этого экрана.
Рис. 3. Сравнение расчетных и экспериментальных величин к.м.с. при удалении экрана от всасывающего проема (пунктирная линия – расчет по формуле, предложенной И.Е. Идельчик; сплошная линия – расчет по формуле (2); кружочки – данные эксперимента, полученные нами на экспериментальной установке):
а) ;
; б)
;
: в)
;
; г)
;
Сравнение результатов расчетных значений к.м.с. (рис. 3) и натурного эксперимента, выполненного на лабораторной установке (рис. 4) показывают достоверность расчетов по формуле (2).
Рис. 4. Схема экспериментальной установки для определения сопротивления при входе воздуха
в экранированное круглое отверстие: 1 - экран с центральным отверстием и выступом; 2 – щит;
3 – направляющая трехгранная призма; 4 – труба; 5 – микроманометр с наклонной трубкой;
6 – пневмометрическая трубка Пито-Прандтля; 7 – стальные стержни-шпильки; 8 – гайки для фиксации экрана
Экспериментальное значение к.м.с. при входе воздуха в экранированный воздуховод определяется в соответствии с общепринятым соотношением , представляющим собой отношение разности полных давлений к динамическому (скоростному) давлению. Средняя скорость u воздуха в трубе определялась с помощью измерения скоростного давления по оси трубы и внесения поправочного коэффициента на неравномерность поля скоростей, величина которого была определена измерениями скорости в четырех точках равновеликих колец. Средняя скорость воздуха в трубе составляла 7 – 9 м/с.
Проведено также экспериментальное исследование отрывных течений на входе в щелевидные всасывающие каналы Экспериментальная установка (рис. 5) – модель аспирационного укрытия представляет собой короб с подведенной к нему вытяжкой. С торцевой стороны короба находится неплотность, высота которой фиксирована. Для измерения давлений внутри укрытия выведены два штуцера: один для замера статического давления, второй, находящийся в трубе вытяжного устройства, – для замера статического и динамического давлений. Стенки укрытия выполнены из прозрачного материала.
Исследовалось влияние механических экранов в виде козырька, пластины и прямого двугранного уголка (рис. 6) на коэффициент местного сопротивления неплотности аспирационного укрытия.
Рис. 5. Модель аспирационного укрытия
Рис. 6. Схема расположения козырька и экрана в виде уголка
Были определены параметры механического экрана в виде прямого двугранного уголка (рис. 6), позволяющие добиться наибольшего повышения кмс на входе в неплотность.
При помощи шлирен-метода были получены снимки эффекта срыва потока воздуха с края козырька (рис. 7).
Рис. 7. Срыв потока воздуха с края козырька, установленного над щелевой неплотностью укрытия, вид сбоку
Для оценки эффективности того или иного приспособления, влияющего на аэродинамику потока на входе в укрытие использовался коэффициент местного сопротивления. Коэффициент сопротивления , приведенный к скорости потока
в сечении неплотности аспирационного укрытия (
) рассчитывался при помощи формулы:
, скорость в том же сечении
, где
– атмосферное давление, Па;
– среднее значение полного статического давления в укрытии, Па;
– динамический напор потока воздуха в сечении трубы, Па,
– площадь потока воздуха в сечении трубы, м2;
– площадь потока воздуха в сечении неплотности, м2,
– скорость воздушного потока в сечении трубы, м/с,
– плотность окружающего воздуха,
, при температуре tизм, оС [16, 17].
Рис. 8. Изменение величины от
при a/B = 1
Рис. 9. Решетка козырьков
В результате экспериментов было выяснено, что при снабжении щелевой неплотности с горизонтальным козырьком длиной 0,5 – 0,7 калибра [11, 13] экраном в виде уголка высотой в 1 калибр и шириной 0,5 калибра на расстоянии до козырька в диапазоне 0,5 – 1 калибра к.м.с. увеличивается на величину более 16,6 %. Максимальная величина отклонения к.м.с. при оборудовании неплотности горизонтальным козырьком и экраном в виде уголка, при этом расход воздуха через неплотность снижается на
(рис. 8). Кривые 1 - 6 соответствуют величинам длины горизонтальной части уголка b/B = 1/6; 1/3; 1/2; 2/3; 5/6; 1.
Обнаружено также, что использование решетки козырьков из металлических пластин размером общей длиной 0,7 калибра, наклоненных под углом
, расстояние между которыми 0,117 калибра (рис. 9) приводит к увеличению к.м.с. на
и снижению расхода воздуха на
(к.м.с. входа в неплотность
,
).
Выводы. На основании стационарных дискретных вихревых колец разработана математическая модель, вычислительный алгоритм и компьютерная программа расчета отрыва потока на входе во всасывающий круглый канал, снабженный механическим кольцевым экраном с выступом. Вычислительный и натурный эксперимент позволили установить размеры и расположение экранов, способствующих наибольшему сопротивлению входа во всасывающий канал за счет эффекта отрыва струи. Предложена эмпирическая формула связи коэффициента местного сопротивления и коэффициента сжатия струи.
На разработанном лабораторном стенде аспирационного укрытия с неплотностью прямоугольной формы проведены экспериментальные исследования влияния горизонтальных, вертикальных, наклонных и двухгранных непроницаемых механических экранов на коэффициент местного сопротивления входа в неплотность. Показано, что оборудование щелевой неплотности горизонтальным козырьком и экраном в виде уголка при определенных параметрах расположения позволяет снизить расход воздуха, поступающего в аспирационное укрытие, на 17 %.
1. Chzhen P. Otryvnye techeniya. M.: Mir. T.1. 1972. 300s.
2. Chzhen P. Otryvnye techeniya. M.: Mir. T.2. 1973. 280s.
3. Chzhen P. Otryvnye techeniya. M.: Mir. T.3. 1973. 334s.
4. Posohin V.N. Aerodinamika ventilyacii. M: AVOK-PRESS. 2008. 209s.
5. Logachev K.I., Posohin V.N. Raschet techeniya vblizi kruglogo vsasyvayuschego patrubka // Izv. vuzov. Aviacionnaya tehnika. 2004. №1. S. 29-32.
6. Posohin V.N., Katkov M.V. Eksperimental'noe izuchenie vihrevyh zon v potokah vblizi vsasyvayuschih schelevyh otverstiy // Izvestiya vuzov. Aviacionnaya tehnika. 2001. №1. S. 61-63.
7. Averkova O.A., Logachev I.N., Logachev K.I.. Otryvnye techeniya v spektrah vytyazhnyh kanalov. Moskva-Izhevsk: IKI. 2012. 288s.
8. Averkova O.A., Logachev I.N., Logachev K.I. Modelirovanie potencial'nyh techeniy s neizvestnymi granicami na osnove stacionarnyh diskretnyh vihrey // Vychislitel'nye metody i programmirovanie. 2011. T.12. №2. C. 213-219.
9. Averkova O.A., Logachev A.K., Logachev I.N., Logachev K.I. Zakonomernosti otryvnogo techeniya pri vhode v vystupayuschiy kanal s ekranami // Uchenye zapiski CAGI. 2013. T.44 №2. S. 33-49.
10. Averkova O.A. Eksperimental'noe issledovanie otryvnyh techeniy na vhode vo vsasyvayuschie otverstiya // Vestnik BGTU im. V.G. Shuhova. 2012. № 1. S. 158-160.
11. Logachev I.N., Logachev K.I., Zorya V.Yu., Averkova O.A. Modelirovanie otryvnyh techeniy vblizi vsasyvayuschey scheli // Vychislitel'nye metody i programmirovanie. 2010. T.11. №1. 43-52.
12. Logachev I.N., Logachev K.I., Averkova O.A. Matematicheskoe modelirovanie otryvnyh techeniy pri vhode v ekranirovannyy ploskiy kanal // Vychislitel'nye metody i programmirovanie. 2010. T.11. №1. 68-77.
13. Averkova O.A., Zorya V.Yu., Logachev K.I., Logachev I.N. Chislennoe modelirovanie vozdushnyh techeniy na vhode v schelevye neplotnosti aspiracionnyh ukrytiy// Novye ogneupory. 2010. №5. S. 31-36.
14. Averkova O.A. Modelirovanie otryvnyh techeniy na vhode vo vsasyvayuschie kanaly s ispol'zovaniem metoda diskretnyh vihrey v nestacionarnoy postanovke // Vestnik BGTU im. V.G. Shuhova. 2012. №2. S. 133-135
15. Lifanov I.K. Metod singulyarnyh integral'nyh uravneniy i chislennyy eksperiment. M.: Yanus. 1995. 520s.
16. Idel'chik I.E. Spravochnik po gidravlicheskim soprotivleniyam. III izdanie. M.: «Mashinostroenie». 1992. 672s.
17. Veleckiy R.K., Grigina N.N. Izmerenie parametrov pylegazovyh potokov v chernoy metallurgii. M.: «Metallurgiya». 1979 g. 77s.
