|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| In [[physics]], the '''acoustic wave equation''' governs the propagation of [[acoustic wave]]s through a material medium. The form of the equation is a second order [[partial differential equation]]. The equation describes the evolution of [[acoustic pressure]] <math>p</math> or [[particle velocity]] '''''u''''' as a function of position '''''r''''' and time <math>t</math>. A simplified form of the equation describes acoustic waves in only one spatial dimension, while a more general form describes waves in three dimensions.
| | Hi there. My title is Sophia Meagher although it is not the name on my birth certification. I've always cherished residing in Alaska. To play lacross is some thing he would by no means give up. My day job is an info officer but I've already utilized for another 1.<br><br>My blog phone psychic, [http://kjhkkb.net/xe/notice/374835 kjhkkb.net], |
| | |
| For lossy media, more intricate models need to be applied in order to take into account frequency-dependent attenuation and phase speed. Such models include acoustic wave equations that incorporate fractional derivative terms, see also the [[acoustic attenuation]] article or the survey paper.<ref name="Nasholm2">S. P. Näsholm and S. Holm, "On a Fractional Zener Elastic Wave Equation," Fract. Calc. Appl. Anal. Vol. 16, No 1 (2013), pp. 26-50, DOI: 10.2478/s13540-013--0003-1 [http://arxiv.org/abs/1212.4024 Link to e-print]</ref>
| |
| | |
| ==In one dimension==
| |
| | |
| === Equation ===
| |
| | |
| Feynman<ref name = "Feynman 1">Richard Feynman, Lectures in Physics, Volume 1, 1969, Addison Publishing Company, Addison</ref> derives the wave equation that describes the behaviour of sound in matter in one dimension (position <math>x</math>) as:
| |
| | |
| :<math> { \partial^2 p \over \partial x ^2 } - {1 \over c^2} { \partial^2 p \over \partial t ^2 } = 0 </math>
| |
| | |
| where <math>p</math> is the [[acoustic pressure]] (the local deviation from the ambient pressure), and where <math>c</math> is the [[speed of sound]].
| |
| | |
| ===Solution===
| |
| | |
| Provided that the speed <math>c</math> is a constant, not dependent on frequency (the dispersionless case), then the most general solution is
| |
| | |
| :<math>p = f(c t - x) + g(c t + x)</math>
| |
| | |
| where <math>f</math> and <math>g</math> are any two twice-differentiable functions. This may be pictured as the [[Superposition principle|superposition]] of two waveforms of arbitrary profile, one (<math>f</math>) travelling up the x-axis and the other (<math>g</math>) down the x-axis at the speed <math>c</math>. The particular case of a sinusoidal wave travelling in one direction is obtained by choosing either <math>f</math> or <math>g</math> to be a sinusoid, and the other to be zero, giving
| |
| | |
| :<math>p=p_0 \sin(\omega t \mp kx)</math>.
| |
| | |
| where <math>\omega</math> is the [[angular frequency]] of the wave and <math>k</math> is its [[wave number]].
| |
| | |
| ===Derivation===
| |
| [[File:Derivation of acoustic wave equation.png|400px|thumbnail|Derivation of the acoustic wave equation]]
| |
| The wave equation can be developed from the linearized one-dimensional continuity equation, the linearized one-dimensional force equation and the equation of state.
| |
| | |
| The equation of state ([[ideal gas law]])
| |
| | |
| :<math>PV=nRT</math>
| |
| | |
| In an [[adiabatic process]], pressure ''P'' as a function of density <math>\rho</math> can be linearized to
| |
| | |
| :<math>P = C \rho \,</math>
| |
| | |
| where ''C'' is some constant. Breaking the pressure and density into their mean and total components and noting that <math>C=\frac{\partial P}{\partial \rho}</math>:
| |
| | |
| :<math>P - P_0 = \left(\frac{\partial P}{\partial \rho}\right) (\rho - \rho_0)</math>.
| |
| | |
| The adiabatic [[bulk modulus]] for a fluid is defined as
| |
| | |
| :<math>B= \rho_0 \left(\frac{\partial P}{\partial \rho}\right)_{adiabatic}</math>
| |
| | |
| which gives the result
| |
| | |
| :<math>P-P_0=B \frac{\rho - \rho_0}{\rho_0}</math>.
| |
| | |
| Condensation, ''s'', is defined as the change in density for a given ambient fluid density.
| |
| | |
| :<math>s = \frac{\rho - \rho_0}{\rho_0}</math>
| |
| | |
| The linearized equation of state becomes
| |
| | |
| :<math>p = B s\,</math> where ''p'' is the acoustic pressure (<math>P-P_0</math>).
| |
| | |
| The [[continuity equation]] (conservation of mass) in one dimension is
| |
| | |
| ::<math>\frac{\partial \rho}{\partial t} + \frac{\partial }{\partial x} (\rho u) = 0</math>.
| |
| | |
| Where ''u'' is the [[flow velocity]] of the fluid.
| |
| Again the equation must be linearized and the variables split into mean and variable components.
| |
| | |
| :<math>\frac{\partial}{\partial t} ( \rho_0 + \rho_0 s) + \frac{\partial }{\partial x} (\rho_0 u + \rho_0 s u) = 0</math>
| |
| | |
| Rearranging and noting that ambient density does not change with time or position and that the condensation multiplied by the velocity is a very small number:
| |
| | |
| :<math>\frac{\partial s}{\partial t} + \frac{\partial }{\partial x} u = 0</math>
| |
| | |
| Euler's Force equation (conservation of momentum) is the last needed component. In one dimension the equation is:
| |
| | |
| :<math>\rho \frac{D u}{D t} + \frac{\partial P}{\partial x} = 0</math>,
| |
| | |
| where <math>D/Dt</math> represents the [[convective derivative|convective, substantial or material derivative]], which is the derivative at a point moving with medium rather than at a fixed point.
| |
| | |
| Linearizing the variables:
| |
| | |
| :<math>(\rho_0 +\rho_0 s)\left( \frac{\partial }{\partial t} + u \frac{\partial }{\partial x} \right) u + \frac{\partial }{\partial x} (P_0 + p) = 0</math>.
| |
| | |
| Rearranging and neglecting small terms, the resultant equation becomes the linearized one-dimensional Euler Equation:
| |
| | |
| :<math>\rho_0\frac{\partial u}{\partial t} + \frac{\partial p}{\partial x} = 0</math>.
| |
| | |
| Taking the time derivative of the continuity equation and the spatial derivative of the force equation results in:
| |
| | |
| :<math>\frac{\partial^2 s}{\partial t^2} + \frac{\partial^2 u}{\partial x \partial t} = 0</math>
| |
| | |
| :<math>\rho_0 \frac{\partial^2 u}{\partial x \partial t} + \frac{\partial^2 p}{\partial x^2} = 0</math>.
| |
| | |
| Multiplying the first by <math>\rho_0</math>, subtracting the two, and substituting the linearized equation of state,
| |
| | |
| :<math>- \frac{\rho_0 }{B} \frac{\partial^2 p}{\partial t^2} + \frac{\partial^2 p}{\partial x^2} = 0</math>.
| |
| | |
| The final result is
| |
| | |
| :<math> { \partial^2 p \over \partial x ^2 } - {1 \over c^2} { \partial^2 p \over \partial t ^2 } = 0 </math>
| |
| | |
| where <math>c = \sqrt{ \frac{B}{\rho_0 }}</math> is the speed of propagation.
| |
| | |
| ==In three dimensions==
| |
| | |
| === Equation ===
| |
| | |
| Feynman<ref name = "Feynman 1"/> derives the wave equation that describes the behaviour of sound in matter in three dimensions as:
| |
| | |
| :<math> \nabla ^2 p - {1 \over c^2} { \partial^2 p \over \partial t ^2 } = 0 </math>
| |
| | |
| where <math>\nabla ^2</math> is the [[Laplace operator]], <math>p</math> is the [[acoustic pressure]] (the local deviation from the ambient pressure), and where <math>c</math> is the [[speed of sound]].
| |
| | |
| ===Solution===
| |
| | |
| The following solutions are obtained by [[separation of variables#Partial differential equations|separation of variables]] in different coordinate systems. They are [[phasor (sine waves)|phasor]] solutions, that is they have an implicit time-dependence factor of <math>e^{i\omega t}</math> where <math>\omega = 2 \pi f</math> is the [[angular frequency]]. The explicit time dependence is given by
| |
| :<math>p(r,t,k) = \operatorname{Real}\left[p(r,k) e^{i\omega t}\right]</math>
| |
| Here <math> k = \omega/c \ </math> is the [[wave number]].
| |
| | |
| ====Cartesian coordinates====
| |
| | |
| : <math>p(r,k)=Ae^{\pm ikr} </math>.
| |
| | |
| ====Cylindrical coordinates====
| |
| | |
| : <math>p(r,k)=AH_0^{(1)}(kr) + \ BH_0^{(2)}(kr)</math>.
| |
| | |
| where the asymptotic approximations to the [[Hankel function]]s, when <math>kr\rightarrow \infty </math>, are
| |
| | |
| :<math> H_0^{(1)}(kr) \simeq \sqrt{\frac{2}{\pi kr}}e^{i(kr-\pi/4)}</math>
| |
| | |
| :<math> H_0^{(2)}(kr) \simeq \sqrt{\frac{2}{\pi kr}}e^{-i(kr-\pi/4)}</math>.
| |
| | |
| ====Spherical coordinates====
| |
| | |
| : <math>p(r,k)=\frac{A}{r}e^{\pm ikr}</math>.
| |
| | |
| Depending on the chosen Fourier convention, one of these represents an outward travelling wave and the other an unphysical inward travelling wave. The inward travelling solution wave is only unphysical because of the singularity that occurs at r=0; inward travelling waves do exist.
| |
| | |
| ==References==
| |
| | |
| {{reflist}}
| |
| | |
| == See also ==
| |
| | |
| * [[Acoustics]]
| |
| * [[Acoustic attenuation]]
| |
| * [[Wave Equation]]
| |
| * [[Differential Equations]]
| |
| * [[Thermodynamics]]
| |
| * [[Fluid Dynamics]]
| |
| * [[Pressure]]
| |
| * [[Ideal Gas Law]]
| |
| | |
| {{DEFAULTSORT:Acoustic Wave Equation}}
| |
| [[Category:Acoustics]]
| |
Hi there. My title is Sophia Meagher although it is not the name on my birth certification. I've always cherished residing in Alaska. To play lacross is some thing he would by no means give up. My day job is an info officer but I've already utilized for another 1.
My blog phone psychic, kjhkkb.net,