In physics, the acoustic wave equation governs the propagation of acoustic waves through a material medium. The form of the equation is a second order partial differential equation. The equation describes the evolution of acoustic pressure p or particle velocity u as a function of position r and time t. A simplified form of the equation describes acoustic waves in only one spatial dimension, while a more general form describes waves in three dimensions.
The acoustic wave equation was an important point of reference in the development of the electromagnetic wave equation in Kelvin's master class at Johns Hopkins University.^{[1]}
For lossy media, more intricate models need to be applied in order to take into account frequencydependent attenuation and phase speed. Such models include acoustic wave equations that incorporate fractional derivative terms, see also the acoustic attenuation article or the survey paper.^{[2]}
Contents

In one dimension 1

Equation 1.1

Solution 1.2

Derivation 1.3

In three dimensions 2

Equation 2.1

Solution 2.2

Cartesian coordinates 2.2.1

Cylindrical coordinates 2.2.2

Spherical coordinates 2.2.3

See also 3

References 4
In one dimension
Equation
Richard Feynman^{[3]} derives the wave equation that describes the behaviour of sound in matter in one dimension (position x) as:

{ \partial^2 p \over \partial x ^2 }  {1 \over c^2} { \partial^2 p \over \partial t ^2 } = 0
where p is the acoustic pressure (the local deviation from the ambient pressure), and where c is the speed of sound.
Solution
Provided that the speed c is a constant, not dependent on frequency (the dispersionless case), then the most general solution is

p = f(c t  x) + g(c t + x)
where f and g are any two twicedifferentiable functions. This may be pictured as the superposition of two waveforms of arbitrary profile, one (f) travelling up the xaxis and the other (g) down the xaxis at the speed c. The particular case of a sinusoidal wave travelling in one direction is obtained by choosing either f or g to be a sinusoid, and the other to be zero, giving

p=p_0 \sin(\omega t \mp kx).
where \omega is the angular frequency of the wave and k is its wave number.
Derivation
Derivation of the acoustic wave equation
The wave equation can be developed from the linearized onedimensional continuity equation, the linearized onedimensional force equation and the equation of state.
The equation of state (ideal gas law)

PV=nRT
In an adiabatic process, pressure P as a function of density \rho can be linearized to

P = C \rho \,
where C is some constant. Breaking the pressure and density into their mean and total components and noting that C=\frac{\partial P}{\partial \rho}:

P  P_0 = \left(\frac{\partial P}{\partial \rho}\right) (\rho  \rho_0).
The adiabatic bulk modulus for a fluid is defined as

B= \rho_0 \left(\frac{\partial P}{\partial \rho}\right)_{adiabatic}
which gives the result

PP_0=B \frac{\rho  \rho_0}{\rho_0}.
Condensation, s, is defined as the change in density for a given ambient fluid density.

s = \frac{\rho  \rho_0}{\rho_0}
The linearized equation of state becomes

p = B s\, where p is the acoustic pressure (PP_0).
The continuity equation (conservation of mass) in one dimension is


\frac{\partial \rho}{\partial t} + \frac{\partial }{\partial x} (\rho u) = 0.
Where u is the flow velocity of the fluid. Again the equation must be linearized and the variables split into mean and variable components.

\frac{\partial}{\partial t} ( \rho_0 + \rho_0 s) + \frac{\partial }{\partial x} (\rho_0 u + \rho_0 s u) = 0
Rearranging and noting that ambient density does not change with time or position and that the condensation multiplied by the velocity is a very small number:

\frac{\partial s}{\partial t} + \frac{\partial }{\partial x} u = 0
Euler's Force equation (conservation of momentum) is the last needed component. In one dimension the equation is:

\rho \frac{D u}{D t} + \frac{\partial P}{\partial x} = 0,
where D/Dt represents the convective, substantial or material derivative, which is the derivative at a point moving with medium rather than at a fixed point.
Linearizing the variables:

(\rho_0 +\rho_0 s)\left( \frac{\partial }{\partial t} + u \frac{\partial }{\partial x} \right) u + \frac{\partial }{\partial x} (P_0 + p) = 0.
Rearranging and neglecting small terms, the resultant equation becomes the linearized onedimensional Euler Equation:

\rho_0\frac{\partial u}{\partial t} + \frac{\partial p}{\partial x} = 0.
Taking the time derivative of the continuity equation and the spatial derivative of the force equation results in:

\frac{\partial^2 s}{\partial t^2} + \frac{\partial^2 u}{\partial x \partial t} = 0

\rho_0 \frac{\partial^2 u}{\partial x \partial t} + \frac{\partial^2 p}{\partial x^2} = 0.
Multiplying the first by \rho_0, subtracting the two, and substituting the linearized equation of state,

 \frac{\rho_0 }{B} \frac{\partial^2 p}{\partial t^2} + \frac{\partial^2 p}{\partial x^2} = 0.
The final result is

{ \partial^2 p \over \partial x ^2 }  {1 \over c^2} { \partial^2 p \over \partial t ^2 } = 0
where c = \sqrt{ \frac{B}{\rho_0 }} is the speed of propagation.
In three dimensions
Equation
Feynman^{[3]} provided a derivation of the wave equation that describes the behaviour of sound in matter in three dimensions as:

\nabla ^2 p  {1 \over c^2} { \partial^2 p \over \partial t ^2 } = 0
where \nabla ^2 is the Laplace operator, p is the acoustic pressure (the local deviation from the ambient pressure), and where c is the speed of sound.
A similar looking wave equation but for the vector field particle velocity is given by

\nabla ^2 \mathbf{u}\;  {1 \over c^2} { \partial^2 \mathbf{u}\; \over \partial t ^2 } = 0 .
In some situations, it is more convenient to solve the wave equation for an abstract scalar field velocity potential which has the form

\nabla ^2 \Phi  {1 \over c^2} { \partial^2 \Phi \over \partial t ^2 } = 0
and then derive the physical quantities particle velocity and acoustic pressure by the equations (or definition, in the case of particle velocity):

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

p = \rho {\partial \over \partial t}\Phi .
Solution
The following solutions are obtained by separation of variables in different coordinate systems. They are phasor solutions, that is they have an implicit timedependence factor of e^{i\omega t} where \omega = 2 \pi f is the angular frequency. The explicit time dependence is given by

p(r,t,k) = \operatorname{Real}\left[p(r,k) e^{i\omega t}\right]
Here k = \omega/c \ is the wave number.
Cartesian coordinates

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

p(r,k)=AH_0^{(1)}(kr) + \ BH_0^{(2)}(kr).
where the asymptotic approximations to the Hankel functions, when kr\rightarrow \infty , are

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)}.
Spherical coordinates

p(r,k)=\frac{A}{r}e^{\pm ikr}.
Depending on the chosen Fourier convention, one of these represents an outward travelling wave and the other an unphysical inward travelling wave. The inward travelling solution wave is only unphysical because of the singularity that occurs at r=0; inward travelling waves do exist.
See also
References

^ William Thomson, Lord Kelvin (1904) Lecture notes taken by A. S. Hathaway, Molecular dynamics and the wave theory of light, page 80, Cambridge University Press, link from Internet Archive

^ S. P. NÃ¤sholm and S. Holm, "On a Fractional Zener Elastic Wave Equation," Fract. Calc. Appl. Anal. Vol. 16, No 1 (2013), pp. 2650, DOI: 10.2478/s1354001300031 Link to eprint

^ ^{a} ^{b} Richard Feynman, Lectures in Physics, Volume 1, 1969, Addison Publishing Company, Addison
This article was sourced from Creative Commons AttributionShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, EGovernment Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a nonprofit organization.