employee
employee
employee
VAC 05.17.00 Химическая технология
VAC 05.23.00 Строительство и архитектура
GRNTI 20.53 Технические средства обеспечения информационных процессов
The study of vortex flow character and parameters is important problem to calculation of values of different technological processes like grinding, classification of solid phase, dust and ash cleaning and others. The mathematical and simulation modeling of processes, occurring in vortex cameras, may be use for calculating of energy carriers aerodynamic parameters in these cameras. Mathematical modeling of processes, occurring in vortex cameras, based on determination of Navier-Stokes equation system for particular initial and boundary conditions. The setting of these conditions is depending of decided problem, dimensions of vortex camera and the energy carrier injection to the camera process specifications. In article there was also given the results of simulation modeling of velocity and pressure profiles in vortex camera. The offered modeling methods make possible the piloting of aerodynamic parameters in vortex cameras.
vortex flow, mathematical modeling, simulation modeling
Введение. Микроструктура закрученного потока в вихревых камерах представляет интерес для понимания физических механизмов проходящих в них процессов, таких, как, например, измельчение дисперсных материалов в вихре-акустических диспергаторах [1–3], классификация твердой фазы [4, 5], пыле- и золоулавливание [6], тепломассообмен [7] и других [8, 9]. Поток энергоносителя должен обладать необходимыми скоростными характеристиками, несущей способностью для организации технологических процессов. Определение аэродинамических параметров энергоносителя (воздуха) в вихревых камерах возможно с использованием математического и имитационного моделирования протекающих в указанных камерах процессов.
В данной статье будут рассмотрены короткие вихревые камеры цилиндрической формы с вертикальной осью, в верхней крышке которых встроены четыре цилиндрических резонатора. Резонаторы создают дополнительные возмущения потока в камере. Подвод энергоносителя осуществляется тангенциально, выход энергоносителя проходит через центральное отверстие в верхней крышке вихревой камеры [10].
При сильной закрутке газовых потоков в коротких вихревых камерах радиальная скорость у боковой стенки камеры значительно превышает среднерасходные значения, а в окрестности выходного патрубка образуется зона обратных течений. Эту область хорошо наблюдать при подаче в вихревую камеру воды. Зона обратных течений располагается не в центре, а смещена относительно центральной оси камеры. Величина смещения при этом зависит от параметров энергоносителя на входе в камеру.
Особенностью закрученных вихревых потоков является возможность формирования в камере областей с активным и консервативным характером воздействия центробежных массовых сил на структуру потока [11–13].
Основная часть. Если в вихревой камере организован один тангенциальный вход энергоносителя, то вихревой поток в ней неоднороден по пространству камеры, а также не является осесимметричным.
В сильно закрученных потоках имеются значительные радиальные градиенты давления, которые в свою очередь приводят к образованию центральной вихревой зоны. Наличие такой зоны с интенсивной завихренностью способствует выполнению ряда требований, предъявляемых к вихревым камерам.
Если скорость потока обозначить как u, а угловую скорость – ω, то, учитывая, что , система уравнений Навье – Стокса для несжимаемой жидкости в цилиндрической системе координат z, r, φ для коротких плоских вихревых камер записывается в виде:
; (1)
; (2)
; (3)
,
,
, (4)
где σrr, σrφ, σφφ – компоненты тензора напряжений, μ – коэффициент вязкости.
Приведенной выше системой определяются поля радиальной ur и окружной uφ составляющих вектора скорости и давление p. В случае вертикального расположения оси короткой вихревой камеры можно принять допущение, что радиальная и окружная составляющие вектора скорости не зависят от координаты z. Коэффициент вязкости считается не зависящим от координат, но зависящим от закрутки потока.
Приведенную систему уравнений следует решать при дополнении начальных и граничных условий. Постановка этих условий определяется конкретной решаемой задачей, геометрическими параметрами, которыми обладает вихревая камера, а также технологическими условиями подачи энергоносителя в вихревую камеру.
На срезе и стенках сопла будут выполняться условия непротекания, поэтому . Подставив это соотношение в систему, можно получить
. (5)
Движение газа из ресивера через сужающее сопло описывается выражением [2]
, (6)
где u – скорость энергоносителя выходящего из сопла в камеру; p0 – начальное давление газа в ресивере; p – давление в рассматриваемом сечении; k – показатель адиабаты.
Массовый расход газа при изоэнтропическом истечении через сопло [14]
, (7)
где ρ0 – начальная плотность газа в ресивере; ρ – плотность газа в рассматриваемом сечении.
Значение объемного расхода воздуха Q зависит от величины давления, м3/мин
, (8)
где dс – диаметр сопла.
Давление pk, при котором произведение ρV имеет максимальное значение, называется критическим. Оно может быть определено из соотношения [11]
, (9)
Для воздуха, при k = 1,41, критическое отношение давлений pk/p0 = 0,528. Расход газа имеет максимальное значение при критическом отношении давлений.
Струя воздуха, попадая в вихревую камеру, взаимодействует с внутренним потоком энергоносителя [15]. Деформация сечения струи обусловлена разностью давлений в различных точках сечения. Со стороны набегающего потока давление всегда оказывается больше, чем с противоположной стороны. Глубина проникновения круглой струи в поток
, (10)
где d0 – диаметр сопла; a – скорость звука; u1 – скорость потока в окрестности входящей струи; u2 – скорость струи на входе в поток; ρ1 – плотность газа в потоке; ρ2 – плотность газа в струе.
Для прямоугольной струи
, (11)
где – эквивалентный диаметр струи; Aс – площадь сопла.
Рис. 1. Схема взаимодействия струи внутренним
потоком в вихревой камере
Угол отклонения струи от плоскости симметрии φс связан с углом косого среза сопла ψ следующим соотношением
, (12)
При ψ = 45° расчетное значение угла φс отклонения струи составляет 21°. При ψ = 45° расчетное значение угла φс составляет 8°, замеры, проведенные по результатам износа, оказались равными 5°. Такое отличие объясняется наличием дисперсной фазы в набегающем потоке.
Путь, который входная струя энергоносителя проходит до встречи со стенкой вихревой камеры зависит от радиуса камеры и от угла и вычисляется по формуле
, (13)
Для оценки влияния трения газа о корпус помольной камеры можно ввести в рассмотрение некоторую «эквивалентную трубу». Для поддержания равномерного движения вязкого газа в трубе необходимо к сечениям трубы длины l приложить движущий перепад давлений [14]
, (14)
где d – диаметр трубы; uср – средняя скорость в трубе; λ – коэффициент гидравлического сопротивления.
Для моделирования вихревых течений уравнения Навье – Стокса осредняются по Рейнольдсу. В проведенных расчетах используется осредненное по малому масштабу времени влияние турбулентности на параметры вихревого потока. При этом крупномасштабные временные изменения осредненных по малому масштабу времени составляющих газодинамических параметров потока (давления, скорости, температуры) учитываются введением соответствующих производных по времени. Для замыкания системы уравнений используются уравнения переноса кинетической энергии турбулентности и ее диссипации в рамках k–ε модели турбулентности.
Проведено численное исследование сформулированной задачи в сечении, проходящем через срез сопла и центр вихревой камеры. Результаты численных расчетов для вихревой камеры с диаметром равным 0,2 м и высотой 0,03 м представлены на рис. 2.
Численные расчеты скорости и давления показали, что поток энергоносителя в вихревой камере с одним тангенциальным патрубком подвода не является осесимметричным. Значения скорости зависят от диаметра выходного отверстия, от скорости выходящего потока из входного сопла в вихревую камеру.
Имитационное моделирование позволяет получать распределение аэродинамических параметров энергоносителя в вихревой камере реального размера. Результаты имитационного моделирования распределения скорости представлены на рис.3.
На рис. 3 видно распределение скорости энергоносителя в вихревой камере, а также поведение энергоносителя в области резонаторов. Полученные результаты расчетов отражают структуру в горизонтальной плоскости, в объеме вихревой камеры и на выходе через центральный патрубок.
Результаты расчетов распределения давления в вихревой камере представлены на рис. 4.
На рисунках видно, что расположение зоны с обратным течением зависит от входных параметров энергоносителя и от геометрических размеров подводящих патрубков.
Рис. 2. Расчет параметров энергоносителя в вихревой камере:
а – скорости; б – давления
Рис. 3. Распределение скоростей энергоносителя в вихревой камере:
а – плоский случай; б – объемный случай
Рис. 4. Распределения давления энергоносителя в вихревой камере при диаметре патрубка подачи
энергоносителя в камеру:
а – d = 0,015 м; б – d = 0,010 м
Выводы. Предложенный метод математического и имитационного моделирования, проведенные исследования и результаты расчетов позволяют осуществлять управление аэродинамическими параметрами в вихревых камерах для организации различных технологических процессов.
1. Ural'skiy A.V., Kolesnikov A.V., Perelygin D.N. Sinica E.V. Tehnologicheskie moduli dlya kompleksnogo izmel'cheniya materialov // Vestnik Belgorodskogo gosudarstvennogo tehnologicheskogo universiteta im. V.G. Shuhova. 2007. №3. S. 80-86.
2. Perelygin D.N., Sevost'yanov V.S., Boychuk I.P., Kuznecova I.A. Issledovanie udarnogo kontakta abrazivnyh chastic s futerovkoy vihre-akusticheskogo dispergatora // Vestnik Belgorodskogo gosudarstvennogo tehnologicheskogo universiteta im. V.G. Shuhova. 2015. №2. 122-125.
3. GorlovA.S., Gubarev A.V. Matematicheskoe opisanie processov udarnogo razrusheniya chastic tverdoy fazy v kamere izmel'cheniya vihreakusticheskogo dispergatora // Fundamental'nye i prikladnye problemy tehniki i tehnologii. 2015. № 5. S. 40-45.
4. Varaksin A.Yu. Turbulentnye techeniya gaza s tverdymi chasticami. M.: FIZMATLIT, 2003. 192 s.
5. Sevost'yanov V.S., Perelygin D.N. Ural'skiy V.I. Gorlov A.S., Glagolev E.S., Babaevskiy A.N. Razrabotka i issledovaniya energosberegayuschego pomol'nogo oborudovaniya dlya vysokodispersnogo izmel'cheniya materialov // Vestnik Belgorodskogo gosudarstvennogo tehnologicheskogo universiteta im. V.G. Shuhova. 2014. №3. 76-80.
6. Shvydkiy V.S., Ladygichev M.G. Ochistka gazov. M.: «Teplotehnik», 2005. 640 s.
7. Zhukauskas A.A. Konvektivnyy perenos v teploobmennikah. M.: Nauka, 1982. 472 s.
8. Il'ina T.N., Boychuk I.P., Emel'yanov D.A. O vzaimodeystvii voloknistyh materialov v vozdushnom potoke // Vestnik Belgorodskogo gosudarstvennogo tehnologicheskogo universiteta im. V.G. Shuhova. 2016. №6. 116-121.
9. Sevost'yanov M.V., Il'ina T.N., Kuznecova I.A., Shkarpetkin E.A., Shinkarev L.I. Oborudovanie dlya kompaktirovaniya tehnogennyh materialov // Vestnik Belorussko-Rossiyskogo universiteta. 2016. № 2. S. 92-101.
10. Pat. 2250138 Rossiyskaya Federaciya, MPK7 B 02 C 19/06. Vihre-akusticheskiy dispergator / A.M. Gridchin, V.S. Sevost'yanov, V.S. Lesovik, A.S. Gorlov, D.N. Perelygin, B.Z. Fedorenko; zayavitel' i patentoobladatel' BGTU im. V.G. Shuhova. № 2003123664/03, zayavl. 25.07. 03; opubl. 20.04. 05 Byul. № 11. 8 s.
11. Gupta A., Lilli D., Sayred N. Zakruchennye potoki. M.: Mir, 1987. 345 s.
12. Gol'dshtik M.A. Vihrevye potoki. Novosibirsk: Nauka, 1981. 367 s.
13. Smul'skiy I.I. Aerodinamika i processy v vihrevyh kamerah. Novosibirsk: VO «Nauka». Sibirskaya izdatel'skaya Firma, 1992 g. 301 s.
14. Al'tshul' A.D. Zhivotovskiy L.S., Ivanov L.P. Gidravlika i aerodinamika. M.: Stroyizdat, 1987. 414 s.
15. Abramovich G.N., Girshkovich T.A., Krasheninnikov S.Yu., Sekundov A.N., Smirnova I.P. Teoriya turbulentnyh struy / Pod red. G.N. Abramovicha. M.: Nauka. Glavnaya redakciya fiziko-matematichkskoy literatury, 1984. 717s.
