3. Затухание звуковых волн

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

Для вычисления диссипируемой в единицу времени энергии Е_мех воспользуемся следующими общими соображениями. Механическая энергия представляет собой не что иное, как максимальную работу, которую можно получить при переходе из данного неравновесного состояния в состояние термодинамического равновесия. Как известно из термодинамики, максимальная работа совершается, если переход происходит обратным образом (т.е. без изменения энтропии), и равна соответственно этому:

Е_мех = Е₀-Е(S),

Где Е₀есть заданное начальное значение энергии тела в состоянии равновесия с той же энтропиейS, которую тело имело вначале. Дифференцируя по времени, получаем:

Е_мех = - Е(S) = -S.

Производная от энергии по энтропии есть температура. Поэтому

– температура, которую имело бы тело, если бы оно находилось в состоянии термодинамического равновесия (с заданным значением энтропии). Обозначая эту температуру как Т₀имеем, следовательно:

Е_мех = Т₀S.

Воспользуемся для Sвыражением:

(1.21)

включающим в себя возрастание энтропии, обусловленное как теплопроводностью, так и вязкостью. Поскольку температура Т мало меняется вдоль жидкости и мало отличается от Т₀, то можно вынести ее из-под знака интеграла и писать Т вместо Т₀:

.

Эта формула представляет собой обобщение формулы

на случай сжимаемой жидкости и наличия теплопроводности.

Пусть ось х совпадает с направлением распространения звуковой волны. Тогда

Два последних члена в (1.21) дают

.

Нас, конечно, интересует среднее по времени значение величин; усреднение дает

.

(V₀– объем жидкости).

Далее, вычислим первый член в (1.21). Отклонение Т′ температуры в звуковой волне от своего равновесного значения связано со скоростью формулой

так что градиент температуры равен

.

Для среднего по времени значения от первого члена в (1.21) получаем:

.

С помощью известных термодинамических формул

можно переписать выражение в виде

.

Собирая полученные выражения, находим среднее значение диссипации энергии в виде

(1.22)

Полная же энергия звуковой волны равна

. (1.23)

Для звука имеем дело с задачей, в которой звуковая волна распространяется вдоль жидкости и ее интенсивность падает с увеличением пройденного расстояния x. Очевидно, что это уменьшение будет происходить по закону, а для амплитуды как –, где коэффициент поглощения γ определяется посредством

.

Подставляем сюда (1.22) и (1.23), находим, таким образом, следующее выражение для коэффициента поглощения звука[10]:

, (1.24)

которое используется для расчёта объёмного эффекта затухания звуковых волн при верификации.

Глава 2 Методика расчета

1. Постановка задачи

Перейдём теперь непосредственно к постановке и решению задачи. Для этого рассмотрим цилиндрическую трубу (ротор), заполненную гексофторидом урана . Газ обладает теплопроводностью и молекулярной вязкостью. Ротор вращается с угловой скоростью ω. Его длина L намного больше, чем радиус r (L>> r), что позволяет считать ротор бесконечным (Рис.6.). Предполагаем, что температура T на внешней стенке постоянна и равна 300 K. Внутри ротора находится генератор звуковых волнc волновым вектором k направленным вдоль оси вращения. 

Рис.6. РоторНеобходимо разработать численный метод расчёта коэффициента затухания звуковых волн для вышеописанной модели и исследовать зависимости глубины проникновения звуковой волны от её волнового вектора, а также от радиуса и скорости вращения ротора.Фундамент исследования составила работа [14] в которой предложен метод верификации, основанный на полуаналитическом решении задачи о циркуляции газа в роторе бесконечной длины. Поставленная задача решается с гармоническим возмущением малой амплитуды во вращающемся газе. В работе также показано, как решение системы уравнений в частных производных сводится к решению системы однородных дифференциальных уравнений, которые могут быть решены почти с любой точностью на персональном компьютере.Запишем основную систему дифференциальных уравнений во вращающейся цилиндрической системе координат, описывающих движение в роторе [10]:

, (2.1)

, (2.2)

, (2.3)

(2.4)

(2.5)

Плотность и давление, подчиняются следующим распределениям:

, (2.6)

, (2.7)

где – давление и плотность на стенке ротора, соответственно,

образуется система уравнений, которая численно решается с помощью Maple при граничных условиях скользящей стенки:

=(2.8)

=0

и граничных условиях трения на стенке:

. (2.9)

Сравнение результатов, полученных с помощью данной полуаналитической модели и результатов численного моделирования, полученных в среде ANSYS CFX, показывает, что результаты эквивалентны[14].