Phương trình sóng âm

Trong vật lý học, phương trình sóng âm là phương trình vi phân riêng phần bậc hai chi phối sự truyền sóng âm thông qua dụng cụ vật chất tương ứng với trường sóng dừng. Phương trình này mô tả sự tiến triển của áp suất âm p hay vận tốc hạt u là hàm của li độ x và thời gian t. Dạng phương trình rút gọn (vô hướng) mô tả sóng âm chỉ trong một chiều không gian, trong khi dạng tổng quát mô tả sóng theo ba chiều. Sóng truyền theo hướng xác định trước cũng có thể được tính bằng phương trình sóng một chiều bậc nhất.

Trong không gian một chiều

Phương trình

Phương trình sóng mô tả một trường sóng dừng trong không gian một chiều (li độ $x$) là:

${\partial ^{2}p \over \partial x^{2}}-{1 \over c^{2}}{\partial ^{2}p \over \partial t^{2}}=0,$

Trong đó:

• $p$ là áp suất âm (độ lệch cục bộ so với áp suất xung quanh)
• $c$ là tốc độ âm thanh.

Nghiệm

Với điều kiện là tốc độ $c$ là một hằng số, không phụ thuộc vào tần số (trường hợp không tán sắc), thì nghiệm tổng quát nhất là:

$p=f(ct-x)+g(ct+x)$

Trong đó $f$ và $g$ là hai hàm đạo hàm hai lần bất kỳ (có thể được hình dung là sự chồng chất của hai dạng sóng có cấu hình tùy ý), một sóng ($f$) truyền lên trục x và sóng ($g$) kia truyền xuống trục dưới tốc độ $c$. Trường hơp cụ thể của sóng hình sin khi đang truyền theo một hướng thu được bằng cách chọn $f$ hoặc $g$ là một sóng hình sin, còn phương kia bằng 0, cho ra

$p=p_{0}\sin(\omega t\mp kx)$.

Trong đó $\omega$ là tần số góc của sóng và $k$ chính là số sóng của nó.

Trong không gian ba chiều

Phương trình

Feynman cung cấp một vi phân của phương trình sóng cho âm thanh ba chiều là:

$\nabla ^{2}p-{1 \over c^{2}}{\partial ^{2}p \over \partial t^{2}}=0,$

Trong đó $\nabla ^{2}$ chính là toán tử Laplace, $p$ là áp suất âm (độ lệch cục bộ so với áp suất xung quanh), và $c$ là tốc độ âm thanh.

Phương trình sóng tương tự nhưng đối với vận tốc hạt trong trường vectơ được cho bởi phương trình:

$\nabla ^{2}\mathbf {u} \;-{1 \over c^{2}}{\partial ^{2}\mathbf {u} \; \over \partial t^{2}}=0$.

Trong một số trường hợp, sẽ thuận tiện hơn khi giải phương trình sóng cho thế vận tốc trong trường vô hướng trừu tượng có dạng:

$\nabla ^{2}\Phi -{1 \over c^{2}}{\partial ^{2}\Phi \over \partial t^{2}}=0$

và sau đó rút ra các đại lượng vật lý về vận tốc hạt và áp suất âm theo các phương trình (hoặc định nghĩa, trong trường hợp vận tốc hạt):

$\mathbf {u} =\nabla \Phi \;$,

$p=-\rho {\partial \over \partial t}\Phi$.

Cách giải

Phương trình này có thể tìm được nghiệm bằng cách tách biến số trong các hệ tọa độ khác nhau. Các nghiệm này là phasơ, tức là chúng có hệ số phụ thuộc thời gian tiềm ẩn là $e^{i\omega t}$ (trong đó $\omega =2\pi f$ chính là tần số góc). Sự phụ thuộc thời gian rõ ràng được đưa ra bởi công thức:

$p(r,t,k)=\operatorname {Real} \left[p(r,k)e^{i\omega t}\right]$

Ở đây $k=\omega /c\$ chính là hệ số sóng.

Trong mặt phẳng tọa độ Descartes

$p(r,k)=Ae^{\pm ikr}$.

Trong mặt phẳng tọa độ trụ

$p(r,k)=AH_{0}^{(1)}(kr)+\ BH_{0}^{(2)}(kr)$.

Trong đó các phép xấp xỉ tiệm cận của hàm Hankel, khi $kr\rightarrow \infty$, là:

$H_{0}^{(1)}(kr)\simeq {\sqrt {\frac {2}{\pi kr}}}e^{i(kr-\pi /4)}$

$H_{0}^{(2)}(kr)\simeq {\sqrt {\frac {2}{\pi kr}}}e^{-i(kr-\pi /4)}$.

Trong mặt phẳng tọa độ cầu

$p(r,k)={\frac {A}{r}}e^{\pm ikr}$.

Tùy thuộc vào quy ước Fourier đã chọn, một trong những sóng này biểu thị sóng truyền ra ngoài và sóng còn lại biểu thị sóng phi vật lý truyền vào trong. Sóng truyền vào trong chỉ là phi vật lý do điểm kỳ dị xảy ra tại r=0; sóng truyền vào trong có tồn tại.