Acoustic theory: Difference between revisions

From Flowpedia
No edit summary
No edit summary
Line 9: Line 9:




\section{Acoustic Theory}
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.


\noindent 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.\\
<math display="block">
\dfrac{D\rho}{Dt}+\rho(\nabla\cdot\mathbf{v})=0
</math>


\[\dfrac{D\rho}{Dt}+\rho(\nabla\cdot\mathbf{v})=0\]
<math display="block">
\rho\dfrac{D\mathbf{v}}{Dt}+\nabla p=0
</math>


\[\rho\dfrac{D\mathbf{v}}{Dt}+\nabla p=0\]\\
Remember that <math>D/Dt</math> denotes the substantial derivative operator defined as follows


\noindent Remember that $D/Dt$ denotes the substantial derivative operator defined as follows\\
<math display="block">
\dfrac{D}{Dt}=\dfrac{\partial }{\partial t}+\mathbf{v}\cdot\nabla
</math>


\[\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.


\noindent where $\partial/\partial t$ is the local temporal derivative and $\mathbf{v}\cdot\nabla$ is the convective derivative.\\
We are going to analyze acoustic waves in one dimension, which means that the equations above reduces to


\noindent We are going to analyze acoustic waves in one dimension, which means that the equations above reduces to\\
<math display="block">
 
\begin{equation}
\dfrac{\partial \rho}{\partial t}+u\dfrac{\partial \rho}{\partial x}+\rho\dfrac{\partial u}{\partial x}=0
\dfrac{\partial \rho}{\partial t}+u\dfrac{\partial \rho}{\partial x}+\rho\dfrac{\partial u}{\partial x}=0
\label{eq:unstady:acoustic:wave:cont}
</math>
\end{equation}\\


\begin{equation}
<math display="block">
\rho\dfrac{\partial u}{\partial t}+\rho u\dfrac{\partial u}{\partial x}+\dfrac{\partial p}{\partial x}=0
\rho\dfrac{\partial u}{\partial t}+\rho u\dfrac{\partial u}{\partial x}+\dfrac{\partial p}{\partial x}=0
\label{eq:unstady:acoustic:wave:mom}
</math>
\end{equation}\\


\noindent 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 ($\rho$) and entropy ($s$).\\
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\]\\
<math display="block">
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
</math>


\noindent Since weak acoustic waves are considered, entropy will be constant and thus $ds=0$, which means that\\
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\]\\
<math display="block">
dp=\left(\dfrac{\partial p}{\partial \rho}\right)_s d\rho=a^2d\rho
</math>


\begin{equation}
<math display="block">
\rho\dfrac{\partial u}{\partial t}+\rho u\dfrac{\partial u}{\partial x}+a^2\dfrac{\partial \rho}{\partial x}=0
\rho\dfrac{\partial u}{\partial t}+\rho u\dfrac{\partial u}{\partial x}+a^2\dfrac{\partial \rho}{\partial x}=0
\label{eq:unstady:acoustic:wave:mom:b}
</math>
\end{equation}\\


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


\[
<math display="block">
\begin{aligned}
\begin{aligned}
&\rho=\rho_\infty+\Delta \rho\\
&\rho=\rho_\infty+\Delta \rho\\
Line 57: Line 62:
&u=u_\infty+\Delta u=\{u_\infty=0\}=\Delta u\\
&u=u_\infty+\Delta u=\{u_\infty=0\}=\Delta u\\
\end{aligned}
\end{aligned}
\]\\
</math>


\noindent 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\\
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


\begin{equation}
<math display="block">
\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
\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
\label{eq:unstady:acoustic:wave:cont:pert}
</math>
\end{equation}


\begin{equation}
<math display="block">
(\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
(\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
\label{eq:unstady:acoustic:wave:mom:pert}
</math>
\end{equation}\\


\noindent In the same way as pressure, being a thermodynamic variable, can be expressed as a function of two other thermodynamic varaibles, so can the speed of sound. Once again we will select density and entropy as the two thermodynamic variables\\
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)\]
<math display="block">
a^2=a^2(\rho,s)
</math>


\noindent and since entropy is constant\\
and since entropy is constant


\[a^2=a^2(\rho)\]\\
<math display="block">
a^2=a^2(\rho)
</math>


\noindent Taylor expansion of $a^2$ around the reference state $a_\infty$ with $\Delta \rho=\rho-\rho_\infty$ gives\\
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\]\\
<math display="block">
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
</math>


\noindent Inserted in Eqn.~\ref{eq:unstady:acoustic:wave:mom:pert}, we get\\
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\]\\
<math display="block">
(\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
</math>


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


\begin{equation}
<math display="block">
\dfrac{\partial}{\partial t}(\Delta \rho)+\rho_\infty\dfrac{\partial}{\partial x}(\Delta u)=0
\dfrac{\partial}{\partial t}(\Delta \rho)+\rho_\infty\dfrac{\partial}{\partial x}(\Delta u)=0
\label{eq:unstady:acoustic:wave:cont:linear}
</math>
\end{equation}


\begin{equation}
<math display="block">
\rho_\infty\dfrac{\partial}{\partial t}(\Delta u)+a_\infty^2\dfrac{\partial}{\partial x}(\Delta \rho)=0
\rho_\infty\dfrac{\partial}{\partial t}(\Delta u)+a_\infty^2\dfrac{\partial}{\partial x}(\Delta \rho)=0
\label{eq:unstady:acoustic:wave:mom:linear}
</math>
\end{equation}\\


\noindent 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.\\
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.


\noindent The temporal derivative of the continuity equation:\\
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\]\\
<math display="block">
\dfrac{\partial^2}{\partial t^2}(\Delta \rho)+\rho_\infty\dfrac{\partial^2}{\partial x\partial t}(\Delta u)=0
</math>


\noindent The divergence of the momentum equation:\\
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\]\\
<math display="block">
\rho_\infty\dfrac{\partial^2}{\partial x \partial t}(\Delta u)+a_\infty^2\dfrac{\partial^2}{\partial x^2}(\Delta \rho)=0
</math>


\noindent 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\\
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


\begin{equation}
<math display="block">
\frac{\partial^2}{\partial t^2}(\Delta \rho)=a^2_\infty\frac{\partial^2}{\partial x^2}(\Delta \rho)
\frac{\partial^2}{\partial t^2}(\Delta \rho)=a^2_\infty\frac{\partial^2}{\partial x^2}(\Delta \rho)
\label{eq:wave}
</math>
\end{equation}\\


<!--
%Combining the linearized continuity and momentum equations we get the wave propagation equation for acoustic waves (small perturbations)\\
%Combining the linearized continuity and momentum equations we get the wave propagation equation for acoustic waves (small perturbations)\\
-->
which is a one-dimensional form of the classic wave equation with the general solution


 
<math display="block">
 
\noindent which is a one-dimensional form of the classic wave equation with the general solution\\
 
\begin{equation}
\Delta \rho = F(x-a_\infty t)+G(x+a_\infty t)
\Delta \rho = F(x-a_\infty t)+G(x+a_\infty t)
\label{eq:wave:solution}
</math>
\end{equation}\\


\noindent $F$ and $G$ are arbitrary functions. The function $F$ describes the shape of a wave traveling in the positive $x$-direction at the speed of sound of the ambient gas and the function $G$ describes the shape of a wave traveling in the negative $x$-direction at the same speed. In Eqn.~\ref{eq:wave} $\Delta \rho$ 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}.\\
<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}.


\begin{equation*}
<math display="block">
\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)=\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}
\end{equation*}\\
</math>


\begin{equation}
<math display="block">
\frac{\partial}{\partial t}(\Delta \rho)=-a_\infty F'+a_\infty G'
\frac{\partial}{\partial t}(\Delta \rho)=-a_\infty F'+a_\infty G'
\end{equation}\\
</math>


\begin{equation*}
<math display="block">
\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\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}
\end{equation*}\\
</math>


\begin{equation}
<math display="block">
\frac{\partial^2}{\partial t^2}(\Delta \rho)=a_\infty^2 F''+a_\infty^2 G''
\frac{\partial^2}{\partial t^2}(\Delta \rho)=a_\infty^2 F''+a_\infty^2 G''
\label{eq:wave:ddt}
</math>
\end{equation}\\


\begin{equation*}
<math display="block">
\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)=\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}
\end{equation*}\\
</math>


\begin{equation}
<math display="block">
\frac{\partial}{\partial x}(\Delta \rho)=F'+G'
\frac{\partial}{\partial x}(\Delta \rho)=F'+G'
\end{equation}\\
</math>


\begin{equation*}
<math display="block">
\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)=\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}
\end{equation*}\\
</math>


\begin{equation}
<math display="block">
\frac{\partial^2}{\partial x^2}(\Delta \rho)= F'' + G''
\frac{\partial^2}{\partial x^2}(\Delta \rho)= F'' + G''
\label{eq:wave:ddx}
</math>
\end{equation}\\


\noindent Eqns. \ref{eq:wave:ddt} and \ref{eq:wave:ddx} inserted Eqn. \ref{eq:wave} gives\\
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'')\]\\
<math display="block">
a_\infty^2 F''+a_\infty^2 G''=a_\infty^2(F'' + G'')
</math>


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


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


\begin{equation}
<math display="block">
\Delta \rho(x,t)=F(x-a_\infty t)
\Delta \rho(x,t)=F(x-a_\infty t)
\label{eq:wave:solution:F}
</math>
\end{equation}\\


\noindent If $\Delta \rho$ is constant, i.e. a wave with constant amplitude, we see from Eqn.~\ref{eq:wave:solution:F} that $(x-a_\infty t)$ is constant and thus\\
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\]\\
<math display="block">
x=a_\infty t + c\Rightarrow \dfrac{dx}{dt}=a_\infty
</math>


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


\[\dfrac{\partial}{\partial t}(\Delta \rho)=-a_\infty F'\]\\
<math display="block">
\dfrac{\partial}{\partial t}(\Delta \rho)=-a_\infty F'
</math>


\[\dfrac{\partial}{\partial x}(\Delta \rho)=F'\]\\
<math display="block">
\dfrac{\partial}{\partial x}(\Delta \rho)=F'
</math>


\noindent and thus\\
and thus


\begin{equation}
<math display="block">
\dfrac{\partial}{\partial x}(\Delta \rho)=-\dfrac{1}{a_\infty}\dfrac{\partial}{\partial t}(\Delta \rho)
\dfrac{\partial}{\partial x}(\Delta \rho)=-\dfrac{1}{a_\infty}\dfrac{\partial}{\partial t}(\Delta \rho)
\label{eq:wave:solution:F:b}
</math>
\end{equation}\\


\noindent which gives a relation between the temporal derivative of $\Delta \rho$ and the spatial derivative of $\Delta \rho$. 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\\
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\]
<math display="block">
\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
</math>


\[\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\]\\
<math display="block">
\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
</math>


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


\[\Delta u-\dfrac{a_\infty}{\rho_\infty}\Delta \rho=0\]\\
<math display="block">
\Delta u-\dfrac{a_\infty}{\rho_\infty}\Delta \rho=0
</math>


\noindent or\\
or


\[\Delta u=\dfrac{a_\infty}{\rho_\infty}\Delta \rho\]\\
<math display="block">
\Delta u=\dfrac{a_\infty}{\rho_\infty}\Delta \rho
</math>


\noindent If instead $F$ is set to zero and $G$ is non-zero, we get\\
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\]\\
<math display="block">
\Delta u=-\dfrac{a_\infty}{\rho_\infty}\Delta \rho
</math>


\[\left(\dfrac{\partial p}{\partial \rho}\right)_s=a^2\Rightarrow \Delta p=a_\infty^2 \Delta \rho\]\\
<math display="block">
\left(\dfrac{\partial p}{\partial \rho}\right)_s=a^2\Rightarrow \Delta p=a_\infty^2 \Delta \rho
</math>


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


\begin{equation}
<math display="block">
\Delta u=\dfrac{a_\infty}{\rho_\infty}\Delta \rho=\dfrac{1}{a_\infty \rho_\infty}\Delta p
\Delta u=\dfrac{a_\infty}{\rho_\infty}\Delta \rho=\dfrac{1}{a_\infty \rho_\infty}\Delta p
\label{eq:wave:solution:F:c}
</math>
\end{equation}\\


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


\begin{equation}
<math display="block">
\Delta u=-\dfrac{a_\infty}{\rho_\infty}\Delta \rho=-\dfrac{1}{a_\infty \rho_\infty}\Delta p
\Delta u=-\dfrac{a_\infty}{\rho_\infty}\Delta \rho=-\dfrac{1}{a_\infty \rho_\infty}\Delta p
\label{eq:wave:solution:G:c}
</math>
\end{equation}\\