Sound waves with attenuation

Just a simple derivation of the role of attenuation in the standard sound wave equation. Original work: Stokes, 1845.

Starting with the Navier-Stokes momentum equation

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

where \lambda is a Lamé viscosity coefficient. The bulk viscosity coeficient  is defined as  \zeta = \lambda + (2/3) \mu. The last term  is often neglected, even in compressible flow, but sound attenuation is one of the few cases where it may have some influence. All viscosities are assumed to be constant, but in this case this is a safe assumption, since we are going to assume small departures about equilibrium values.


Let us not forget mass continuity:
\nabla \cdot \mathbf{\rho u} + \frac{\partial \rho}{\partial t} =0

Now, let us suppose small pressure and density fluctuations about a constant background. I.e:
\mathbf{u} \approx \mathbf{u} \qquad p\approx p_0+\delta p \qquad \rho\approx \rho_0+\delta \rho .

The first equation means the medium onto which sounds travel is not moving (should be modified for e.g. wind). In this case, the fluctuations in pressure and density, being small, should be related by a constant which, with some foresight, we will call c^2:
\delta p \approx c^2 \delta \rho .

The viscosities were already assumed to have small variations about their equilibrium values. Then, neglecting second and higher order perturbation terms we obtain linearized equations, involving only the velocity and the pressure (not the density):

\frac{\partial }{\partial t} \mathbf{u} = - \frac{1}{\rho_0} \nabla p + \nu \nabla^2 \mathbf{u} + \left(\frac{\lambda}{\rho_0}+\nu\right) \nabla (\nabla\cdot\mathbf{u})  

\rho_0 c^2 \nabla \cdot \mathbf{\rho u} + \frac{\partial p}{\partial t} =0

Here, we have dropped the “\deltas” (it’s confusing, but the expressions are so much cleaner). The kinematic viscosity is defined as \nu=\mu/\rho_0.

If we differentiate with respect to time in the first, and with respect to space (applying \nabla) in the second, we can eliminate the pressure term \nabla(\partial p/\partial t) (I know, exchanging both derivatives, which makes mathematicians nervous.) The following equation for the velocity results:

\frac{\partial^2 }{\partial t^2 } \mathbf{u} =c^2\nabla (\nabla\cdot\mathbf{u})  + \nu  \nabla^2\frac{\partial\mathbf{u}}{\partial t}  + \left(\frac{\lambda}{\rho_0}+\nu\right) \nabla (\nabla\cdot \frac{\partial\mathbf{u}}{\partial t}   )  .

Clearly, just a wave equation if there was no viscosity. Let’s try a wave solution of the form

\mathbf{u} = \mathbf{u}_0 e^{i(\omega t - \beta x)} .

Here, \omega=2\pi f is the angular frequency of the sound wave, but  \beta may be a complex number. If we find (as we will)

 \beta =  k - \alpha i ,

that would mean

\mathbf{u} = \mathbf{u}_0e^{- \alpha x} e^{i(\omega t - k x)} ,

which clearly identifies \omega=2\pi /\lambda as the (real) wave-number for a wave-lenght \lambda, and \alpha=1 /L , as a sound attenuation coefficient, with L the penetration length.

Now, second time derivatives just yield

\frac{\partial^2 }{\partial t^2 } \mathbf{u} = -\alpha^2 \mathbf{u}

But, it is quite interesting that these two second space derivatives are not quite the same:

\nabla^2 \mathbf{u} =  -k^2 \mathbf{u}

\nabla (\nabla\cdot \mathbf{u} ) = -k^2 (\mathbf{u}\cdot\mathbf{i})\mathbf{i} ,

where \mathbf{i} is the unit vector in the x direction. Clearly, the second version of the derivative produces a longitudinal wave, with a vector component just along the direction of propagation!

Transverse waves

I know, sound waves are longitudinal. But, what happens if we plug these derivatives for the y Cartesian component. Well:

-\omega^2 = - i \nu \omega \beta^2,


\beta^2 = \frac{\omega}{i \nu} = - \frac{\omega i}{\nu} .

Now, this is a very classic complex analysis problem. Recall  -i=e^{3\pi/2+2\pi}  You do need the 2\pi to get the two solutions. On of the solutions has the signs reversed, and corresponds to a wave propagating and attenuating in the –x direction. The one we are looking for is:

\beta =\sqrt{ \frac{\omega}{ \nu}}e^{ 7 \pi/4}.


k=\alpha =\sqrt{ \frac{\omega}{ 2 \nu}}.

This is an interesting wave, since the attenuation coefficient is equal to the wave number. It looks like a simple function: e^{-x}\sin(x), which would be called simplistic by a physicist. If you plot it you’ll see it decays very fast, with just one maximum of minimum of importance. (You can just type “e^-x*cos(x) from 0 to 10” in google, it will plot it!). So, yes, sound waves are basically longitudinal, since their transverse components get attenuated very fast. How fast? As we said before, the attenuation length is the inverse of \alpha, so L=\sqrt{2\nu/\omega}. This means that for everyday “sound”, i.e. audible frequencies, that length is quite small. With the numerical data in the Table, that length will be only be about 0.7 mm for a very low frequency of 10 Hz (just below the hearing range), and will decrease as the inverse of the square root of the frequency.

c (m/s) \nu (m{}^2/s) \zeta/\rho_0 (m{}^2/s) f_c (GHz)
water 1480 1.0 \times 10^{-6} 3.1 \times 10^{-6} 79
air 340 1.48 \times 10^{-5} 0.5 \times 10^{-5} ? 0.7

Table: Numerical values for two important substances. Question marks are speculative, since in Cramer air is said to have a negligible bulk viscosity… but in a graph this quantity’s value is seen to be about half the shear viscosity value, and air is mostly nitrogen.

Longitudinal waves

Plugging the derivatives for the x Cartesian component:

-\omega^2 = - c^2\beta^2 - i (4\nu/3+\zeta/\rho_0)  \omega \beta^2.

We see that the viscosity now features the bulk viscosity, as seems fitting for a longitudinal disturbance, which involves compression. Also, a new important term appears. If viscosity were negligible, the solution is just

 \omega =  c \beta ,

the usual dispersion relation for sound. In general, though:

\beta^2 = \frac{\omega^2}{c^2+i(4\nu/3+\zeta/\rho_0)\omega}=\frac{\omega^2}{c^2}\frac{1}{1+ i\omega/\omega_c},

where we define the important crossover angular frequency

\omega_c := \frac{c^2}{4\nu/3+\zeta/\rho_0}.

Numerical values can be found in the table (we provide the linear frequency). The value is really high for hearing range, which goes up to about 20 kHz for humans, 160kHz for porpoises, which hold the record. Medical ultrasound goes as high as 16 MHz only, and only acoustic microscopy reaches a few GHz, the range at which our value for air sits.

Low frequencies

At frequencies much below the crossover frequency, we may expand the term in the denominator in a Taylor series, then again for the square root. The end result is

\beta = \frac{\omega}{c} (1-i   \omega/(2 \omega_c)).


Therefore the wave number is

\beta = \frac{\omega}{c},

as if there were no viscosity. The attenuation coefficient is

\alpha = \frac{\omega^2}{2 c\omega_c}=\frac{\omega^2}{2 c^3}(4\nu/3+\zeta/\rho_0).

It therefore grows as the square of the frequency.  This agrees with the expression in wikipedia (but for a factor of two in the general formula, which I have just corrected – let’s see if it stays). This would mean that for that very high porpoise’s pitch at 160Hz the attenuation length would be about 1.5 kilometers (in water, of course), which may be important for long-range communication of these animals. For medical ultrasound at 16MHz this length is 14cm, which can clearly have an impact for human tissues (I am not sure if this attenuation may be used to our advantage). For the human bodies I have used water values, which is a fair approximation.

High frequencies

At frequencies much higher than the crossover, we may neglect the “1” in the denominator, to obtain

\beta^2 = \frac{\omega\omega_c}{i c^2} .

Now, this looks familiar, especially if we substitute the crossover frequency:

\beta^2 = \frac{\omega}{i(4\nu/3+\zeta/\rho_0)} .

Basically, the same expression as for transverse waves, with a different combination of viscosities. The conclusion is similar:

k=\alpha =\sqrt{ \frac{\omega}{ 2(4\nu/3+\zeta/\rho_0)}},

and the resulting waves are heavily attenuated.

The crossover

A fast check on the above approximations is to see what happens when the frequency is exactly the crossover frequency. At this point, the growth of the attenuation coefficient as the square of frequency crosses over to a growth as the square root. This would mean a bend in a log-log plot, between two straight lines with different slopes.

The extrapolation of the low frequency expression yields

\alpha_c\approx \omega_c/(2 c) ,

whereas the high frequency expression yields

\alpha_c\approx \omega_c/(\sqrt{2} c) ,

not so different at all.

The exact expression can be shown to be, after some complex algebra,

\alpha_c\approx \omega_c \sin(\pi/8) / (2^{1/4} c) ,

a value just a bit below the other two. This means the approximations remain quite fair up to the limit of their respective ranges.

By the way, for water this value corresponds to an attenuation length of about 58 nm, which is really short. For air, it is about 1.5 micrometers, the size of a really small cell.




Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s