Ричард Фейнман - 4a. Кинетика. Теплота. Звук
Мы, естественно, хотим описать поведение газа в масштабе, большем, чем длина свободного пробега, так что свойства газа не будут определяться поведением отдельных молекул. Например, смещение есть смещение центра инерции небольшого объема газа, а давление или плотность относятся к этому же объему. Мы обозначим давление через Р, а плотность через r, причем обе величины будут функциями от х и t. Необходимо помнить, что наше описание приближенное и справедливо лишь, когда свойства газа не слишком быстро меняются с расстоянием.
§ 3. Волновое уравнение
Итак, физические явления, происходящие в звуковой волне, обладают следующими тремя свойствами:
I. Газ движется, и плотность его меняется. II. При изменении плотности меняется и давление. III. Неравномерное распределение давления вызывает движение газа.
Рассмотрим сначала свойство П. Для любого газа, жидкости или твердого тела давление является функцией плотности. До прихода звуковой волны мы имели равновесное состояние с давлением Р0 и плотностью r. Давление Р зависит от плотности среды: Р=f(r), и в частности равновесное давление Р0=f(r0). Отклонения величины давления от равновесного в звуковой волне очень малы. Давление удобно измерять в барах (1 бар=105н/м2). Давление в одну стандартную атмосферу приблизительно равно 1 бар (1 атм=1,0133 бар). Для звука обычно используется логарифмическая шкала интенсивности, так как восприятие уха, грубо говоря, растет логарифмически. В этой децибельной шкале уровень звукового давления I связан с амплитудой звукового давления:
I=20log10(P/Pотн) дб, (47.1)
где давление отнесено к некоторому стандартному давлению Ротн=2·10-10 бар.
Звуковое давление Р=103 Ротн=2·10-7 бар соответствует довольно сильному звуку в 60 дб. Мы видим, что давление меняется в звуковой волне на очень малую величину по сравнению с равновесным или средним, равным 1 атм. Смещение и перепады плотности также очень малы. При взрывах, однако, изменения уже не столь малы; избыточное звуковое давление может превышать 1 атм. Такие большие перепады давления приводят к новым явлениям, которые мы рассмотрим позже. В звуковых волнах уровень силы звука выше 100 дб встречается редко; уровень силы звука в 120 дб уже вызывает боль в ушах. Поэтому, написав для звуковой волны
Р=Р0+Рu, r = r0+ru, (47.2)
можно считать, что изменение давления Puочень мало по сравнению с P0, а изменение плотности ru очень мало по сравнению с r0. Тогда
P0+Рu=f(r0+ru)=f(r0)+ ruf'(r0), (47.3)
где P0 = f(r0) и f'(r0) — производная от f(r), взятая при значении r =r0. Второе равенство здесь возможно только потому, что ru очень мало. Таким образом, мы находим, что избыточное давление Puпропорционально избыточной плотности ru; коэффициент пропорциональности обозначается через к:
(II) Рu=cru, где c=f'(r0)=(dP/dr)0. (47.4)
Это весьма простое соотношение и составляет точное содержание свойства II.
Перейдем теперь к свойству I. Предположим, что положение элемента объема воздуха, не возмущенного звуковой волной, есть х, а звук смещает его в момент времени t на величину c(х,t), так что его новое положение есть x+c(x,t), как показано на фиг. 47.3.
Фиг. 47.3. Смещение воздуха в точке х есть c (х,t), а в точке х+Dx равно c(x+Dx,t).
Первоначальный объем, приходящийся на единицу площади в плоской звуковой волне, есть Dx, а окончательный объем равен Dx+c(x+Dx,t)-c(x,t).
Далее, положение соседнего элемента объема есть х+Dx, и его смещенное положение есть х+Dx+c(х+Dx,t). Теперь можно найти изменение плотности. Поскольку мы рассматриваем плоскую волну, удобно взять единичную площадку, перпендикулярную оси х, т. е. направлению распространения волны. Количество воздуха, приходящееся на единичную площадку в интервале Dx, есть r0Dx, где r0 — невозмущенная, или равновесная, плотность воздуха. Эта порция воздуха, смещенная звуковой волной, будет находиться теперь между x+c (x,t) и x+Dх+c (х+Dx,t), причем количество воздуха в этом интервале то же самое, что в интервале Dx до прихода волны. Если через r обозначить новую плотность, то
r0Dx=r [x+Dx+c (x+Dx,t)-x-c (x,t)]. (47.5)
Поскольку Dx мало, можно написать c (x+Dx,t)-c (x,t)=(дc/дx)Dx. Здесь уже появляется частная производная, потому что c зависит и от x, и от времени. Наше уравнение принимает вид
r0Dx =r ((дc/дx) Dx +Dx), (47.6)
или
r0=(r0+ru)дc/дx+r0+ru. (47.7)
Но в звуковой волне все изменения малы, так что ru мало, c мало и дc/дх тоже мало. Поэтому в уравнении, которое мы только что написали,
ru=-r0(дc/дx)- ru(дc/дx), (47.8)
можно пренебречь ru(дc/дх) по сравнению с r0(дc/дх). Так мы приходим к соотношению, которое требовалось согласно свойству I:
(I) ru=-r0дc/дx. (47.9)
Именно такой вид уравнения можно было ожидать из чисто физических соображений. Если смещение различно для разных х, плотность будет изменяться. Знак тоже правильный: если смещение c растет с ростом х, так что воздух расширяется, плотность должна уменьшаться.
Теперь нам нужно найти третье уравнение — уравнение движения, производимого избытком давления. Зная соотношение между силой и давлением, можно получить уравнение движения. Возьмем объем воздуха толщиной Dx и с единичной площадью грани, перпендикулярной х, тогда масса воздуха в этом объеме есть r0Dx, а ускорение воздуха есть д2c/дt2, так что масса, умноженная на ускорение для этого слоя, есть r0Dx(д2c/дt2). (Если Dx; мало, то безразлично, где брать ускорение — на краю слоя или где-нибудь посредине.) Сила, действующая на единичную площадку нашего слоя, перпендикулярную оси x, должна быть равна r0Dx(д2хc/дt2). В точке х мы имеем силу Р(х,t), действующую на единицу площади в направлении +х, а в точке x+Dx; возникает сила в обратном направлении, по величине равная Р(x;+ Dx, t) (фиг. 47.4):
Фиг. 47.4. Результирующая сила в направлении оси х, возникающая за счет давления на единичную площадку, перпендикулярную к оси х, есть — (дР/дх)Dx.
Р(х, t)-P(x+Dx, t)=-(дP/дx) Dx=(дPu/дx) Dx. (47.10)
Мы учли, что Dx; мало и что только избыточное давление Ри меняется в зависимости от х. Итак, согласно свойству III мы получаем
(III) r0=д2c/дt2=-дPu/дx. (47.11)
Теперь уже уравнений достаточно, чтобы увязать все величины и привести к одной переменной, скажем х. Можно выразить Рuв (47.11) с помощью (47.4):
r0д2c/дt2-cдru/дx (47.12)
а затем исключить ru с помощью (I). Тогда r0 сократится и у нас останется
д2c/дt2=xд2c/дx2. (47.13)
Обозначим с2s =x, тогда можно написать
Это и есть волновое уравнение, которое описывает распространение звука в среде.
§ 4. Решения волнового уравнения
Посмотрим теперь, действительно ли волновое уравнение описывает основные свойства звуковых волн в среде. Прежде всего мы хотим вывести, что звуковое колебание, или возмущение, движется с постоянной скоростью. Кроме того, нам нужно доказать, что два различных колебания могут свободно проходить друг через друга, т. е. принцип суперпозиции. Мы хотим еще доказать, что звук может распространяться и вправо и влево. Все эти свойства должны содержаться в нашем одном уравнении.
Раньше мы отмечали, что любое возмущение, имеющее вид плоской волны и движущееся с постоянной скоростью, записывается в виде f(x-vt). Посмотрим теперь, является ли f(x-vt) решением волнового уравнения. Вычисляя дc/дх, получаем производную функции dcldx=f'(x-vt). Дифференцируя еще раз, находим