Nope, I'm not Dinosaur. I'm the Chuck that wrote the website Chuck's Chimes or perhaps better known as
I'd have to dig out my notes from the class but off the top of head NASTRAN (?) can use two different formulations for a fluid. For air, other than at resonance with the structure, a light (less complex) formulation can be used. I would assume this formulation would require less computer time to solve.
I had not thought about temperature. If I could do this at work this would be rather easy to do; however, I don't think my employer would be willing to support my hobby.
I would point out again that author of the acoustic book I contacted didn't appear to have a problem with coupling of the axial air column and transverse vibration. Here is what he wrote.
I'm afraid my book was the wrong place to go! It deals with the production of sound by flow, i.e. where the energy is extracted from the flow.
>From my brief look at the your website (and another linked on your site)
I gather that the excitation mechanism is a striker. In that case I
should think sound generation is certainly via the resonant excitation
of a resonant acoustic mode in the tube by a structural mode of the
tube. The resistance of the air is so small that it seems to me that the
method of calculation you appear to have been using is quite adequate to
determine the frequency.
If you want to express the acoustic amplitude in terms of the tube
motion, however, you will have to calculate the tube mode shape and use
it to determine the normal velocity distribution on the inner surface of
the tube as a function of time (the amplitude will decay, predominantly
because of structural damping). You would then have to solve the
acoustic wave equation within the tube (say with zero pressure open end
conditions, in a first approximation) for the sound generated by a point
source. Then imagine the inner surface divided into infinitesimal
surface elements, the normal velocity of each of which determines its
effective source strength, and calculate the sound produced by each of
these sources using your point source solution. Then integrate over the
inner area to get the net effect of all the sources. The result (i.e.
the solution of the wave equation due to surface forcing) will have a
large peak at the matching resonant frequency of the air in the tube,
the magnitude being dependent on the decay rate of the tube vibration.
Next use the solution to calculate the fluctuating mean axial volume
velocity of the air at the open ends, and use this as the effective
source strength of the two monopole sources (one at each end) of the
sound radiating into the ambient free pace.
chuck