Thursday, August 31, 2017

waves - How is the relationship between the end correction and the pipe diameter related to acoustic impedance?


I know that Levine and Schwinger calculated the exact value of the end correction by doing something with the acoustic impedance but I don't understand their calculations.


I've looked at paisanco's answer to a similar question, but several points confused me.




  • How can a wave have mechanical impedance? I thought impedance was a property of the medium.





  • Why does the imaginary part of the "radiation impedance of the pipe end" lead to the end correction?




I'm looking for a conceptual explanation rather than a mathematics based one, but any help would be appreciated.



Answer



The reason the end correction impedance is imaginary (not real) is because it represents energy dissipation from the pipe, by radiation into the infinite region of fluid (air) outside the pipe.


Since it dissipates energy, you can think of it physically as a type of damping. All types of damping is represented by imaginary terms in the transfer function - see any textbook or web site on damped harmonic motion for examples.


The impedance is a property of the medium, not of "the wave". The relevant properties include the geometry of the region - i.e. the diameter of the pipe, but not the properties of the wave itself. For example the end correction is independent of the wavelength of sound waves in the pipe.


In a first course (high school level) on sound waves in pipes, you were probably told that at the "open end" of a pipe the pressure is zero and the velocity is non-zero. In fact that is only approximately correct. If it was exactly true, there would be no work done on the air outside the pipe, and you would not hear any sound produced by the pipe. (Work = $\int P\, dV$, which is $0$ if $P = 0$).



Levine and Schwinger produced a series approximation which is a better approximation to the true end condition, but it is still not exact in real world situations. For example, it ignores the viscosity of the air, which creates a boundary layer in the flow inside the pipe and therefore means that the axial flow velocity over a cross-section through the pipe is not constant.


If the wavelength of the wave in the pipe short, with the same order of magnitude as L&S's end correction formula, the approximation breaks down and the L&S formula becomes "completely wrong," but that situation only arises for a very short pipe where the length is similar to the diameter, or for a high harmonic of the oscillation in a longer pipe.


Physically, you can think of the end correction as representing a cylindrical "lump" of air outside the end of the pipe (with length = the end correction, diameter = the pipe diameter) which is forced to oscillate by the vibrating air inside the pipe, and then radiates all its sound energy away into free space.


Note, the math in the L&S derivation is quite "advanced", and their overall approach to solving the problem predates the use of numerical methods (in 1948, computers had recently been invented, but had not yet become commonplace or powerful enough to attack problems like this) - so don't worry about the fact that you don't understand it, unless you are currently studying "theoretical physics" at graduate level! If you set up a computer model to calculate the 3-D flow pattern in and around the pipe using current methods (e.g. the Boundary Element method for acoustic problems) you will get the "end correction" included automatically because of the "no energy reflection" boundary conditions you specify for the flow field infinitely far from the pipe, without having to do anything "special" to include it in the model.


You might find this paper interesting and perhaps more accessible than the original L&S paper: https://arxiv.org/pdf/0811.3625.pdf


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...