Friday, 22 March 2019

fluid dynamics - assumptions about sound waves


When deriving the sound wave equation: 1c22pt2=Δ2p by linearizing the Euler equation:


ρdvdt=p and the continuity equation:


ρt+(ρv)=0 using an approach of small deviations ρ,v,pfrom an equilibrium ρ0,v0,p0 with v0=0. So p=p0+pv=0+vρ=ρ0+ρ


Why can we just neglect the convective term in the euler equation? Meaning why can we use: ρ0vt=p?


Why can we just assume (v)v0 ?



Answer



The assumption when linearizing is that the deviations/perturbations are very small compared to the reference (averaged) values.



Typically the derivatives of the deviations are of the same order as the deviations themselves. Consider the deviations having this functional form in 1D: u=Δusinkxxu=kΔucoskx The deviation and its derivative are of order O(Δu)O(1). The convective terms then are: uxu=12Δu2ksin2kx which is order O(Δu2)O(Δu) and therefore negligible compared to other terms.


An analysis of the equations without evaluating the derivatives can be done if we start more general with the continuity and Cauchy momentum equation (neglecting viscous stresses): \partial_{t}\rho+\boldsymbol{\nabla}\cdot\rho\boldsymbol{v}=0


\partial_{t}\rho\boldsymbol{v}+\boldsymbol{\nabla}\cdot\left[\rho\boldsymbol{v}\otimes\boldsymbol{v}+p\boldsymbol{I}\right]=0


linearizing:


\partial_{t}\left(\rho_{0}+\rho'\right)+\boldsymbol{\nabla}\cdot\left(\rho_{0}+\rho'\right)\boldsymbol{v}'=0


\partial_{t}\left(\rho_{0}+\rho'\right)\boldsymbol{v}'+\boldsymbol{\nabla}\cdot\left[\left(\rho_{0}+\rho'\right)\boldsymbol{v}'\otimes\boldsymbol{v}'+\left(p_{0}+p'\right)\boldsymbol{I}\right]=0


i am sure you will agree that \rho' \boldsymbol{v}' \ll \rho_0 \boldsymbol{v}' and \rho'\boldsymbol{v}'\otimes\boldsymbol{v}'\ll\rho_0\boldsymbol{v}'\otimes\boldsymbol{v}'\ll p'\boldsymbol{I}, which yields the linearized equations you are looking for:


\partial_{t}\rho'+\rho_{0}\boldsymbol{\nabla}\cdot\boldsymbol{v}'=0


\rho_{0}\partial_{t}\boldsymbol{v}'=-\boldsymbol{\nabla} p'


Appendix: Using the identity: \boldsymbol{\nabla}\cdot\left(\boldsymbol{A}\otimes\boldsymbol{B}\right)=\boldsymbol{B}\left(\boldsymbol{\nabla}\cdot\boldsymbol{A}\right)+\left(\boldsymbol{A}\cdot\boldsymbol{\nabla}\right)\boldsymbol{B}



we can rewrite the momentum equation in simplified form: \begin{align}\partial_{t}\rho\boldsymbol{v}+\boldsymbol{\nabla}\cdot\left(\rho\boldsymbol{v}\otimes\boldsymbol{v}\right) &=\boldsymbol{v}\partial_{t}\rho+\rho\partial_{t}\boldsymbol{v}+\boldsymbol{v}\boldsymbol{\nabla}\cdot\left(\rho\boldsymbol{v}\right)+\left(\rho\boldsymbol{v}\cdot\boldsymbol{\nabla}\right)\boldsymbol{v} \\ &=\boldsymbol{v}\left[\partial_{t}\rho+\boldsymbol{\nabla}\cdot\left(\rho\boldsymbol{v}\right)\right]+\rho\left[\partial_{t}\boldsymbol{v}+\left(\boldsymbol{v}\cdot\boldsymbol{\nabla}\right)\boldsymbol{v}\right]\end{align}


where the first term on the last line is identically zero due to the continuity equation.


No comments:

Post a Comment

Understanding Stagnation point in pitot fluid

What is stagnation point in fluid mechanics. At the open end of the pitot tube the velocity of the fluid becomes zero.But that should result...