digplanet beta 1: Athena
Share digplanet:


Applied sciences






















Acoustic theory is the field relating to mathematical description of sound waves. It is derived from fluid dynamics. See acoustics for the engineering approach.

The propagation of sound waves in a fluid (such as water) can be modeled by an equation of continuity (conservation of mass) and an equation of motion (conservation of momentum) . With some simplifications, in particular constant density, they can be given as follows:

     \frac{\partial p}{\partial t} + \kappa~\nabla \cdot \mathbf{u} & = 0  \qquad \text{(Mass balance)} \\
     \rho_0 \frac{\partial \mathbf{u}}{\partial t} + \nabla p & = 0  \qquad \text{(Momentum balance)}

where p(\mathbf{x}, t) is the acoustic pressure and \mathbf{u}(\mathbf{x}, t) is the flow velocity vector, \mathbf{x} is the vector of spatial coordinates x, y, z, t is the time, \rho_0 is the static mass density of the medium and \kappa is the bulk modulus of the medium. The bulk modulus can be expressed in terms of the density and the speed of sound in the medium (c_0) as

\kappa = \rho_0 c_0^2 ~.

If the flow velocity field is irrotational,  \nabla \times \mathbf{u}=\mathbf{0}, then the acoustic wave equation is a combination of these two sets of balance equations and can be expressed as [1]

   \cfrac{\partial^2 \mathbf{u}}{\partial t^2} - c_0^2~\nabla^2\mathbf{u} = 0
   \qquad \text{or} \qquad
   \cfrac{\partial^2 p}{\partial t^2} - c_0^2~\nabla^2 p = 0,

where we have used the vector Laplacian,  \nabla^2 \mathbf{u} = \nabla(\nabla \cdot \mathbf{u}) - \nabla \times (\nabla \times \mathbf{u}) . The acoustic wave equation (and the mass and momentum balance equations) are often expressed in terms of a scalar potential \varphi where \mathbf{u} = \nabla\varphi. In that case the acoustic wave equation is written as

  \cfrac{\partial^2 \varphi}{\partial t^2} - c_0^2~\nabla^2 \varphi = 0

and the momentum balance and mass balance are expressed as

   p + \rho_0~\cfrac{\partial\varphi}{\partial t} = 0 ~;~~
   \rho + \cfrac{\rho_0}{c_0^2}~\cfrac{\partial\varphi}{\partial t} = 0 ~.

Derivation of the governing equations[edit]

The derivations of the above equations for waves in an acoustic medium are given below.

Conservation of momentum[edit]

The equations for the conservation of linear momentum for a fluid medium are

   \rho \left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + \nabla \cdot\boldsymbol{\tau} + \rho\mathbf{g}

where \mathbf{g} is the body force per unit mass, p is the pressure, and \boldsymbol{\tau} is the deviatoric stress. If \boldsymbol{\tau} is the Cauchy stress, then

   p := -\tfrac{1}{3}~\text{tr}(\boldsymbol{\tau}) ~;~~
   \boldsymbol{\tau} := \boldsymbol{\tau} + p~\boldsymbol{\mathit{1}}

where  \boldsymbol{\mathit{1}} is the rank-2 identity tensor.

We make several assumptions to derive the momentum balance equation for an acoustic medium. These assumptions and the resulting forms of the momentum equations are outlined below.

Assumption 1: Newtonian fluid[edit]

In acoustics, the fluid medium is assumed to be Newtonian. For a Newtonian fluid, the deviatoric stress tensor is related to the flow velocity by

\boldsymbol{\tau} = \mu~\left[\nabla\mathbf{u} + (\nabla\mathbf{u})^T\right] +
       \lambda~(\nabla \cdot \mathbf{u})~\boldsymbol{\mathit{1}}

where \mu is the shear viscosity and \lambda is the bulk viscosity.

Therefore, the divergence of \boldsymbol{\tau} is given by

   \nabla\cdot\boldsymbol{\tau} \equiv \cfrac{\partial s_{ij}}{\partial x_i} & = 
     \mu \left[\cfrac{\partial}{\partial x_i}\left(\cfrac{\partial u_i}{\partial x_j}+\cfrac{\partial u_j}{\partial x_i}\right)\right] + \lambda~\left[\cfrac{\partial}{\partial x_i}\left(\cfrac{\partial u_k}{\partial x_k}\right)\right]\delta_{ij} \\
    & = \mu~\cfrac{\partial^2 u_i}{\partial x_i \partial x_j} + \mu~\cfrac{\partial^2 u_j}{\partial x_i\partial x_i} + \lambda~\cfrac{\partial^2 u_k}{\partial x_k\partial x_j} \\
    & = (\mu + \lambda)~\cfrac{\partial^2 u_i}{\partial x_i \partial x_j} + \mu~\cfrac{\partial^2 u_j}{\partial x_i^2} \\
    & \equiv (\mu + \lambda)~\nabla(\nabla\cdot\mathbf{u}) + \mu~\nabla^2\mathbf{u} ~.

Using the identity \nabla^2\mathbf{u} = \nabla(\nabla\cdot\mathbf{u}) - \nabla\times\nabla\times\mathbf{u}, we have

  \nabla\cdot\boldsymbol{\tau} = (2\mu + \lambda)~\nabla(\nabla\cdot\mathbf{u}) -

The equations for the conservation of momentum may then be written as

   \rho \left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + (2\mu + \lambda)~\nabla(\nabla\cdot\mathbf{u}) -
      \mu~\nabla\times\nabla\times\mathbf{u} + \rho\mathbf{g}

Assumption 2: Irrotational flow[edit]

For most acoustics problems we assume that the flow is irrotational, that is, the vorticity is zero. In that case

   \nabla\times\mathbf{u} = 0

and the momentum equation reduces to

   \rho \left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + (2\mu + \lambda)~\nabla(\nabla\cdot\mathbf{u}) + \rho\mathbf{g}

Assumption 3: No body forces[edit]

Another frequently made assumption is that effect of body forces on the fluid medium is negligible. The momentum equation then further simplifies to

   \rho \left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + (2\mu + \lambda)~\nabla(\nabla\cdot\mathbf{u})

Assumption 4: No viscous forces[edit]

Additionally, if we assume that there are no viscous forces in the medium (the bulk and shear viscosities are zero), the momentum equation takes the form

   \rho \left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p

Assumption 5: Small disturbances[edit]

An important simplifying assumption for acoustic waves is that the amplitude of the disturbance of the field quantities is small. This assumption leads to the linear or small signal acoustic wave equation. Then we can express the variables as the sum of the (time averaged) mean field (\langle\cdot\rangle) that varies in space and a small fluctuating field (\tilde{\cdot}) that varies in space and time. That is

   p = \langle p\rangle + \tilde{p} ~;~~
   \rho = \langle\rho\rangle + \tilde{\rho} ~;~~
   \mathbf{u} = \langle\mathbf{u}\rangle + \tilde{\mathbf{u}}


   \cfrac{\partial\langle p \rangle}{\partial t} = 0 ~;~~ \cfrac{\partial\langle \rho \rangle}{\partial t} = 0 ~;~~
   \cfrac{\partial\langle \mathbf{u} \rangle}{\partial t} = \mathbf{0} ~.

Then the momentum equation can be expressed as

   \left[\langle\rho\rangle+\tilde{\rho}\right] \left[\frac{\partial\tilde{\mathbf{u}}}{\partial t}  + \left[\langle\mathbf{u}\rangle+\tilde{\mathbf{u}}\right] \cdot \nabla \left[\langle\mathbf{u}\rangle+\tilde{\mathbf{u}}\right]\right] = -\nabla \left[\langle p\rangle+\tilde{p}\right]

Since the fluctuations are assumed to be small, products of the fluctuation terms can be neglected (to first order) and we have

   \langle\rho\rangle~\frac{\partial\tilde{\mathbf{u}}}{\partial t} & +
   \left[\langle\rho\rangle+\tilde{\rho}\right]\left[\langle\mathbf{u}\rangle\cdot\nabla \langle\mathbf{u}\rangle\right]+
   \langle\rho\rangle\left[\langle\mathbf{u}\rangle\cdot\nabla\tilde{\mathbf{u}} +
     \tilde{\mathbf{u}}\cdot\nabla\langle\mathbf{u}\rangle\right] \\
   & = -\nabla \left[\langle p\rangle+\tilde{p}\right]

Assumption 6: Homogeneous medium[edit]

Next we assume that the medium is homogeneous; in the sense that the time averaged variables \langle p \rangle and \langle \rho \rangle have zero gradients, i.e.,

  \nabla\langle p \rangle =  0 ~;~~ \nabla\langle \rho \rangle = 0 ~.

The momentum equation then becomes

   \langle\rho\rangle~\frac{\partial\tilde{\mathbf{u}}}{\partial t} +
   \left[\langle\rho\rangle+\tilde{\rho}\right]\left[\langle\mathbf{u}\rangle\cdot\nabla \langle\mathbf{u}\rangle\right]+
   \langle\rho\rangle\left[\langle\mathbf{u}\rangle\cdot\nabla\tilde{\mathbf{u}} +
   = -\nabla\tilde{p}

Assumption 7: Medium at rest[edit]

At this stage we assume that the medium is at rest which implies that the mean flow velocity is zero, i.e. \langle\mathbf{u}\rangle = 0 . Then the balance of momentum reduces to

    \langle\rho\rangle~\frac{\partial\tilde{\mathbf{u}}}{\partial t} = -\nabla\tilde{p}

Dropping the tildes and using \rho_0 := \langle\rho\rangle, we get the commonly used form of the acoustic momentum equation

    \rho_0~\frac{\partial\mathbf{u}}{\partial t} + \nabla p = 0 ~.

Conservation of mass[edit]

The equation for the conservation of mass in a fluid volume (without any mass sources or sinks) is given by

\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0

where \rho(\mathbf{x},t) is the mass density of the fluid and \mathbf{u}(\mathbf{x},t) is the flow velocity.

The equation for the conservation of mass for an acoustic medium can also be derived in a manner similar to that used for the conservation of momentum.

Assumption 1: Small disturbances[edit]

From the assumption of small disturbances we have

   p = \langle p\rangle + \tilde{p} ~;~~
   \rho = \langle\rho\rangle + \tilde{\rho} ~;~~
   \mathbf{u} = \langle\mathbf{u}\rangle + \tilde{\mathbf{u}}


   \cfrac{\partial\langle p \rangle}{\partial t} = 0 ~;~~ \cfrac{\partial\langle \rho \rangle}{\partial t} = 0 ~;~~
   \cfrac{\partial\langle \mathbf{u} \rangle}{\partial t} = \mathbf{0} ~.

Then the mass balance equation can be written as

  \frac{\partial\tilde{\rho}}{\partial t} +
    \left[\langle\rho\rangle+\tilde{\rho}\right]\nabla \cdot\left[\langle\mathbf{u}\rangle+\tilde{\mathbf{u}}\right] +
   \nabla\left[\langle\rho\rangle+\tilde{\rho}\right]\cdot \left[\langle\mathbf{u}\rangle+\tilde{\mathbf{u}}\right]= 0

If we neglect higher than first order terms in the fluctuations, the mass balance equation becomes

  \frac{\partial\tilde{\rho}}{\partial t} +
    \left[\langle\rho\rangle+\tilde{\rho}\right]\nabla \cdot\langle\mathbf{u}\rangle+
    \langle\rho\rangle\nabla\cdot\tilde{\mathbf{u}} +
    \nabla\langle\rho\rangle\cdot\tilde{\mathbf{u}}= 0

Assumption 2: Homogeneous medium[edit]

Next we assume that the medium is homogeneous, i.e.,

  \nabla\langle \rho \rangle = 0 ~.

Then the mass balance equation takes the form

  \frac{\partial\tilde{\rho}}{\partial t} +
    \left[\langle\rho\rangle+\tilde{\rho}\right]\nabla \cdot\langle\mathbf{u}\rangle+
    \langle\rho\rangle\nabla\cdot\tilde{\mathbf{u}} +
    = 0

Assumption 3: Medium at rest[edit]

At this stage we assume that the medium is at rest, i.e., \langle\mathbf{u}\rangle = 0 . Then the mass balance equation can be expressed as

  \frac{\partial\tilde{\rho}}{\partial t} +
    \langle\rho\rangle\nabla\cdot\tilde{\mathbf{u}} = 0

Assumption 4: Ideal gas, adiabatic, reversible[edit]

In order to close the system of equations we need an equation of state for the pressure. To do that we assume that the medium is an ideal gas and all acoustic waves compress the medium in an adiabatic and reversible manner. The equation of state can then be expressed in the form of the differential equation:

    \cfrac{dp}{d\rho} = \cfrac{\gamma~p}{\rho} ~;~~ \gamma := \cfrac{c_p}{c_v} ~;~~ c^2 = \cfrac{\gamma~p}{\rho} ~.

where c_p is the specific heat at constant pressure, c_v is the specific heat at constant volume, and c is the wave speed. The value of \gamma is 1.4 if the acoustic medium is air.

For small disturbances

   \cfrac{dp}{d\rho} \approx \cfrac{\tilde{p}}{\tilde{\rho}}  ~;~~
   \cfrac{p}{\rho} \approx \cfrac{\langle p \rangle}{\langle \rho \rangle} ~;~~
   c^2 \approx c_0^2 = \cfrac{\gamma~\langle p\rangle}{\langle \rho \rangle} ~.

where c_0 is the speed of sound in the medium.


   \cfrac{\tilde{p}}{\tilde{\rho}} = \gamma~\cfrac{\langle p \rangle}{\langle \rho \rangle}
     = c_0^2 \qquad \implies \qquad
   \cfrac{\partial\tilde{p}}{\partial t} = c_0^2 \cfrac{\partial\tilde{\rho}}{\partial t}

The balance of mass can then be written as

  \cfrac{1}{c_0^2}\frac{\partial\tilde{p}}{\partial t} +
    \langle\rho\rangle\nabla\cdot\tilde{\mathbf{u}} = 0

Dropping the tildes and defining \rho_0 :=  \langle\rho\rangle gives us the commonly used expression for the balance of mass in an acoustic medium:

  \frac{\partial p}{\partial t} + \rho_0~c_0^2~\nabla\cdot\mathbf{u} = 0 ~.

Governing equations in cylindrical coordinates[edit]

If we use a cylindrical coordinate system (r,\theta,z) with basis vectors \mathbf{e}_r, \mathbf{e}_\theta, \mathbf{e}_z, then the gradient of p and the divergence of \mathbf{u} are given by

     \nabla p & = \cfrac{\partial p}{\partial r}~\mathbf{e}_r + \cfrac{1}{r}~\cfrac{\partial p}{\partial \theta}~\mathbf{e}_\theta + \cfrac{\partial p}{\partial z}~\mathbf{e}_z \\
     \nabla\cdot\mathbf{u} & = \cfrac{\partial u_r}{\partial r} + \cfrac{1}{r}\left(\cfrac{\partial u_\theta}{\partial \theta} + u_r\right) + \cfrac{\partial u_z}{\partial z}

where the flow velocity has been expressed as \mathbf{u} = u_r~\mathbf{e}_r+u_\theta~\mathbf{e}_\theta+u_z~\mathbf{e}_z.

The equations for the conservation of momentum may then be written as

   \rho_0~\left[\cfrac{\partial u_r}{\partial t}~\mathbf{e}_r+\cfrac{\partial u_\theta}{\partial t}~\mathbf{e}_\theta+\cfrac{\partial u_z}{\partial t}~\mathbf{e}_z\right] +
\cfrac{\partial p}{\partial r}~\mathbf{e}_r + \cfrac{1}{r}~\cfrac{\partial p}{\partial \theta}~\mathbf{e}_\theta + \cfrac{\partial p}{\partial z}~\mathbf{e}_z = 0

In terms of components, these three equations for the conservation of momentum in cylindrical coordinates are

    \rho_0~\cfrac{\partial u_r}{\partial t} + \cfrac{\partial p}{\partial r}  = 0 ~;~~
    \rho_0~\cfrac{\partial u_\theta}{\partial t} + \cfrac{1}{r}~\cfrac{\partial p}{\partial \theta}  = 0 ~;~~
    \rho_0~\cfrac{\partial u_z}{\partial t} + \cfrac{\partial p}{\partial z}  = 0 ~.

The equation for the conservation of mass can similarly be written in cylindrical coordinates as

   \cfrac{\partial p}{\partial t} + \kappa\left[\cfrac{\partial u_r}{\partial r} + \cfrac{1}{r}\left(\cfrac{\partial u_\theta}{\partial \theta} + u_r\right) + \cfrac{\partial u_z}{\partial z}\right] = 0 ~.

Time harmonic acoustic equations in cylindrical coordinates[edit]

The acoustic equations for the conservation of momentum and the conservation of mass are often expressed in time harmonic form (at fixed frequency). In that case, the pressures and the flow velocity are assumed to be time harmonic functions of the form

   p(\mathbf{x}, t) = \hat{p}(\mathbf{x})~e^{-i\omega t} ~;~~
   \mathbf{u}(\mathbf{x}, t) = \hat{\mathbf{u}}(\mathbf{x})~e^{-i\omega t} ~;~~ i := \sqrt{-1}

where \omega is the frequency. Substitution of these expressions into the governing equations in cylindrical coordinates gives us the fixed frequency form of the conservation of momentum

   \cfrac{\partial\hat{p}}{\partial r} = i\omega~\rho_0~\hat{u}_r ~;~~
   \cfrac{1}{r}~\cfrac{\partial\hat{p}}{\partial \theta} = i\omega~\rho_0~\hat{u}_\theta ~;~~
    \cfrac{\partial\hat{p}}{\partial z} = i\omega~\rho_0~\hat{u}_z

and the fixed frequency form of the conservation of mass

   \cfrac{i\omega \hat{p}}{\kappa} = \cfrac{\partial \hat{u}_r}{\partial r} + \cfrac{1}{r}\left(\cfrac{\partial \hat{u}_\theta}{\partial \theta} + \hat{u}_r\right) + \cfrac{\partial \hat{u}_z}{\partial z} ~.

Special case: No z-dependence[edit]

In the special case where the field quantities are independent of the z-coordinate we can eliminate u_r, u_\theta to get

   \frac{\partial^2 p}{\partial r^2} + \frac{1}{r}\frac{\partial p}{\partial r} +
     \frac{1}{r^2}~\frac{\partial^2 p}{\partial\theta^2} + \frac{\omega^2\rho_0}{\kappa}~p = 0

Assuming that the solution of this equation can be written as

    p(r,\theta) = R(r)~Q(\theta)

we can write the partial differential equation as

   \cfrac{r^2}{R}~\cfrac{d^2R}{dr^2} + \cfrac{r}{R}~\cfrac{dR}{dr} + \cfrac{r^2\omega^2\rho_0}{\kappa} = -\cfrac{1}{Q}~\cfrac{d^2Q}{d\theta^2}

The left hand side is not a function of \theta while the right hand side is not a function of r. Hence,

   r^2~\cfrac{d^2R}{dr^2} + r~\cfrac{dR}{dr} + \cfrac{r^2\omega^2\rho_0}{\kappa}~R = \alpha^2~R ~;~~ \cfrac{d^2Q}{d\theta^2} = -\alpha^2~Q

where \alpha^2 is a constant. Using the substitution

   \tilde{r} \leftarrow \left(\omega\sqrt{\cfrac{\rho_0}{\kappa}}\right) r = k~r

we have

   \tilde{r}^2~\cfrac{d^2R}{d\tilde{r}^2} + \tilde{r}~\cfrac{dR}{d\tilde{r}} + (\tilde{r}^2-\alpha^2)~R = 0 ~;~~ \cfrac{d^2Q}{d\theta^2} = -\alpha^2~Q

The equation on the left is the Bessel equation which has the general solution

   R(r) =  A_\alpha~J_\alpha(k~r) + B_\alpha~J_{-\alpha}(k~r)

where J_\alpha is the cylindrical Bessel function of the first kind and A_\alpha, B_\alpha are undetermined constants. The equation on the right has the general solution

   Q(\theta) = C_\alpha~e^{i\alpha\theta} + D_\alpha~e^{-i\alpha\theta}

where C_\alpha,D_\alpha are undetermined constants. Then the solution of the acoustic wave equation is

   p(r,\theta) = \left[A_\alpha~J_\alpha(k~r) + B_\alpha~J_{-\alpha}(k~r)\right]\left(C_\alpha~e^{i\alpha\theta} + D_\alpha~e^{-i\alpha\theta}\right)

Boundary conditions are needed at this stage to determine \alpha and the other undetermined constants.


  1. ^ Douglas D. Reynolds. (1981). Engineering Principles in Acoustics, Allyn and Bacon Inc., Boston.

See also[edit]

Original courtesy of Wikipedia: http://en.wikipedia.org/wiki/Acoustic_theory — Please support Wikipedia.
This page uses Creative Commons Licensed content from Wikipedia. A portion of the proceeds from advertising on Digplanet goes to supporting Wikipedia.
1 videos found


https://youtube.com/devicesupport http://m.youtube.com

1 videos found

13 news items


Tue, 14 Apr 2015 12:02:46 -0700

You heard Theory of a Deadman rock before, but the band is about to give you a more intimate look at some of their songs with a brand new acoustic collection. The Angel Acoustic EP will be released on April 28 through Roadrunner/604. As the title would ...

Ottawa Citizen

Ottawa Citizen
Wed, 20 May 2015 07:08:35 -0700

Lockheed Martin Canada Mission Systems & Training (LM Canada MST) announced it awarded research and development funding to Ottawa-based Maritime Way Scientific Inc. The funding will support continued development of Maritime Way's advanced ...


Mon, 02 Feb 2015 08:01:41 -0800

Discover Beijing III: A Modern City In Motion highlights Beijing's transformation from an ancient city to a modern capital with high rises, with global stores and businesses established. This draws a contrast with the quiet life of the Hutong, which is ...


Tue, 27 May 2014 07:30:22 -0700

At Rothenberg, students studied sound wave and acoustic theory with the help of the Cincinnati Symphony Orchestra and Cincinnati Pops, whose conductor John Morris Russell in an acoustics Q&A with students. SCPA developed a program for offering ...
Thu, 16 Feb 2012 22:53:18 -0800

... the mid-summer sunrise and mid-winter sunset and there is widespread agreement that it was used for cremation burials. “However, I don't think you'll find many archaeologists who know about Stonehenge giving this particular acoustic theory a lot of ...
The Engineer
Tue, 26 Feb 2013 16:31:44 -0800

Those measurements, and the acoustic theory we developed to interpret the data, provided exactly the foundation we needed to undertake this critically important study that will be relevant to the seabed in somewhat deeper waters. 'As a greenhouse gas, ...
Fri, 07 Dec 2012 14:41:15 -0800

Warren wrote in 2005 that Marshall was developing an acoustic theory of lateral reflections for auditoria. The town hall became "his most successful guinea pig". Hamish Hay told his fellow councillors that the expense of Warren's trip was justified by ...
Stereophile Magazine
Sat, 08 Dec 2007 00:00:00 -0800

Audiophilia's underestimation of a room's effect on musical fidelity can be traced to two causes: 1) the science of acoustics is quite mathematical and arcane, and 2) acoustic theory has not been sufficiently popularized to translate into practical ...

Oops, we seem to be having trouble contacting Twitter

Support Wikipedia

A portion of the proceeds from advertising on Digplanet goes to supporting Wikipedia. Please add your support for Wikipedia!

Searchlight Group

Digplanet also receives support from Searchlight Group. Visit Searchlight