CC BY
21
Эволюция акустического излучения ансамбля вихревых колец
в воздухе
Д. Ю. Черкасовa, Ф.В. Шугаевb
Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра квантовой статистики и теории поля.
Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2. E-mail: a [email protected], b [email protected] Статья поступила 20.03.2017, подписана в печать 19.06.2017.
Исследована эволюция акустического излучения от ансамбля вихревых колец в воздухе на основе нестационарных уравнений Навье-Стокса. Использованы разложения искомых функций в ряд по степеням начальной завихренности, которая считается малой величиной. Система уравнений Навье-Стокса сводится к параболической системе с постоянными коэффициентами при старших производных. Задача ставится следующим образом. Завихренность определена внутри тороида при t = 0. Остальные параметры газа полагаются постоянными во всем пространстве в начальный момент времени. Решение выражается через кратные интегралы, которые рассчитываются с помощью сеток Коробова. Исследованы колебания плотности. Результаты показывают, что спектр частот зависит от времени, а именно при малых временах наблюдаются высокочастотные колебания, затем появляются низкочастотные колебания. В то же время амплитуда высокочастотных колебаний уменьшается по сравнению с низкочастотными. Таким образом, происходит переход энергии от высокочастотного спектра к низкочастотному. Представленные результаты могут быть полезны для моделирования затухающей сеточной турбулентности.
Ключевые слова: ансамбль вихревых колец, акустическое излучение, уравнения Навье-Стокса, переход энергии от высокочастотных колебаний к низкочастотным.
УДК: 534-13. PACS: 47.32.-y.
Введение
Вихревые структуры (вихревые кольца и цилиндрические вихри) играют важную роль в излучении звука в газовых потоках. Излучение звука вихревыми структурами было детально рассмотрено в [1, 2]. Исследование акустического излучения от системы двух вихревых колец выполнено в [3, 4]. Влияние вихревого кольца на свойства турбулентного потока подтверждено многочисленными экспериментами [5, 6] и расчетами [7, 8]. Акустическое излучение изолированного вихревого кольца в воздухе проанализировано в [9]. Мы исследовали эволюцию акустического излучения, испускаемого ансамблем вихревых колец в воздухе. Анализ основан на нестационарных уравнениях Навье-Стокса. Используются разложения искомых функций в ряд по степеням малого параметра. Уравнения Навье-Стокса сводятся к параболической системе с постоянными коэффициентами. В результате мы получаем решение в виде степенного ряда с кратными интегралами при постоянных коэффициентах. Первый член ряда определяет свойства акустического излучения при малой начальной завихренности, которая служит малым параметром. Далее анализируется эволюция излучения с течением времени.
1. Основные уравнения
Используется разложение Гельмгольца поля скоростей на потенциальную и соленоидальную части.
Уравнения Навье-Стокса в безразмерной форме могут быть представлены в виде [9] дО-
-О = ехр[да + 0.75А] АО +
Ш
. 3 г . / дvk дvm\ д2Н
+ 4в,. ехр[. + °.75^ — + —) дддХт
da dvi - + "m dx— stti+/1i,
dxj
dw dw
Ht = -vi Щ + s' (1)
ds = — Aw + 4 exp[w + 0.75h] As -dt y 3
- ^1 eh + 0.5sexp[w + 0.75h]^ Ah +
nnr-n dvi d2h ds + 1.5 exp[w + °.75h] ^f - vj — + /2,
д Xj dxidxj
dh y
dt = Pr exp[w + 0.75h] Ah - (y - 1) s - Vj — -2
- 3y(y - 1) exp[w - 0.25h] s2 + /3, " = rot v, s = div v, w = - ln p, h = ln T,
V = м/p, Pr = ^p, A = д2
j dxj dh
Л° д%1 д%1
I = 1,2,3, , = 1,2,3, . = 1,2,3, т = 1,2,3.
Здесь VI — компонента скорости, — антисимметричный тензор, р, Т, v — безразмерные вели-
чины плотности, температуры, скорости (деленные на р0 = 1.225 кг/м3, Т0 = 300 К, с0 соответственно), р, V, Л — коэффициенты динамической вязкости, кинематической вязкости и теплопроводности, с — низкочастотная скорость звука, с0 — скорость звука в воздухе, 7 — показатель адиабаты, Рг — число Прандтля, ср — теплоемкость при постоянном давлении. Значения времени и координат отнесены к соответствующим характеристическим величинам: характеристической длине 10 = ^/с0 и характеристическому времени ¿0 = ^/с0.
При выводе системы (1) учтено уравнение Менделеева-Клапейрона, а также считается, что зависимости вязкости р и теплопроводности Л от температуры выражаются степенным законом:
М/М0 = Т0 75, Л/Л0 = Т0 75.
Функции /ц, /2, [3 — нелинейные члены по отношению к первым производным по координатам. Индекс «0» относится к начальному состоянию, которое должно быть однородным.
1.1. Начальные данные задачи
Как известно, сеточная турбулентность в максимальной степени близка к однородной изотропной турбулентности. Сеточная турбулентность возникает при прохождении газового потока через систему отверстий в плоской пластине. В результате в потоке образуются вихревые кольца. В наших расчетах ансамбль вихревых колец состоит из двух слоев по четыре кольца в каждом и в некоторой мере моделирует сеточную турбулентность.
имеет вид [10]
x' = (r00u1 sin «2 + Rc) cos «3,
= (r00u' sin «2 + Rc) sin «3, (2)
x3 = r00«' cos «2, 0 < «1 <1, 0 < «2, «3 < 2n, где r00, Rc — начальные размеры вихревого кольца, r00 — радиус сечения, Rc — радиус кольца. Считается, что безразмерная начальная завихренность мала: w0 с 1, w0 = const. Начальные условия имеют вид
О1 = -w0 sin «3, О2 = w0 cos «3, О3 = 0, (3)
внутри начального тороида и
w(x,0) = s(x, 0) = h(x, 0) = 0
во всем остальном пространстве.
(4)
1.2. Решение задачи
Уравнения (1) являются нелинейной параболической системой. Мы ищем решение параболической системы в виде разложения в степенной ряд. Учитывая уравнение (1) и начальные условия, получим
О, (х, 0 = £ 0£1)(х, ^ + £2 0£2)(х, ъ + £3 0£3)(х, ^ + ..., О (х, $ = £2 О32) (х, *) + £3 О33) (х, £4 О34) (х, $ + . . . , v¡(x, *) = £ 0,(1)(х, *) + £2 v(2)(x, *) + £3у(3)(х, *) + . . . ,
£ = ш0, а = 1,2, I = 1,2,3. (5)
Подставляя (5) в уравнения (1), имеем ш(х, *) = £2 Ш(1)(х, *) + £3 Ш(2)(х, *) + £4 Ш(3)(х, *) + . . . , 5(х, *) = £2 5(1)(х, *) + £3 5(2)(х, *) + £4 5(3)(х, *) + . . . ,
Й(х, *) = £2 х, *) + £3 х, *) + £4 х, *) + ... .
(6)
Для функций наинизшего порядка уравнения таковы:
дО(1)
dt dw(1) dt
ds(1)
= ДО
= s(1),
(1)
= 1 дш(1) + 4 Д5(1) - 1 Дй.(1) +
dt y 3 y 3
(7)
dh(1)
(1)
_ = J_ Дй(1) - (Y - 1) s(1)+ ^
^ = '-i- j ^41) = 2(y - 1) j D^1*
Dj = dx/
4 (1)
+
dx„
Из уравнения (1) следует:
Рис. 1. Ансамбль из восьми вихревых колец
В начальный момент завихренность имеет ненулевое значение только внутри восьми газовых торои-дов (рис. 1). Уравнение отдельного вихревого кольца
v(1)(x, t) = -
0.03125
n5/2t3/2
di
dr'
sin 0' dd' x
r3
2n
d4>' 0(1)(i,0) x nexp -|x + x' - i|
0
0.25
x
n = {sin в' cos ф', sin в' cos ф', cos в'} , x' = {r' sin в' cos ф', r' sin в' cos ф', r' cos в'} . (8)
Система (7) состоит из трех однородных параболических уравнений относительно Ц(1) и неоднородной параболической подсистемы. Все уравнения имеют постоянные коэффициенты при старших производных. Решение системы может быть получено с помощью преобразования Фурье. Первое уравнение из (8) дает
2п
a(1)(x, t) =
0.125
(nt)3/2
R3
(-0t5|x - о
Ц(1)(|,0)ехр í-^|x -||2 )d|.
(9)
Уравнения (9) позволяют определить члены ^41). Преобразование Фурье однородной параболической подсистемы (8) дает
йй(1)
dt
= 7(1),
ds(1) k2~ (1) 4,2~(1) k2 ~(1)
—— =--w(1) - ттk2s +--/г(1),
dt y 3
Y
(10)
^=-k2Pr*(1) - (y -1)
/3++£) /2+k
w(1)(x, t) =
1
(2n)3/2
dr
dk
d| exp{ik(x - I) x
R3
R3
x {a^^I-t) + a13^(I, r)} . (13)
Здесь a12, a13 — элементы фундаментальной матрицы решений системы (10). Соответствующие формулы даны в приложении.
Введем новую переменную:
X = I - x, (14)
X = {R1 sinв1 cos^i, Ri sinв1 sinR1 cosв[} .
Замена переменных дает
w(1)= f-
©
1/2
dr
k dk
sin(kR1) R1 dR1
sin в1 dв1 x
d^1 {a12^31)(x + X, r) + a13^41)(x + X, r)} . (15)
Аналогичные выражения можно написать для функций 5(1)^, ¿), Н(1)(х,
Отклонение плотности от ее начального значения может быть записано в виде
(16)
Здесь индекс й обозначает размерную величину.
Функция ш(1) не зависит от ш0, как и частота колебаний плотности. Функции w(")(x, ¿), s(")(x, £), h(")(x, £), п >1 могут быть получены аналогично:
(Pd - Р0)/Р0 ~ -w « w(1).
(!)
1/2
dr
k dk
R1 dR1 sin в1 dв1 d<p1 x
x sin kR1 c,(n-1), (17)
c^n) = a¡1(k, t-r) Q2 (x+X, r) +
+ a¡2(k, t-r) Q3n)(x+X, r) + a¡3(k, t-r)Q4n)(x+X, r),
(18)
w(n)(x, t)= y{n), s(n)(x, t)= y2n), h(n)(x, t)= y3n),
(19)
Тильда означает преобразование Фурье с волновым числом к. Характеристическим уравнением для системы (10) будет
(& +■)/ 4=
(11)
При 0 < к < к*, к* « 1 для воздуха, корни уравнения (11) таковы:
/1 = (п(к:), /2,3 = (г(к) ± I ш (к), (Ь ( < 0. (12)
При к > к* все корни быстро затухают, поэтому этот случай не рассматривается. Дисперсионная кривая шг(к) при 0< к < к* имеет две ветви. Мы рассматриваем ветвь, относящуюся только к меньшим величинам коэффициентов затухания ( 1, ( 2 (0 < к < к1 < к*). Для функции ш(1) имеем
q22) =
q32) = -2
dw(1) dx¡
dv,(1) dvj2) _dv([) d2h(1) (1) ds(1)
dx¡ dx¡
+ 1.5
dx¡ 3x¡ dx¡
v
(1)
dxn
'dw(1 дП;(1) + 15 dh(1) dÜ(kiy
dx¡ dxn
dx¡ dxn
(20)
q42) = -v
(1)
dh(1) dxn
+
+ y(y -1)
l dx¡ dx¡ I l dx¡ dx¡ I
dvf
Значения Q2"), \ 04п) при п >2 имеют громоздкий вид и поэтому не приводятся.
Таким образом, уравнения Навье-Стокса были сведены к параболической системе с постоянными коэффициентами при производных. В результате мы получаем степенной ряд по ш0. Коэффициенты ряда являются известными функциями x, I (кратные интегралы). Первые члены ряда могут быть использованы для анализа спектра частот колебаний плотности для случая малой завихренности.
Уравнение (15) использовалось для изучения ос-цилляций плотности. Кратный интеграл был рассчитан с помощью сеток Коробова [12]. Поле плотности было рассмотрено вблизи ансамбля вихревых колец. Параметры каждого кольца таковы: Яс = 0.15 см, г00 = 0.03 см, отношение радиуса поперечного сечения кольца к радиусу кольца равно 0.2, расстояние между центрами соседних колец равно 0.72 см. При вычислениях использовалось значение ш0 = 0.0001. Значения плотности рассматривались в точке х1 = х2 = х3 = 0.
х
2. Результаты
Расчеты выполнены для ансамбля из восьми вихревых колец. Схема их расположения показана на рис. 1. Исследовались осцилляции плотности в воздухе. Значения параметров таковы: 7 = 1.4, Pr = 0.71, число Рейнольдса Ке = 171. В качестве характеристической длины в выражении для числа Ке взят начальный диаметр кольца.
Численные результаты представлены на рис. 2-8.
Колебания плотности на ранней стадии представлены на рис. 2. Видно, что возникают высокочастотные колебания, модулированные низкочастотным сигналом. Высокая частота равна 117 кГц, низкая низкая частота равна 3.5 кГц (рис. 3). Как можно видеть, имеются только два значения частоты,
т, с
к ч
с
§
<
50 ООО 100 ООО 150 000
Частота, Гц
Рис. 3. Спектр частот колебаний плотности вблизи ансамбля вихрей, I = 0.19-1.2 мс (ранняя стадия)
г, с
Рис. 2. Колебания плотности вблизи ансамбля вихревых колец (восемь колец, Rc = 0.15 см, гю = 0.03 см, ранняя стадия)
Рис. 4. Колебания плотности вблизи ансамбля вихревых колец (восемь колец, Rc = 0.15 см, поздняя стадия)
20000
40000
60000 80000 Частота, Гц
100000
120000
140000
Рис. 5. Спектр частот колебаний плотности вблизи ансамбля вихревых колец, I = 3.84-4.93 мс
(поздняя стадия)
Т с
Рис. 6. Дальнейшая стадия осцилляций плотности вблизи ансамбля вихревых колец
0 20000 40000 60000 80000 100000 120000 140000 160000
Частота, Гц
Рис. 7. Спектр частот осцилляций плотности вблизи ансамбля вихревых колец, * = 5.78-6.98 мс
0(0/0(0)
1.0
0.8
0.6
0.4
0.2
0
_|_I_1_
_|_I_I_1_
1.0
1.5 Т, мс
Рис. 8. Затухание завихренности в центре поперечного сечения одиночного вихревого кольца
амплитуда низкочастотных колебаний меньше, чем амплитуда высокочастотных. Картина существенно меняется с течением времени (рис. 4, 5). Наблюдается высокая частота и помимо этого низкочастотные колебания появляются при 23 кГц и при 53 кГц (рис. 5). Их амплитуда превышает амплитуду высокочастотных колебаний. Таким образом, имеется переход энергии от высокочастотных колебаний к низкочастотным. Высокую частоту можно считать собственной частотой ансамбля. С течением времени происходит дальнейшее изменение спектра осцил-ляций (рис. 6). Как можно заметить, амплитуда
высокочастотных колебаний уменьшается еще больше, а низкочастотные колебания становятся выраженными все более ярко (рис. 7). График затухания завихренности в центре поперечного сечения кольца представлен на рис. 8. Более быстрое затухание высокочастотных осцилляций по сравнению с низкочастотными имеет простое объяснение. Как известно, величина диссипации зависит от производных скорости по координатам. Соответственно в первом случае диссипация больше, чем во втором, и, следовательно, высокочастотные осцилляции затухают быстрее низкочастотных.
Заключение
При малой начальной завихренности построено решение системы нестационарных уравнений На-вье-Стокса для анализа вихревых колец в вязком теплопроводном газе, представимое в виде интеграла. Решение использовано для исследования эволюции ансамбля. Показано, что имеет место передача энергии от высокочастотных осцилляций к низкочастотным. Одновременно происходит хаотизация течения. Сравнение осцилляций плотности в случае одиночного кольца [9] и ансамбля показывает, что во втором случае частотный спектр имеет более сложный характер, а хаотизация наступает быстрее.
Приложение
Ниже приведены значения элементов фундаментальной матрицы решений системы (10):
all = сП 21 + с2122 + с31 Zз,
al2 = с1221 - с1222 + с3223,
а13 = с13 21 - с13 22 + с33 23,
а21 = с1121 + (с2102 + с31 Шг) 22 + (с31 02 - с21 Шг) 23,
а22 = с1221 + ( —с1202 + с32Шг) 22 + (с3202 - Й2Шг) 23,
а23 = с1321 + (-с13О2 + с33Шг) 22 + (с3302 - с13Шг) 23,
а31 = с11 ^221 + (с21 gз + с31 ^4) 22 + (-с21 ^4 + с31 gз) 23, а32 = с12^221 + (-с12^3 + C32g4) 22 + (с^ + с32^3) 23, а33 = С13Я221 + (-с13^3 + с33^4) 22 + (с^ + с33^3) 23,
с11 = (7 (о2 + Шг2) - £2) /gl,
с12 = -2 (02 + 2к2/3) /Ц0, с13 = ¿2/&1,
с21 = (701 (01 - 202) + £2) /ёи
с31 = (701 (02 - 01 02 - Ш^ + &2(01 - 02)) /(Шг^1), с32 = (0? - 02 + Ш2 + 4£2(01 - 02)/3) /(Шг£0), с33 = (01 - 02)/(Шг£1),
gо = (02 - 01 )2 + Ш2, gl = 7£0, g2 = 1 + 701 (4/£2) , gз = 1 + ¿2 (02 - Шг2) + 4702/3, = 27Шг (2/3 + 02/£2),
21 =ехр[01 (* - т)],
22 = ехр[02(* - т)] соб[шг(* - т)],
23 = ехр[02(* - т)] б1п[шг(* - т)].
Список литературы
1. Lighthill M.J. // Proc. Roy. Soc. A 1952. 211. Р. 564.
2. Howe M.S. Theory of vortex sound. Cambridge University Press: Cambridge, UK, 2003.
3. Verzicco R., Ia/rati A., Riccardi G., Fatica M. // J. Sound Vibr. 1997. 200. Р. 347.
4. Inoue O. // Phys. Fluids 2002. 14. Р. 3361.
5. Maxworthy T.J. // J. Fluid Mech. 1974. 64. Р. 227.
6. Maxworthy T.J. // J. Fluid Mech. 1977. 81. Р. 465.
7. Liu C, Yan Y, Lu P. // Comput. Fluids. 2014. 102. Р. 353.
8. Yan Y, Chen C, Fu H, Liu C. // J. Turbul. 2014. 15. Р. 1.
9. Shugaev F.V., Cherkasov D.Y., Solenaya O.A. // Aerospace 2015. 2. Р. 627.
10. Morton T.S. // J. Fluid Mech. 2004. 503. Р. 247.
11. Truesdell C. // J. Ration. Mech. Anal. 1953. 2. Р. 643.
12. Коробов Н.М. Теоретико-числовые методы в приближенном анализе. М.: Физматгиз, 1963.
The evolution of acoustic radiation by an ensemble of vortex rings in air D. Yu. Cherkasova, F. V. Shugaevb
Department o/ Quantum Statistics and Filed Theory, Faculty o/ Physics, Lomonosov Moscow State
University. Moscow 119991, Russia.
E-mail: a [email protected], b [email protected].
The evolution of acoustic radiation emitted by an ensemble of vortex rings in air is studied on the basis of non-stationary Navier-Stokes equations. We use the expansions of required functions into a power series of the initial vorticity which is a small value. The Navier-Stokes equation system reduces to a parabolic system with constant coefficients for the higher derivatives. The problem is posed as follows. The vorticity is defined inside the toroid at t = 0. The other parameters of the gas are assumed to be constant throughout the space at the initial instant of time. The solution is expressed in terms of multiple integrals, which are calculated using Korobov grids. The density oscillations were investigated. The results show that the frequency spectrum depends on time; high-frequency oscillations are observed at small times and low-frequency oscillations then occur. At the same time, the amplitude of high-frequency oscillations decreases in comparison with low-frequency oscillations. Thus, a transition of energy from the high-frequency spectrum to the lowfrequency spectrum occurs. These results can be useful for modeling decaying grid turbulence.
Keywords: ensemble of vortex rings, acoustic radiation, Navier-Stokes equations, transition of energy from high-frequency oscillations to low-frequency oscillations. PACS: 47.32.-y. Received 20 March 2017.
English version: Moscow University Physics Bulletin. 2018. 72, No. 2. Pp. 173-178.
Сведения об авторах
1. Черкасов Дмитрий Юрьевич — аспирант; e-mail: [email protected].
2. Шугаев Федор Васильевич — доктор физ.-мат. наук, профессор; e-mail: [email protected].
