КАТЕГОРИИ: Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748) |
Поглощение звука в воздухе
Задачи, связанные с акустическим зондированием, требуют определения потерь энергии распространяющегося излучения. Наряду со сферической расходимостью вклад в ослабление излучения вносят процессы классического и молекулярного (релаксационного) поглощений, рассеяния на турбулентности и на частицах атмосферных образований. Для плоской монохроматической волны интенсивность I на некотором расстоянии R от источника может быть записана в следующем виде: I=I0exp(-aR), где I0 - начальная интенсивность звука, a- коэффициент ослабления звука, в общем случае равный, как уже говорилось, сумме коэффициентов поглощения и рассеяния. Физика процессов классического и молекулярного поглощений звука в воздухе достаточно хорошо разработана и освещена в литературе (см., например, [1, 73, 74]). Тем не менее имеющиеся эмпирические соотношения для коэффициентов молекулярного поглощения, полученные различными исследователями, часто имеют существенные расхождения, что приводит к неопределенности в выборе модели поглощения. Это становится особенно актуальным в задачах многочастотного акустического зондирования атмосферы, где малые ошибки в задании модели поглощения могут приводить к большим ошибкам в определении значения метеовеличин. Ниже для диапазона звуковых частот 0,1 < f < 30 кГц проводится сравнение имеющихся формул и экспериментальных данных, а также оценка вкладов различных механизмов поглощения звука в реальных атмосферных условиях [78]. Классическое поглощение звука связано с вязкими и термическими потерями. Вязкие потери обусловлены сдвиговыми напряжениями, возникающими между соседними участками газа с различными скоростями движения. Возникшие напряжения стремятся устранить неравенство скоростей. Процесс выравнивания заключается в обмене количеством движения между участками путем диффузии молекул, что приводит к потере энергии звуковой волны. Термические механизмы поглощения звука обусловлены процессом частичного выравнивания температур между участками сжатия с повышенной температурой и участками разряжения с пониженной температурой, а также с окружающей средой в целом, что приводит к потере волной некоторой доли энергии. Выравнивание температур между сжатиями и разряжениями происходит путем теплопроводности, а между участками с изменившейся температурой и средой - путем теплоизлучения. При распространении звука в атмосфере вклад теплоизлучения в ослабление незначителен и в дальнейшем не рассматривается. Влияние теплопроводности на поглощение звука сравнимо с влиянием сдвиговых напряжений. Классическое поглощение звука было впервые теоретически рассмотрено Стоксом [79] и развито Редеем [69] и Кирхгофом [80]. Коэффициент классического поглощения звука, обусловленного вязкостью и теплопроводностью воздуха, определяется известной формулой Стокса - Кирхгофа, которая после некоторых упрощений принимает вид [73, 75]
где h - вязкость воздуха, р - атмосферное давление, с - скорость звука, f - частота звука. Для обычных атмосферных условий [76, 77]
Подробные экспериментальные исследования поглощения звука в атмосфере осуществил Харрис [83, 84], проведя измерения в третьоктавных частотных интервалах в диапазоне частот от 2 до 12,5 кГц и температур от -0,5 до 25,1°С при нормальном атмосферном давлении. Результаты были экстраполированы на частотный диапазон 125 Гц Согласно термодинамическому принципу равномерного распределения энергии по степеням свободы возбужденной системы звуковая волна при большом числе межмолекулярных столкновений теряет часть своей кинетической энергии за счет насыщения внутренних степеней свободы многоатомных молекул. Скорость насыщения, т. е. установление термодинамического равновесия, определяется величиной t - временем релаксации.
Степень влияния релаксационного процесса на величину кинетической энергии звуковой волны зависит от соотношения между периодом волны t и временем релаксации t данного многоатомного газа. Максимум потерь кинетической энергии достигается при t = t [74]. Изложенная модель релаксационного процесса относится к однотипным многоатомным молекулам с одной внутренней степенью свободы. Реальные газы могут иметь сразу несколько степеней свободы, скорости насыщения которых характеризуются определенными временами релаксации t. В частности, для воздуха, являющегося смесью различных многоатомных газов, следует учитывать колебательные и вращательные степени свободы молекул [73], причем величина t колебательной релаксации гораздо больше, чем вращательной. Рассмотрим вращательную релаксацию. Обусловленный ею коэффициент поглощения с учетом малости времени релаксации записывается для диапазона частот f < 30 кГц в виде [75]
а суммарный коэффициент классического поглощения и поглощения за счет вращательной релаксации в виде [75, 85]
где p0=1014 мбар, Т0=293,15 К. Более простая зависимость приведена в [86].
Несмотря на кажущееся отличие выражений (1.18) и (1.19), расчеты по ним в диапазоне температур от - 40 до +40°С различаются не более чем на 10%. В случае колебательной релаксации для указанного диапазона частот основной вклад в поглощение звука в воздухе вносят механизмы колебательной релаксации молекул 02, N2, C02 и Н2О [73], причем главной особенностью молекул Н20 является то, что величина их концентрации существенно влияет на время релаксации каждого газа. Как следует из [86], суммарный коэффициент поглощения звука за счет колебательной релаксации смеси различных газов записывается в виде
где
Несколько иной вид коэффициента aКОЛ приведен в [77]:
где mm, lm, fm относятся к кислороду. Рассмотрим составляющие суммы (1.20). Коэффициент mm, определяющий максимальное поглощение звука, наступающее при f= fm, для двухатомных газов записывается в виде [86]
где N - молярная концентрация газа, св - удельная колебательная теплоемкость, R - универсальная газовая постоянная. В свою очередь,
где Q - характеристическая температура газа, Т - абсолютная температура воздуха. Колебательная релаксация кислорода. Учитывая, что в атмосфере 02 содержится в количестве N»0,21, а его характеристическая температура Q=2239 К, на основании (1.23) и (1-24) получим
Сведение mm из [77] к форме (1.25) дает
т. е. отклонение (1.25) от (1.26) составляет примерно 4%. Коэффициент
Здесь следует заметить, что в формуле (3.43) работы [74]
де fm - 3,06h1/3×104 Гц, h — молярная концентрация водяного пара (отношение парциального давления водяного пара к полному давлению, %), функциональная связь aмол с (f/fm) ошибочна, так как из нее следует уменьшение поглощения с увеличением частоты. Однако это противоречит экспериментальным результатам, в частности приведенным в [74] на той же странице, что и формула. Предполагая, что в знаменателе формулы сделана опечатка и следует читать (fm/f)2, ее можно привести к виду, аналогичному формулам других авторов. Проделав преобразования, получим
где mm= 0,625×10-5 × с. При Т = 293 К и соответственно с = 343,2 м/с коэффициент mm= 2,15× 10-3. Для этой же температуры по формуле (1.25) mm-2,11 × 10-3. Совпадение хорошее. Следует, однако, отметить, что mmиз [74] зависит от температуры как Т1/2 в то время как у других авторов mmимеет несколько иные зависимости. Сравнивая значения Положение максимума ослабления на частотной оси, т. е. выполнение равенства f=fm, приводящего к j (f/fm)=1, зависит от времени колебательной релаксации кислорода, величина которого в наибольшей степени определяется концентрацией паров воды [73]. Анализу tо2 и соответственно
из [83], где а =2,17рh/Т;
из [87];
из [1], где Т0=293 К, р0=1014 мбар;
из [77], где Т0=288 К;
из [92] при Т=298 К и р=1014 мбар.
расширение границ ее частотного диапазона приводит к весьма существенному отклонению от других экспериментальных данных. Аналогично ведет себя кривая, отвечающая формуле (1.32), по крайней мере в высокочастотной области. Определим величину расхождения в оценке
и вычислим для различных h и f величину
где a берется в виде (1.20), (1.21) с той или иной формой Применение экспериментальных результатов работы [85] приводит к значительному уменьшению величины Колебательная релаксация азота. Учитывая, что концентрация азота в воздухе N»0,78, а его характеристическая температура 9=3353 К, из уравнений (1.23) и (1.24) получим
Коэффициент максимального поглощения звука молекулами N2 примерно на порядок меньше аналогичного коэффициента для O2 (см. рис. 1.3). Результаты экспериментальных исследований [88], скорректированные на концентрацию азота в воздухе и приведенные к энергетическим величинам, имеют хорошее согласие с выражением (1.35). В то же время результаты экспериментальных исследований релаксационной частоты Релаксационная частота, согласно [87],
а по результатам работы [88]
Из [89], после учета концентрации N2 в воздухе и влияния других газовых компонентов, получаем
для диапазона температур от 311 до 418 К. Сравнивая значения релаксационных частот Для оценки различия в значениях
где a[87] и a[88] - коэффициенты поглощения за счет колебательной релаксации азота при
Колебательная релаксация паров воды. Для паров воды в воздухе коэффициент максимального поглощения за счет колебательной релаксации [86] имеет вид
Характеристические температуры Н20 и 02 совпадают, и на основании (1.24) и (1.39) получаем
Зависимость Поскольку релаксационные частоты Н20 и 02 совпадают [86], проведенный ранее анализ f о2 полностью применим и к f н2о Оценки величины Колебательная релаксация углекислого газа. В связи с тем, что концентрация СO2 в атмосфере существенно меньше других газов (порядка 3 × 10-2%), то и вклад его в ослабление звука пренебрежим. Однако с уменьшением концентрации H2O (h £ 0,1%) наличие СО2 заметно сказывается на величинах релаксационных частот основных газов: О2 и N2 (особенно О2). Так, при h = 0 % и отсутствии СО2 релаксационная частота кислорода
Рис. 1.9. Области пренебрежения компонентом aКВ с погрешностью dА1. Т=293 К, р=ро=1014 мбар, — — — релаксационная частота
Зависимость релаксационных частот О2 и Na от углекислого газа, а также от реакций обмена энергией при межмолекулярных столкновениях типа О2/О2, N2/N2, О2/N2 и др. учитывается первыми слагаемыми уравнений (1.34), (1.36), (1.37). Оценим существенность вкладов в общее поглощение отдельных его составляющих для реальных атмосферных условий. Прежде всего на координатной сетке { f, h }выделим области, где с той или иной точностью можно пренебречь вкладом от классического поглощения и поглощения за счет вращательной релаксации, т. е. пренебречь коэффициентом aКВ Для этого оценим величину
как функцию частоты и молярной концентрации паров воды h, где a - коэффициент полного поглощения, A1=aКВ/aкол, aкол= Установим теперь области частот и влажностей, в которых можно пренебречь суммарным эффектом колебательных релаксаций N2 и Н20. Вычислим для этого d А2=(1- 1 / А2)×100%, где А2 = aкол/
Проводя дальнейший анализ, можно выделить области f и h, где преобладают те или иные механизмы поглощения звука. Так, на рис. 1.11 выделены области по уровню 64=5%, а на рис. 1.12 - по уровню 64=10%. Очевидно, что с увеличением допустимой ошибки dA расширяются области, связанные с
Рис. 1.14. Нормализованная кривая (———), полученная на основе экспериментальных данных [83]; о о о - вычисления согласно (1.41); — — — вычисления на основе (1.22). Для сравнения влияния температуры на изменения указанных областей приведен рис. 1.13, на котором нанесены уровни dА=5% при Т=273 К. Сужение диапазона h обусловлено понижением порога насыщения; h» 0,6 % соответствует относительной влажности 100%. В данном случае изменение Т существенно сказывается на ошибке dА1 расширяя диапазон необходимого учета компонента aКВ в области h > 0,2% и f > 0,5 кГц. Проведем сравнение с экспериментальными данными Харриса [83]. Для этого рассмотрим коэффициент ослабления на длине волны, определяемый как m= a×l, и на основе уравнений (1.25), (1.34), (1.35), (1.37), с учетом вкладов в поглощение только O2 и N2, составим отношение
Нанесем на нормализованную кривую, приведенную в [83], рассчитанные по (1.41) значения Q = Рассмотрим одиночный релаксационный процесс. Тогда согласно (1.20), (1.21)
где am=mm/lm - максимальное поглощение на единицу расстояния на релаксационной частоте fm. Из (1.42) следует, что
Для сравнения с нормализованной кривой [83] преобразуем выражение (1.43) к виду
При fm® 0 j (f / fm)® 2, а при fm ® ¥ j (f / fm)® 0. Предположим, что оценка поглощения проводится на фиксированной частоте f и меняется молярная концентрация паров воды h. Тогда при некоторой влажности hm выполняется равенство f = f m. Поскольку fm ~ h, то fm / f ~ h / hm. Функциональная зависимость j от fm / f приведена на рис. 1.15. Если же составить отношение m/mm, то получим уравнение
Функциональная зависимость y от fm / f также приведена на рис. 1.15. очевидно, что именно эта форма, т. е. m/mm=y (fm / f), отвечает поведению экспериментальной кривой, а не уравнение (1.43), указывая на неточность, допущенную в [83] при построении нормализованной кривой. Рассматривая в дальнейшем зависимость m/mm от h / hm аналогична зависимости m/mm от u / u m (приведенной в [83]), где u и u m – относительные влажности, поскольку u / u m= h / h m
Вернемся к уравнению (1.41). Анализ этого уравнения показывает, что существуют два максимума функции Q(h/hm). Первый максимум обусловлен резонансом при f = f о2 (h = hm), а второй - при f = Проведенные оценки позволяют сделать вывод, что при h/hm <20 и 0,3 < f < 15 кГц экспериментальная кривая [83] хорошо описывает процесс поглощения звука в воздухе. Однако продолжение этой кривой в область hlhm > 20, как сделано в работе [84], некорректно. Кроме того, расширение частотного диапазона применимости экспериментальной кривой (измерения проводились в диапазоне от 2 до 12,5 кГц), как сделано в работе [83], может привести к существенным ошибкам в оценке поглощения звука, обусловленным влиянием колебательной релаксации N2. Для сравнения на рис. 1.14 (при f =1 кГц) нанесена кривая, отвечающая уравнению (1.22). Она асимптотически стремится к Корректный учет колебательной релаксации N2 обусловливает изменение температурных зависимостей поглощения звука, приведенных в работе [83]. Проведем качественный анализ aкол (Т) на основе уравнений (1.20), (1.25) и (1.35). При постоянной относительной влажности и молярная концентрация паров воды h растет с ростом температуры, что приводит к прямой зависимости f m от Т. Для одиночного релаксационного процесса существует один экстремум-максимум функции aкоп(T), который достигается при f = f m. В случае двух релаксирующих газов и аддитивности их вкладов в поглощение, а также несовпадении релаксационных частот при одинаковой температуре, должны наблюдаться два максимума функции aкол(Т). В частности, максимум
имеют отличное совпадение. Однако с ростом температуры наблюдается явное отличие функциональной связи коэффициента a(Т) из работы [83] от полученной выше. Это отличие обусловлено асимптотическим продолжением нормализованной кривой и область высоких значений h/hm без учета максимума В заключение отметим, что проведенные исследования позволили уточнить формы основных соотношений, необходимых для расчета коэффициентов поглощения звука в воздухе. Сравнение различных механизмов поглощения показало возможность упрощения оценок коэффициента поглощения в широком диапазоне звуковых частот для реальных атмосферных условий. Например, для приближенных вычислений, не требующих высокой точности учета поглощения, можно пользоваться соотношением (1.22). Однако более тонкие измерения, проводимые, в частности, с целью восстановления температуры и влажности на основе оценки коэффициента ослабления, требуют детального учета различных механизмов поглощения. Такое требование обусловило необходимость анализа приведенных в работах [83, 84] экспериментальных результатов и полученных на их основе функциональных связей коэффициента поглощения с температурой и влажностью, поскольку большое число исследователей использует именно эти результаты, что может привести и приводит к ряду существенных неточностей.
Дата добавления: 2015-07-02; Просмотров: 3772; Нарушение авторских прав?; Мы поможем в написании вашей работы! |