КАТЕГОРИИ: Архитектура-(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]
Исследования Кнудсена [81] продемонстрировали сильную зависимость измеренной величины ослабления звука от температуры и влажности воздуха. Он определил, что экспериментально измеренные величины на один-два порядка превышают предсказываемые теорией значения для коэффициента классического поглощения, обусловленного вязкостью и теплопроводностью воздуха. Кнезер [82] объяснил это избыточное поглощение процессами колебательной релаксации молекул воздуха. Это предположение подтвердилось дальнейшими экспериментами, и поглощение звука, обусловленное релаксационными процессами в многоатомных газах, получило название молекулярного поглощения. Оно формально связано вязкими потерями, поскольку вызвано периодическими изменениями сжатия среды при распространении акустической волны. Поэтому его иногда называют объемной или второй вязкостью. Подробные экспериментальные исследования поглощения звука в атмосфере осуществил Харрис [83, 84], проведя измерения в третьоктавных частотных интервалах в диапазоне частот от 2 до 12,5 кГц и температур от -0,5 до 25,1°С при нормальном атмосферном давлении. Результаты были экстраполированы на частотный диапазон 125 Гц 12,5 кГц и на диапазон изменения температур - 10 30°С. На рис. 1.2 приведены данные, описанные выше. На основе большого числа экспериментальных результатов были сделаны следующие выводы: 1) коэффициент поглощения звука растет с ростом температуры; 2) максимум на кривой, изображающей зависимость коэффициента поглощения, сдвигается в сторону уменьшения относительной влажности с ростом температуры. Согласно термодинамическому принципу равномерного распределения энергии по степеням свободы возбужденной системы звуковая волна при большом числе межмолекулярных столкновений теряет часть своей кинетической энергии за счет насыщения внутренних степеней свободы многоатомных молекул. Скорость насыщения, т. е. установление термодинамического равновесия, определяется величиной 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], суммарный коэффициент поглощения звука за счет колебательной релаксации смеси различных газов записывается в виде
где — коэффициент максимального поглощения звука в i -м газе на длине волны = c/fmi, fmi =l/(2pt i) - релаксационная частота, определяемая временем релаксации t i. Функция находится из уравнения
Несколько иной вид коэффициента 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%. Коэффициент определяемый из соотношений работы [74], имеет иной вид:
Здесь следует заметить, что в формуле (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имеет несколько иные зависимости. Сравнивая значения рассчитанные по уравнениям (1.25) - (1.27), с экспериментальными данными работы [83], можно сделать вывод, что уравнения (1.25) и (1.26) достаточно точно определяют зависимость коэффициента максимального поглощения от температуры, а (1.27) лишь при Т» 293 К совпадает с экспериментом (рис. 1.3). Положение максимума ослабления на частотной оси, т. е. выполнение равенства f=fm, приводящего к j (f/fm)=1, зависит от времени колебательной релаксации кислорода, величина которого в наибольшей степени определяется концентрацией паров воды [73]. Анализу tо2 и соответственно посвящено большое число как экспериментальных, так и теоретических работ, достаточно обширный перечень которых приведен в [85]. Кроме того, имеется стандарт [87] на вычисление (далее будем обозначать просто ). Однако, несмотря на это, нельзя дать однозначный ответ относительно абсолютно надежного выбора Функциональной связи с концентрацией водяного пара, поскольку наблюдаются большие различия, особенно в областях низких и повышенных влажностей (рис. 1.4). Рисунок графически отражает зависимости:
из [83], где а =2,17рh/Т;
из [87];
из [1], где Т0=293 К, р0=1014 мбар;
из [77], где Т0=288 К;
из [92] при Т=298 К и р=1014 мбар.
расширение границ ее частотного диапазона приводит к весьма существенному отклонению от других экспериментальных данных. Аналогично ведет себя кривая, отвечающая формуле (1.32), по крайней мере в высокочастотной области. Определим величину расхождения в оценке (далее будем писать просто , и т. д.) при использовании различных форм связи с молярной концентрацией водяного пара h. Для этого возьмем за основу стандарт [87]:
и вычислим для различных h и f величину где a берется в виде (1.20), (1.21) с той или иной формой - по стандарту [87]. На рис. 1.5 представлены результаты вычислений при использовании в форме (1.32). Здесь большие , особенно при h> 1%. Применение экспериментальных результатов работы [85] приводит к значительному уменьшению величины не превышающей ±10% при 0,03 < h < 2,5% и +10 -20% в диапазоне 0,01 <h< 6%, что хорошо видно на рис. 1.6. Таким образом, анализ показал, что для оценки следует использовать в форме (1.25), в форме (1.34). Колебательная релаксация азота. Учитывая, что концентрация азота в воздухе N»0,78, а его характеристическая температура 9=3353 К, из уравнений (1.23) и (1.24) получим
Коэффициент максимального поглощения звука молекулами N2 примерно на порядок меньше аналогичного коэффициента для O2 (см. рис. 1.3). Результаты экспериментальных исследований [88], скорректированные на концентрацию азота в воздухе и приведенные к энергетическим величинам, имеют хорошее согласие с выражением (1.35). В то же время результаты экспериментальных исследований релаксационной частоты , приведенные в работах [88—90], показывают, что имеющийся стандарт [87] на вычисление как функции h дает завышенные значения частоты. Как и в случае колебательной релаксации кислорода, релаксационная частота азота сильно зависит от степени концентрации паров воды в воздухе. Релаксационная частота, согласно [87],
а по результатам работы [88]
Из [89], после учета концентрации N2 в воздухе и влияния других газовых компонентов, получаем
для диапазона температур от 311 до 418 К. Сравнивая значения релаксационных частот (h), вычисленные по (1.36) и (1.37), с экспериментальными результатами работ [89, 90], можно сделать вывод, что соотношение (1.37) более точно определяет связь с h (рис. 1.7). Результаты работы [90] практически совпадают с таковыми из [88] и на рисунок не нанесены. Для оценки различия в значениях определяемых при разных выражениях для вычислим величину
где a[87] и a[88] - коэффициенты поглощения за счет колебательной релаксации азота при соответственно из работ [87] и [88]. На рис. 1.8 нанесены результаты этих вычислений. Видно, что несовпадение релаксационных частот в [87] и [88] приводит к очень большому отклонению в оценке во всем диапазоне влажностей. При частоте звука f >2 кГц вычисления по стандарту [87] дают завышенные значения . При меньших частотах возможны как завышенные, так и заниженные значения .
Колебательная релаксация паров воды. Для паров воды в воздухе коэффициент максимального поглощения за счет колебательной релаксации [86] имеет вид
Характеристические температуры Н20 и 02 совпадают, и на основании (1.24) и (1.39) получаем
Зависимость для нескольких значений h приведена на рис. 1.3. Поскольку релаксационные частоты Н20 и 02 совпадают [86], проведенный ранее анализ f о2 полностью применим и к f н2о Оценки величины свидетельствуют о том, что она может меняться от пренебрежимо малой до сравнимой с поглощением на O2 при высоких значениях h. Колебательная релаксация углекислого газа. В связи с тем, что концентрация СO2 в атмосфере существенно меньше других газов (порядка 3 × 10-2%), то и вклад его в ослабление звука пренебрежим. Однако с уменьшением концентрации H2O (h £ 0,1%) наличие СО2 заметно сказывается на величинах релаксационных частот основных газов: О2 и N2 (особенно О2). Так, при h = 0 % и отсутствии СО2 релаксационная частота кислорода 6Гц, а при наличии СО2 может достигать нескольких Десятков и даже сотен герц [73].
Рис. 1.9. Области пренебрежения компонентом aКВ с погрешностью dА1. Т=293 К, р=ро=1014 мбар, — — — релаксационная частота согласно (1.34).
Зависимость релаксационных частот О2 и Na от углекислого газа, а также от реакций обмена энергией при межмолекулярных столкновениях типа О2/О2, N2/N2, О2/N2 и др. учитывается первыми слагаемыми уравнений (1.34), (1.36), (1.37). Оценим существенность вкладов в общее поглощение отдельных его составляющих для реальных атмосферных условий. Прежде всего на координатной сетке { f, h }выделим области, где с той или иной точностью можно пренебречь вкладом от классического поглощения и поглощения за счет вращательной релаксации, т. е. пренебречь коэффициентом aКВ Для этого оценим величину как функцию частоты и молярной концентрации паров воды h, где a - коэффициент полного поглощения, A1=aКВ/aкол, aкол= + + . Воспользуемся уравнениями (1.19), (1.20), (1.25), (1.34), (1.35), (1.37), (1.40) при Т = 293 К и р= 1014 мбар и составим сетку значений d А1 в координатах f и h. На рис. 1.9 приведено несколько уровней d А1, которые определяют искомые области. Так, для рабочих частот, меньших 1 кГц, пренебрежение a КВприводит к ошибкам менее 5%. Максимум функции f (h, d А1 = const) обусловлен тем, что при заданных условиях частота звуковых колебаний совпадает с релаксационной частотой . Далее идет спад, обусловленный уменьшением . Затем функция f (h, d А1)вновь начинает расти, поскольку возрастает вклад колебательной релаксации N2 в aкол. Таким образом, при заданной рабочей частоте f и известной молярнойконцентрации паров воды h можно оценить ошибку неучета классического поглощения и вращательной релаксации для данной температуры, используя рис. 1.9. Установим теперь области частот и влажностей, в которых можно пренебречь суммарным эффектом колебательных релаксаций N2 и Н20. Вычислим для этого d А2=(1- 1 / А2)×100%, где А2 = aкол/ =1+ / + / . Фиксируя уровни d А2 и нанося их на сетку координат { f, h }, определяем области возможного пренебрежения вкладами и (рис. 1,10).
Проводя дальнейший анализ, можно выделить области f и h, где преобладают те или иные механизмы поглощения звука. Так, на рис. 1.11 выделены области по уровню 64=5%, а на рис. 1.12 - по уровню 64=10%. Очевидно, что с увеличением допустимой ошибки dA расширяются области, связанные с и , уменьшается влияние aКВ и становится возможным пренебрежение компонентом в рассмотренном диапазоне f и h. Рис. 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 = / в зависимости от h / hm где hm подразумевает равенство при h= hm (рис. 1.14). Однако предварительно следует сделать замечание, что по оси ординат нормализованной кривой необходимо откладывать именно m/mm, а не a/a т, как указано в [83] (фактически m/mm там и используется). Это вытекает из проведенного ниже анализа зависимостей aкол и m кол. Рассмотрим одиночный релаксационный процесс. Тогда согласно (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 = . Вычисляя Qи сравнивая с экспериментальной кривой работы [83] (с учетом экстраполяции на более высокие значения h/hm [84]), убеждаемся, что существует различие в поведении сравниваемых зависимостей. В пределах разброса экспериментальных данных [83] результаты хорошо совпадают при h/hm <20 и относительно высоких частотах. Однако с увеличением h/hm экспериментальная кривая и результаты вычислений начинают расходиться, причем расхождение усиливается с уменьшением частот (см. рис. 1.14). Проведенные оценки позволяют сделать вывод, что при h/hm <20 и 0,3 < f < 15 кГц экспериментальная кривая [83] хорошо описывает процесс поглощения звука в воздухе. Однако продолжение этой кривой в область hlhm > 20, как сделано в работе [84], некорректно. Кроме того, расширение частотного диапазона применимости экспериментальной кривой (измерения проводились в диапазоне от 2 до 12,5 кГц), как сделано в работе [83], может привести к существенным ошибкам в оценке поглощения звука, обусловленным влиянием колебательной релаксации N2. Для сравнения на рис. 1.14 (при f =1 кГц) нанесена кривая, отвечающая уравнению (1.22). Она асимптотически стремится к = 0,18 при h/hm> 8, хорошо совпадая как с экспериментальными данными [83], так и с вычислениями по уравнению (1.20) при h/hm < 8. Корректный учет колебательной релаксации N2 обусловливает изменение температурных зависимостей поглощения звука, приведенных в работе [83]. Проведем качественный анализ aкол (Т) на основе уравнений (1.20), (1.25) и (1.35). При постоянной относительной влажности и молярная концентрация паров воды h растет с ростом температуры, что приводит к прямой зависимости f m от Т. Для одиночного релаксационного процесса существует один экстремум-максимум функции aкоп(T), который достигается при f = f m. В случае двух релаксирующих газов и аддитивности их вкладов в поглощение, а также несовпадении релаксационных частот при одинаковой температуре, должны наблюдаться два максимума функции aкол(Т). В частности, максимум относительно максимума сдвинут к более высоким значениям температуры. Качественная картина поведения (Т) дана на рисунке 1.16. Если же учесть и эффекты, формирующие коэффициент aкв, то зависимость суммарного коэффициента поглощения звука a=aкол + aкв, особенно при высоких частотах, будет еще больше отличаться от зависимостей a(Т), приведенных в работе [83] (см. рис. 1.16). Результаты количественного сравнения для и = 40% и f =1 кГц (рис. 1.17) показывают, что при Т < 277 К экспериментальные и расчетные данные
имеют отличное совпадение. Однако с ростом температуры наблюдается явное отличие функциональной связи коэффициента a(Т) из работы [83] от полученной выше. Это отличие обусловлено асимптотическим продолжением нормализованной кривой и область высоких значений h/hm без учета максимума при (построенные в работе [83] кривые a(Т)получены с использованием нормализованной кривой m/mm). В заключение отметим, что проведенные исследования позволили уточнить формы основных соотношений, необходимых для расчета коэффициентов поглощения звука в воздухе. Сравнение различных механизмов поглощения показало возможность упрощения оценок коэффициента поглощения в широком диапазоне звуковых частот для реальных атмосферных условий. Например, для приближенных вычислений, не требующих высокой точности учета поглощения, можно пользоваться соотношением (1.22). Однако более тонкие измерения, проводимые, в частности, с целью восстановления температуры и влажности на основе оценки коэффициента ослабления, требуют детального учета различных механизмов поглощения. Такое требование обусловило необходимость анализа приведенных в работах [83, 84] экспериментальных результатов и полученных на их основе функциональных связей коэффициента поглощения с температурой и влажностью, поскольку большое число исследователей использует именно эти результаты, что может привести и приводит к ряду существенных неточностей.
Дата добавления: 2015-07-02; Просмотров: 3678; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |