# Nonlinear acoustics

Nonlinear acoustics is a branch of physics and acoustics dealing with sound waves of sufficiently large amplitudes. Large amplitudes require using full systems of governing equations of fluid dynamics (for sound waves in liquids and gases) and elasticity (for sound waves in solids). These equations are generally nonlinear, and their traditional linearization is no longer possible. The solutions of these equations show that, due to the effects of nonlinearity, sound waves are being distorted as they travel.

## Introduction

A sound wave propagates through a material as a localized pressure change. Increasing the pressure of a gas, or fluid, increases its local temperature and also increases the local speed of sound in a compressible material increases with temperature; as a result, the wave travels faster during the high pressure phase of the oscillation than during the lower pressure phase. This affects the wave's frequency structure; for example, in an initially plane sinusoidal wave of a single frequency, the peaks of the wave travel faster than the troughs, and the pulse becomes cumulatively more like a sawtooth wave. In other words, the wave self-distorts. In doing so, other frequency components are introduced, which can be described by the Fourier series. This phenomenon is characteristic of a non-linear system, since a linear acoustic system responds only to the driving frequency. This always occurs but the affects of geometric spreading and of absorption usually overcome the self distortion, so linear behavior usually prevails and nonlinear acoustic propagation occurs only for very large amplitudes and only near the source.

Additionally, waves of different amplitudes will generate different pressure gradients, contributing to the non-linear effect.

## Physical analysis

The pressure changes within a medium cause the wave energy to transfer to higher harmonics. Since attenuation generally increases with frequency, a counter effect exists that changes the nature of the nonlinear effect over distance. To describe their level of nonlinearity, materials can be given a nonlinearity parameter, $B/A$. The values of $A$ and $B$ are the coefficients of the first and second order terms of the Taylor series expansion of the equation relating the material's pressure to its density. The Taylor series has more terms, and hence more coefficients (C, D, .. etc) but they are seldom used. Typical values for the nonlinearity parameter in biological mediums are shown in the following table.[1]

Material $B/A$
Blood 6.1
Brain 6.6
Fat 10
Liver 6.8
Muscle 7.4
Water 5.2

In a liquid usually a modified coefficient is used known as $\beta = 1 + \frac{B}{2A}$.

## Mathematical model

### Governing Equations to Derive Westervelt Equation

Continuity:

$\frac{\partial \rho}{\partial t} + \nabla (\rho \textbf{u}) = 0$

Conservation of momentum:

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

Equation of state with $\delta \rho = \rho - \rho_0$:

$p= c_0^2 \delta \rho \left( 1+ \frac{B}{2!A}\frac{\delta \rho}{\rho_0} + ... \right)$

If the second term in the Taylor expansion of pressure is dropped, the viscous wave equation can be derived. If it is kept, the non-linear term in pressure appears in the Westervelt equation.

### Westervelt equation

The general wave equation that accounts for nonlinearity up to the second-order is given by the Westervelt equation[2]

$\, \nabla^{2} p - \frac{1}{c_{0}^{2}} \frac{\partial^{2} p}{\partial t^{2}} + \frac{\delta}{c_{0}^{4}} \frac{\partial^{3} p}{\partial t^{3}} = - \frac{\beta}{\rho_{0} c_{0}^{4}} \frac{\partial^{2} p^{2}}{\partial t^{2}}$

where $p$ is the sound pressure, $c_0$ is the small signal sound speed, $\delta$ is the sound diffusivity, $\beta$ is the non-linearity coefficient and $\rho_0$ is the ambient density.

The sound diffusivity is given by

$\, \delta = \frac{1}{\rho_{0}} \left(\frac{4}{3}\mu+\mu_{B}\right) + \frac{k}{\rho_{0}} \left(\frac{1}{c_{v}} - \frac{1}{c_{p}}\right)$

where $\mu$ is the shear viscosity, $\mu_{B}$ the bulk viscosity, $k$ the thermal conductivity, $c_{v}$ and $c_{p}$ the specific heat at constant volume and pressure respectively.

### Burgers' equation

The Westervelt equation can be simplified to take a one-dimensional form with an assumption of strictly forward propagating waves and the use of a coordinate transformation to a retarded time frame. This is known as the Burgers' equation[3]

$\frac{\partial p}{\partial z} = \frac{\beta p}{\rho_{0} c_{0}^{3}}\frac{\partial p}{\partial \tau} + \frac{\delta}{2 c_{0}^{3}}\frac{\partial^{2} p}{\partial \tau^{2}}$

where $\tau = t-z/c_0$ is retarded time. The Burgers equation is the simplest model that describes the combined effects of nonlinearity and losses on the propagation of plane progressive waves.

### KZK equation

An augmentation to the Burgers equation that accounts for the combined effects of non-linearity, diffraction and absorption in directional sound beams is described by the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation.[4] Solutions to this equation are generally used to model non-linear acoustics.

If the $z$ axis is in the direction of the sound beam path and the $(x,y)$ plane is perpendicular to that, the KZK equation can be written[5]

$\, \frac{\partial^2 p}{\partial z \partial \tau} = \frac{c_0}{2}\nabla^2_{\perp}p + \frac{\delta}{2c^3_0}\frac{\partial^3 p}{\partial \tau^3} + \frac{\beta}{2\rho_0 c^3_0}\frac{\partial^2 p^2}{\partial \tau^2}$

The equation can be solved for a particular system using a finite difference scheme. Such solutions show how the sound beam distorts as it passes through a non-linear medium.

## Common occurrences

### Sonic boom

The nonlinear behavior of the atmosphere leads to change of the wave shape in a sonic boom. Generally, this makes the boom more 'sharp' or sudden, as the high-amplitude peak moves to the wavefront.

### Acoustic levitation

The practice of acoustic levitation would not be possible without understanding nonlinear acoustic phenomena.[6] The nonlinear effects are particularly evident due to the high-powered acoustic waves involved.

### Ultrasonic waves

Because of their relatively high amplitude to wavelength ratio, ultrasonic waves commonly display nonlinear propagation behavior. For example, nonlinear acoustics is a field of interest for medical ultrasonography because it can be exploited to produce a better image quality.

### Musical acoustics

The physical behavior of musical acoustics in mainly nonlinear. Many attempts are made to model their sound generation from physical modeling of emulating their sound from measurements of their non-linearity.[7]

## References

1. ^ Wells, P. N. T. (1999). "Ultrasonic imaging of the human body". Reports on Progress in Physics 62 (5): 671. doi:10.1088/0034-4885/62/5/201. edit
2. ^ Hamilton, M.F.; Blackstock, D.T. (1998). Nonlinear Acoustics. Academic Press. p. 55. ISBN 0-12-321860-8.
3. ^ Hamilton, M.F.; Blackstock, D.T. (1998). Nonlinear Acoustics. Academic Press. p. 57. ISBN 0-12-321860-8.
4. ^ Anna Rozanova-Pierrat. Mathematical analysis of Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation (PDF). Laboratoire Jacques-Louis Lions, Université Pierre et Marie Curie. Retrieved 2008-11-10.
5. ^ V. F. Humphrey. Non-linear propagation for medical imaging (PDF). Department of Physics, University of Bath, Bath, UK. Retrieved 2008-11-10.
6. ^ http://science.howstuffworks.com/acoustic-levitation.htm
7. ^ The Emulation of Nonlinear Time-Invariant Audio Systems with Memory by Means of Volterra Series, JAES Volume 60 Issue 12 pp. 984-996; December 2012.