Acoustic theory: Difference between revisions

In the following we are going to derive the linear acoustic wave equation starting from the continuity and momentum equations on non-conservation differential form. The equations are repeated here for convenience.

\dfrac{D\rho}{Dt}+\rho(\nabla\cdot\mathbf{v})=0

\rho\dfrac{D\mathbf{v}}{Dt}+\nabla p=0

Remember that <math>D/Dt</math> denotes the substantial derivative operator defined as follows

\dfrac{D}{Dt}=\dfrac{\partial }{\partial t}+\mathbf{v}\cdot\nabla

where <math>\partial/\partial t</math> is the local temporal derivative and <math>\mathbf{v}\cdot\nabla</math> is the convective derivative.

We are going to analyze acoustic waves in one dimension, which means that the equations above reduces to

\dfrac{\partial \rho}{\partial t}+u\dfrac{\partial \rho}{\partial x}+\rho\dfrac{\partial u}{\partial x}=0

\rho\dfrac{\partial u}{\partial t}+\rho u\dfrac{\partial u}{\partial x}+\dfrac{\partial p}{\partial x}=0

Pressure is a thermodynamic property and thus it can be expressed as a function of two other thermodynamic properties. Let's express pressure as a function of density (<math>\rho</math>) and entropy (<math>s</math>).

p=p(\rho,s)\Rightarrow dp=\left(\dfrac{\partial p}{\partial \rho}\right)_s d\rho+\left(\dfrac{\partial p}{\partial s}\right)_\rho ds

Since weak acoustic waves are considered, entropy will be constant and thus <math>ds=0</math>, which means that

dp=\left(\dfrac{\partial p}{\partial \rho}\right)_s d\rho=a^2d\rho

\rho\dfrac{\partial u}{\partial t}+\rho u\dfrac{\partial u}{\partial x}+a^2\dfrac{\partial \rho}{\partial x}=0

The acoustic perturbations can be described as small deviations around a reference state

\begin{aligned}

&\rho=\rho_\infty+\Delta \rho\\

&u=u_\infty+\Delta u=\{u_\infty=0\}=\Delta u\\

\end{aligned}

Inserted in Eqns.~\ref{eq:unstady:acoustic:wave:cont} and \ref{eq:unstady:acoustic:wave:mom:b} and using the fact that derivatives of the constant reference state flow quantities are zero, we get

\dfrac{\partial}{\partial t}(\Delta \rho)+\Delta u\dfrac{\partial}{\partial x}(\Delta \rho)+(\rho_\infty+\Delta \rho)\dfrac{\partial}{\partial x}(\Delta u)=0

(\rho_\infty+\Delta \rho)\dfrac{\partial}{\partial t}(\Delta u)+(\rho_\infty+\Delta \rho) \Delta u\dfrac{\partial}{\partial x}(\Delta u)+a^2\dfrac{\partial}{\partial x}(\Delta \rho)=0

In the same way as pressure, being a thermodynamic variable, can be expressed as a function of two other thermodynamic variables, so can the speed of sound. Once again we will select density and entropy as the two thermodynamic variables

a^2=a^2(\rho,s)

and since entropy is constant

a^2=a^2(\rho)

Taylor expansion of <math>a^2</math> around the reference state <math>a_\infty</math> with <math>\Delta \rho=\rho-\rho_\infty</math> gives

a^2=a^2_\infty+\left(\dfrac{\partial}{\partial \rho}(a^2)\right)_\infty\Delta \rho+\left(\dfrac{\partial^2}{\partial \rho^2}(a^2)\right)_\infty(\Delta \rho)^2+\ \cdots

Inserted in Eqn.~\ref{eq:unstady:acoustic:wave:mom:pert}, we get

(\rho_\infty+\Delta \rho)\dfrac{\partial}{\partial t}(\Delta u)+(\rho_\infty+\Delta \rho) \Delta u\dfrac{\partial}{\partial x}(\Delta u)+\left[a^2_\infty+\left(\dfrac{\partial}{\partial \rho}(a^2)\right)_\infty\Delta \rho+\ \cdots\right]\dfrac{\partial}{\partial x}(\Delta \rho)=0

The perturbations <math>\Delta u</math> and <math>\Delta \rho</math> are small, which implies that <math>\Delta u \ll a_\infty</math> and <math>\Delta \rho \ll \rho_\infty</math>. This means that products of perturbations can be canceled and so can higher-order terms in the Taylor expansion of <math>a^2</math>. This means that the continuity and momentum equations reduces to

\dfrac{\partial}{\partial t}(\Delta \rho)+\rho_\infty\dfrac{\partial}{\partial x}(\Delta u)=0

\rho_\infty\dfrac{\partial}{\partial t}(\Delta u)+a_\infty^2\dfrac{\partial}{\partial x}(\Delta \rho)=0

Before making the assumption that the perturbations are small compared to the corresponding reference state flow quantities and thus justifying the cancelation of products of perturbations from the equations, the flow equations were still the exact fully non-linear equations. Eqns.~\ref{eq:unstady:acoustic:wave:cont:linear}. and \ref{eq:unstady:acoustic:wave:mom:linear}, however, are approximations as several terms has been removed. The equations are linear and are good approximations as long as the perturbations are small. The smaller the perturbations, the better the approximation are the linear equations. Eqns.~\ref{eq:unstady:acoustic:wave:cont:linear} and \ref{eq:unstady:acoustic:wave:mom:linear} describes the motion induced in a gas by the passage of a sound wave. By combining the temporal derivative of Eqn.~\ref{eq:unstady:acoustic:wave:cont:linear} with the divergence of Eqn.~\ref{eq:unstady:acoustic:wave:mom:linear}, it is possible to obtain a wave equation describing the propagation of acoustic waves in a quiescent surrounding.

The temporal derivative of the continuity equation:

\dfrac{\partial^2}{\partial t^2}(\Delta \rho)+\rho_\infty\dfrac{\partial^2}{\partial x\partial t}(\Delta u)=0

The divergence of the momentum equation:

\rho_\infty\dfrac{\partial^2}{\partial x \partial t}(\Delta u)+a_\infty^2\dfrac{\partial^2}{\partial x^2}(\Delta \rho)=0

The second term in the first equation is the same as the first term in the second equation. Substituting the term, the two equations reduces to one single equation

\frac{\partial^2}{\partial t^2}(\Delta \rho)=a^2_\infty\frac{\partial^2}{\partial x^2}(\Delta \rho)

<!--

which is a one-dimensional form of the classic wave equation with the general solution

\Delta \rho = F(x-a_\infty t)+G(x+a_\infty t)

<math>F</math> and <math>G</math> are arbitrary functions. The function <math>F</math> describes the shape of a wave traveling in the positive <math>x</math>-direction at the speed of sound of the ambient gas and the function <math>G</math> describes the shape of a wave traveling in the negative <math>x</math>-direction at the same speed. In Eqn.~\ref{eq:wave} <math>\Delta \rho</math> appears with second derivatives in space and time. Let's differentiate the proposed solution (Eqn.~\ref{eq:wave:solution}) two times in time and space, respectively, and check that it is actually a valid solution to Eqn.~\ref{eq:wave}.

\frac{\partial}{\partial t}(\Delta \rho)=\frac{\partial F}{\partial (x-a_\infty t)}\frac{\partial (x-a_\infty t)}{\partial t}+\frac{\partial G}{\partial (x+a_\infty t)}\frac{\partial (x+a_\infty t)}{\partial t}

\frac{\partial}{\partial t}(\Delta \rho)=-a_\infty F'+a_\infty G'

\frac{\partial^2}{\partial t^2}(\Delta \rho)=-a_\infty\frac{\partial F'}{\partial (x-a_\infty t)}\frac{\partial (x-a_\infty t)}{\partial t}+a_\infty\frac{\partial G'}{\partial (x+a_\infty t)}\frac{\partial (x+a_\infty t)}{\partial t}

\frac{\partial^2}{\partial t^2}(\Delta \rho)=a_\infty^2 F''+a_\infty^2 G''

\frac{\partial}{\partial x}(\Delta \rho)=\frac{\partial F}{\partial (x-a_\infty t)}\frac{\partial (x-a_\infty t)}{\partial x}+\frac{\partial G}{\partial (x+a_\infty t)}\frac{\partial (x+a_\infty t)}{\partial x}

\frac{\partial}{\partial x}(\Delta \rho)=F'+G'

\frac{\partial^2}{\partial x^2}(\Delta \rho)=\frac{\partial F'}{\partial (x-a_\infty t)}\frac{\partial (x-a_\infty t)}{\partial x}+\frac{\partial G'}{\partial (x+a_\infty t)}\frac{\partial (x+a_\infty t)}{\partial x}

\frac{\partial^2}{\partial x^2}(\Delta \rho)= F'' + G''

Eqns. \ref{eq:wave:ddt} and \ref{eq:wave:ddx} inserted Eqn. \ref{eq:wave} gives

a_\infty^2 F''+a_\infty^2 G''=a_\infty^2(F'' + G'')

which shows that Eqn. \ref{eq:wave:solution} is a valid solution to the wave equation.

<math>F</math> and <math>G</math> are arbitrary functions and thus <math>G=0</math> is a valid solution, which gives

\Delta \rho(x,t)=F(x-a_\infty t)

If <math>\Delta \rho</math> is constant, i.e. a wave with constant amplitude, we see from Eqn.~\ref{eq:wave:solution:F} that <math>(x-a_\infty t)</math> is constant and thus

x=a_\infty t + c\Rightarrow \dfrac{dx}{dt}=a_\infty

From Eqn.~\ref{eq:wave:solution:F}, we get

\dfrac{\partial}{\partial t}(\Delta \rho)=-a_\infty F'

\dfrac{\partial}{\partial x}(\Delta \rho)=F'

and thus

\dfrac{\partial}{\partial x}(\Delta \rho)=-\dfrac{1}{a_\infty}\dfrac{\partial}{\partial t}(\Delta \rho)

which gives a relation between the temporal derivative of <math>\Delta \rho</math> and the spatial derivative of <math>\Delta \rho</math>. With Eqn.~\ref{eq:wave:solution:F:b}, the linearized momentum equation Eqn.~\ref{eq:unstady:acoustic:wave:mom:linear} can be rewritten as follows

\dfrac{\partial}{\partial t}(\Delta u)=-\dfrac{a_\infty^2}{\rho_\infty}\dfrac{\partial}{\partial x}(\Delta \rho)=\left\{\dfrac{\partial}{\partial x}(\Delta \rho)=-\dfrac{1}{a_\infty}\dfrac{\partial}{\partial t}(\Delta \rho)\right\}=\dfrac{a_\infty}{\rho_\infty}\dfrac{\partial}{\partial t}(\Delta \rho)\Rightarrow

\dfrac{\partial}{\partial t}\left(\Delta u-\dfrac{a_\infty}{\rho_\infty}\Delta \rho\right)=0\Rightarrow \Delta u-\dfrac{a_\infty}{\rho_\infty}\Delta \rho = const

In an undisturbed gas <math>\Delta u=\Delta \rho=0</math> and thus

\Delta u-\dfrac{a_\infty}{\rho_\infty}\Delta \rho=0

or

\Delta u=\dfrac{a_\infty}{\rho_\infty}\Delta \rho

If instead <math>F</math> is set to zero and <math>G</math> is non-zero, we get

\Delta u=-\dfrac{a_\infty}{\rho_\infty}\Delta \rho

\left(\dfrac{\partial p}{\partial \rho}\right)_s=a^2\Rightarrow \Delta p=a_\infty^2 \Delta \rho

Acoustic wave traveling in the positive <math>x</math>-direction:

\Delta u=\dfrac{a_\infty}{\rho_\infty}\Delta \rho=\dfrac{1}{a_\infty \rho_\infty}\Delta p

Acoustic wave traveling in the negative <math>x</math>-direction:

\Delta u=-\dfrac{a_\infty}{\rho_\infty}\Delta \rho=-\dfrac{1}{a_\infty \rho_\infty}\Delta p