There are all sorts of resonances around us, in the world, in our culture, and in our technology. A tidal resonance causes the 55 foot tides in the Bay of Fundy. Mechanical and acoustical resonances and their control are at the center of practically every musical instrument that ever existed. Even our voices and speech are based on controlling the resonances in our throat and mouth. Technology is also a heavy user of resonance. All clocks, radios, televisions, and gps navigating systems use electronic resonators at their very core. Doctors use magnetic resonance imaging or MRI to sense the resonances in atomic nuclei to map the insides of their patients. In spite of the great diversity of resonators, they all share many common properties. In this blog, we will delve into their various aspects. It is hoped that this will serve both the students and professionals who would like to understand more about resonators. I hope all will enjoy the animations.

For a list of all topics discussed, scroll down to the very bottom of the blog, or click here.

Origins of Newton's laws of motion

Non-mathematical introduction to relativity

Three types of waves: traveling waves, standing waves and rotating waves new

History of mechanical clocks with animations
Understanding a mechanical clock with animations
includes pendulum, balance wheel, and quartz clocks

Water waves, Fourier analysis

Sunday, January 31, 2016

The Lagrangian approach to simple waves - several common waves that lack momentum

All postings by author previous: Spherical harmonics up: Contents of this set of postings next: Why only J’s in a standard Bessel expansion?

This posting includes flash animations showing the physics discussed. Most computers have a flash player already installed, but if yours does not, download the free Adobe flash player here.
Flash animations:

The Lagrangian approach to simple waves - several common waves that lack momentum


In this article we demonstrate use of the Lagrangian for a mass on a spring and then for wave propagation in one dimension in several different media. For most of the wave media, the generalized coordinate s is used, attached to the mass elements in place of the usual Cartesian coordinatex. For waves on an elastic cord and for sound waves, we found that the average momentum transported by the waves is zero, in contrast to the usual belief that these waves do possess net momentum. We verified our results with computer simulations which relied on extremely basic equations of physics: F = ma for motion of the mass elements and F = −k Δℓ for the springs (Newton's second law and Hooke's law). This article has several Flash animations as listed above.

The Lagrangian is a mathematical construct used to solve classical mechanics problems in a very general way. It allows use of a very odd coordinate systems chosen for convenience for a particular problem at hand. We demonstrate its use here for several wave media.
Flash animation of mass on a spring
Fig. 1.1. Mass on a spring. Mouse over to see the action. Mouse off it to suspend it and click on it to restart it. If your computer does not support flash click on the video to the left.

1. Lagrangian for a mass on a spring

The Lagrangian is defined to be the difference between the kinetic energy and potential energy in a system. For example for a mass supported on the end of a spring (a classical resonator - see Fig. 1.1) we write:

    ,   (1.1)

where η is the vertical displacement from equilibrium of the mass and k is the spring constant. Note the dot over the η in the first term, standing for a time derivative of η ,  i.e.  .  One oddity of using the Lagrange equation is that we are to treat η and as separate, distinct and unrelated variables for the purpose of differentiation.
We insert our L into the standard Lagrange equation:

    .   (1.2)

Another oddity of this equation is that it specifies that we take derivatives with respect to the unknowns. Most differential equations used in physics are written in terms of derivatives with respect to time and position, not in terms of derivatives with respect to the unknowns such as η and  .  Plugging ahead and inserting L we get:

    ,   (1.3)

which is the standard differential equation governing the motion of a mass on a spring. It has standard oscillatory sinusoidal solutions. For example see this reference.

2. Lagrangian for one dimensional waves in a medium of masses and linear springs

Flash animation of longitudinal waves - mouse over to activate. Mouse off it to suspend it and click on it to restart it.
Fig. 2.1. Animation of longitudinal waves in an array of masses interconnected by linear springs. This is a good model for a continuous stretchy string or spring as long as there are many masses in a wavelength of the wave being observed. Color is used to enhance the visibility of the regions of spring compression and stretching.

Video of longitudinal waves - click on to start

Lumped parameter model for longitudinal waves
Fig. 2.2. Mass and spring model for waves on a stretched slinky. The springs are considered massless and linear, i.e. spring force = −k Δx where k ,  the spring constant which does not change with stretching of the spring and Δx is the amount of stretching of the spring. The numbers on the masses (the rectangles) indicate the value of the generalized coordinate s at those points. The generalized coordinate s is attached to the masses and not to the space as is usual for the Cartesian coordinate x. This model could continue ad infinitum or be finite in length.

The first type of wave we discuss is a longitudinal wave on a stretchy cord or slinky (see Figs. 2.1 and 2.2 above as well as a nice video of longitudinal and transverse waves). As we start out, the math below assumes that there is zero tension on the string of masses and springs when they are in equilibrium. Observing this with a normal slinky would require permentantly prestretching it so that it is not collapsed on itself at zero tension.

We assume that this system can be modeled as masses connected by linear springs as shown in Fig. 2.2. We need to solve for the displacement from equilibrium of each mass.

We use as our ‘‘generalized'' coordinate system a number s equal to the mass number as indicated in Fig. 2.2. As shown in the figure, we see that s = 0 is the location of the left most mass, s = 1 is the location of the 2nd mass, s = 2 is the third mass, etc. Note that this coordinate system is dynamic, the point s = 2 will move around as a passing wave moves the masses around. Furthermore, note that the distortion, ξ, at s = 2 will be the movement from equilibrium of the second mass from the left. In general, ξ will be a function of generalized coordinate s and time t, so we can write ξ(s,t). One of the big advantages of using this generalized coordinate system is that the mass per unit s remains unchanged by passing waves. If we are using a stretchy string, this coordinate could be marked with dots on the string itself. It is also true that the spring constant for a unit s is unaffected by wave motion.

In the Fig. 2.2 we have set s equal to the mass number; however if desired, s could be merely proportional to the mass number. In the following work we assume that the proportionality constant is chosen such that when the string is unstretched s and x are identical. That is for any point along the string, the generalized coordinate s equals the Cartesian coordinate x when the string is in a relaxed state.

Note that the variable for the previously solved mass-on-a-spring η , was only a function of time t, i.e. we might write η(t). Because the linear wave system we are now considering has two independent variables, s and t, we need to broaden Equation (1.2) to include derivatives with respect to both variables [reference]:

    ,      (2.1)

where is defined to be (i.e with respect to time) and ξ' is similarly defined to be the other derivative,  ,  i.e. it is the derivative with respect to the generalized coordinate s.
In this wave medium of masses and linear springs the Lagrangian density is:

    ,      (2.2)

where m1 is the mass per unit length (in coordinate s) of the string assembly, and k1 is the spring constant of a unit length (in s) of the assembly.

Also Δℓ(s,t)/Δs is the change in length of the spring at position s at time t from its non-stretched length x0 , per non-stretched length (i.e. per unit s). We can relate this to the distortion in an incremental system as:

    ,      (2.3)

or in a continuum system we can write this as:

    ,      (2.4)

where, again, η' is defined to be  . This makes our Lagrangian density become:

    .      (2.5)

In (2.2) we measure the displacement η of the masses in regular Cartesean cooordinate increments, instead of the generalized coordinate s. We do this because we know the energy formulas needed in (2.2) in Cartesian form. What is important is that we accurately formulate how the two energies (T1 and V1), and vary with time and s, which we have done.

The Lagrange equation now becomes:

    .      (2.6)

This is the simple wave equation having d'Alembert's solutions:

    .      (2.7)

Equation (2.7) is the standard (i.e. well accepted) solution to the wave equation with one exception: it is written in terms of s, the coordinate attached to the masses, instead of the more usual coordinate x .  The best known solutions are sinusoidal, such as:

    .        (2.8)

Note that this is an exact solution, not just a small amplitude, linearized solution.

Using (2.8) we get the kinetic energy in the wave:

    ,      (2.8a)

which we see has an average value of  .

Similarly,using (2.7), the potential energy is given by:

    ,      (2.9)

which has an average value of  . A little algebra will show that the kinetic energy equals the potential energy (use v = vk1/m1 = ω/κ  for the wave velocity with respect to the coordinate s). The momentum density in a sinusoidal wave is given by:

    .      (2.10)

Because the coordinate s is attached to the masses themselves, the density per unit s does not change due to a passing wave. The momentum is thus sinusoidal and averages to zero. That is  .

Effect of changing from s generalized coordinate to x cartesian coordinate
Fig. 2.3. Graph showing velocity of the masses as a function of the two horizontal coordinate systems: of the generalized coordinate s (in blue) and of the Cartesian coordinate x (shown in green). While the blue trace is a normal sinusoid, the green trace is a distorted sinusoid. The green trace is narrower in the positive half cycle (where the density is greater) and broader in the negative half cycle where the density is less. With respect to average momentum the density variations cancel the distortion effect and produce zero average momentum in these waves.
The density versus x equals m1/[1 - /dx] which is plotted in red. While the density versus x varies over the wave cycle, the density versus s is constant and equals m1

We might take a moment to ponder the difference between a solution in terms of the generalized coordinate s and one in terms of the regular spacial coordinate x. If we have a nice sinusoidal solution in terms of s, we would expect a graph of it to appear as shown in blue in Fig. 2.3. However because s is in general displaced from x by a passing wave, we might expect the wave, when plotted versus x, to appear as plotted in green in the same graph. The equation that relates x to s and ξ is:

x = s + ξ     ,     (2.11)

where for a particular mass, s equals the unstretched x position of that mass and xi; is the amount of stretching.

In terms of momentum, we have an equal mass in the half cycles labeled a and b. In terms of the graph having x as its axis ("b"), we see that the upward half cycle is compacted. Thus, we would expect its mass density to be greater than in the lower half cycle; however at the same time its width is less. The mass in each half cycle are still equal in spite of the mass having been pushed around. Thus the two effects cancel leaving equal and opposite momentums in the two half cycles.

A reader might also note that a given physical point, at a specific time, has a displacement ξ and velocity v = /dt independent of which coordinate system (s or x) we work in.

It would be interesting to show algebraically that the momentum averages to zero in Cartesian coordinates. Such showing involves converting (2.11) from x = f(s) = s + A cos(κs - ωt)  into s = finv(x) . Such inversion is very difficult algebraically. The most straightforward way to accomplish this is numerically. We do this in the Octave simulation section below.

The cancelation of the positive and negative momentum cycles is more obvious in s space where equal Δs always mean equal masses, since the coordinate s is anchored to the mass itself. And since we have perfect sinusoids as solutions in the s space, we can expect exactly zero average momentum.

With static tension: We might wonder how things change if there is a static tension, like when a slinky is stretched between two posts. In this case we would expect a solution (for longitudinal waves) for ξ(s,t) to be something like:

ξ(s,t) = A cos(κs - ωt) + T0s/k1    ,     (2.12)

where T0 is the time average (or static) tension.

The differential equation (2.6) as derived above for this medium involves second derivatives with respect to s and time, t. The new stretching term  (T0s/k1)  in the above solution does not survive the second derivative and so has no effect on the waves in terms of the generalized coordinate s and we can take the above solution (2.12) as an exact solution in generalized coordinates.

At the same time there is an important effect once we convert back to the Cartesian coordinate x. In this converting back we shall used a small amplitude approximation.

The relationship between x and s (shown earlier in (2.11) )  is:

x = s + ξ    ,    (2.13)

where for a particular mass, s equals the unstretched x position of that mass and ξ is the amount of stretching.

Assuming a reasonably small sinusoidal term for the purposes of (2.13), converting between the two coordinate systems s and x, we can approximate ξ as being dominated by the last (static stretching) term in (2.12):

    ,        (2.14)

which means the relation between s and x becomes:

    .        (2.15)

The two media parameters are the spring constant for a unit length in s and the mass m1 in a unit length of s. These equal the spring constant for an unstretched unit length of string and the mass in an unstretched length.

What about these constants in Cartesian coordinates? We will call the spring constant for a unit length in the Cartesian coordinates kx and the mass in a unit length in x as mx . When the string is stretched kx and mx change while k1 and m1 do not change. The relationships between the above parameters are:

    .         (2.16)

The wave velocity in the x coordinate system equals:

    .        (2.17)

In the limit where T0 » k1 these become:


where in the last step of vx we have used  . Note that this result, that  ,  is the standard cited velocity for transverse waves on a string (remembering that we are discussing longitudinal waves here).

When is T0 » k1 ? The stretched length of the string is    so the condition that T0 » k1  will hold if the stretched length is multiple times larger than the unstretched length. If there is preset tension in the string (like many elastic cords have) we must plot the tension versus string length and linearly interpolate back to zero tension to find the mathematical "unstretched length".

Another similar condition (that we will not mathematically deal with here) is where the spring constant k1 is very large (as in the case of using a piano wire or a very stiff fishing line) but the setup has some sort of "elasticity" in the end conditions of the ‘‘string''. This would be the case where a pulley and weight constitute one of the ends or one of the end posts is quite flexible. In either case active potential energy is stored in the weight or flexible end post.

If T0 « k1  then from the equations just above  mx ≈ m1 ,  kx ≈ k1  and  vx ≈ vp . That is to say, the parameters see little change from their unstretched values.

Octave simulation … energy and momentum

The figures below show the setups for two Octave simulations. The springs between the masses (black rectangles) were assumed to be linear. In both Octave simulations we use Newton's second law and Hooks law for a linear spring to calculate the dynamics of wave propagation for the particular setup. Because the simulations rely on extremely basic equations of physics, we believe the simulations are good alternate checks on the results of the Lagrange algebraic methods above.

It is important that we try to achieve single frequency results in the simulations. We do this by exciting the setups in a simple way and watching the resulting waves. The ring resonator is somewhat better because the excitation used there is single frequency and gradual and really lets the medium determines the wave pattern.

Progression of longitudinal waves in a linear array Slow build up of longitudinal rotating waves in a ring resonator

⇑ Fig. 2.4. In the first simulation, a bust of waves were introduced on the left side and observed to propagate towards the right. The simulation ceased when the waves reached the right side. The calculations showed that the potential and kinetic energies were equal. They also demonstrated that the momentum was orders of magnitude smaller than these energies divided by their velocity (the commonly held relationship) and certainly consistent with zero as derived above. Some numbers from this simulation are reported below in Table 2.1. ⇑ Fig. 2.5. In the second simulation, we used a ring resonator setup, where waves traveled round and round the ring, slowly building up. This simulation also showed the potential and kinetic energies to be equal while the momentum was orders of magnitude smaller than the energies divided by the wave velocity. The actual simulation used hundreds of masses around the ring. Some of the results from this simulation are reported below in Table 2.2.
This is a flash animation. To see the action mouse over the image, off to suspend it and click on it to restart it. This animation does not show the slow build up of rotating waves but only the propagation of constant amplitude rotating waves around the ring resonator.
⇑ Fig. 2.6. Progression of small amplitude longitudinal waves in linear simulation setup (shown above). This graph shows the wave at a progression of times, during and after its launch. The length of the simulation medium is 10. An animation of this is shown below in Fig. 2.8. ⇑ Fig. 2.7. Snap shots of longitudinal waves taken at equal time intervals showing the slow build up (towards the end of the simulation) and progression of the rotating waves around the ring resonator. These waves propagate similar to traveling waves and should have similar momentum. The x axis is labeled in terms of distance around the circumference (total circumference equals 10) while the y axis is the wave amplitude. The advantage of this simulation setup is that it avoids the effect of short sinusoidal bursts as opposed to continuous sinusoidal excitation which we used in this ring resonator. An animation of this is shown below in Fig. 2.9.
Fig. 2.8.⇑ Flash animation of a wave burst propagating on a linear medium, such as shown in Fig. 2.4 above. The x displacement is plotted on the y axis and does not mean to imply that the displacement is in the y direction. Mouse over the graph to see the wave propagation. Mouse off it to suspend it and click on it to restart it.
Fig. 2.9.⇑ Flash animation of a rotating wave mode in a ring resonator, shown in Fig. 2.5 above. The f displacement (longitudinal - in the direction around the ring) is plotted on the y axis and does not mean to imply that the displacement is transverse (in the y direction). Mouse over the graph to see the wave propagation. Mouse off it to suspend it and click on it to restart it.
Note that the wave at 2p radians has the same value as that at 0 radians as is required for a circular resonator (since 2p radians is physically the same place as 0 radians). Note also that the waves propagate in the positive direction around the resonator, much the same as traveling waves would. Note also the slow steady build up of the waves due to the exciting source. In the actual Octave simulation, the build up took even more cycles of the resonance.

Table 2.1. Octave calculated values of energy and momentum at four different times during the wave propagation in a linear array. Wave velocity was set equal to 1. These results show that the average momentum is much less than the wave energy divided by the wave velocity, consistent with zero momentum. As in many numerical simulations, the results are unitless, i.e. normalized to some arbitrary value.
Energy in the wave 207, 213, 213, 213
Energy over wave velocity 207, 213, 213, 213
Momentum in the wave -6.1172e-001, -3.2102e-004, -7.6408e-007, 5.3848e-005

Table 2.2. Octave calculated values of energy and momentum at four different times during the wave propagation in a circular array. Wave velocity was set equal to 1. These results show that the average momentum is much less less than the wave energy divided by the wave velocity. The wave amplitude was different than in Table 2.1.
Energy in the wave 5.53, 5.59, 5.61, 5.77
Energy over wave velocity 5.53, 5.59, 5.61, 5.77
Momentum in the wave 1.40e-014, 1.40e-014, 1.40e-014, 1.40e-014

Agreement with other references:
Elmore and Heald and other references imply that all waves transport momentum equal to the wave energy divided by their velocities. This is in disagreement with our findings.

3. Mixed longitudinal and transverse waves on a linear media

Flash animation of mixed wave propagation involving both x and y movement of the masses. Mouse over the figure to activate it. Mouse off it to suspend it and click on it to restart it.
Fig. 3.1. Animation of mixed waves in an array of masses interconnected by linear springs. This is a good model for a continuous stretchy string or spring as long as there are many masses in a wavelength of the wave being observed. If your browser does not support Flash, use the video just below.

Video of mixed wave propagation - click to start and move the cursor off the video.

We add to the above work the possibility that the waves are both longitudinal and transverse. One example of mixed waves are water waves, although they are somewhat more complicated than waves on a stretchy string.

The kinetic energy density for waves having both longitudinal and transverse components is:

    .        (3.1)

The length change of the springs from a horizontal unstretched length is:

    .        (3.1a)

This makes the potential energy density:

    ,        (3.2)

where ξ and η are the x and y displacements of the masses (x and y distances from their equilibrium positions),   and   are their time derivatives, their velocities, and ξ' and η' are their spacial derivates ( and . The length change Δℓs  is the stretching of the springs per unit s of the springs.

Using (3.1) and (3.2) our Lagrangian density can be written:

    .        (3.3)

Because the displacement is now a vector with x and y components, our Lagrange equation is also vectorial (having an equation for each direction). The first Lagrange equation is:

    .        (3.4a)

The second Lagrange equation (for the y direction) is:

    .        (3.4b)

We see that the two components of the Lagrange equation are both rather difficult to unravel and separate out the x and y components (ξ and η) of the motion due to waves.

Two special cases:

1. Pure longitudinal waves will mean (see Equation (3.1a) ):

    ,        (3.5)

and we get the same potential energy density as in Section 2 above ... see (2.4) above for example. The results of calculations with this potential energy density are shown in Section 2 above.

2. Pure transverse waves will mean (again see Equation (3.1a) ):

    .        (3.6)

The Lagrange equation for the vertical motion becomes:

.        (3.7)

For small amplitude waves, the last term can be ignored, leaving the standard waves equation with a wave velocity of

    ,        (3.8)

which is the commonly derived velocity for transverse waves on a string.

3. Strings with a zero unstretched length have particularly simple wave equations. For these (3.2) becomes:

    ,        (3.9)

which makes the Lagrangian:

    .        (3.10)

The two Lagrange equations for x and y directions become:

    .        (3.11)

Note that in contrast to (3.4a) and (3.4b), in Equation (3.11) the x and y components are nicely separated. That is, the first equation has only the x component of the distortion (ξ) as an unknown, while the second equation has only the y component η as an unknown. Thus the longitudinal waves and transverse waves are independent of each other. Furthermore, each of these two equations is a simple wave equation, the same as (2.6) above. As we stated for (2.6) these equations have simple sinusoidal solutions and convey zero average momentum. Just as for (2.6) this equation is written in terms of the generalized coordinate s instead of the usual Cartesian coordinate x.

In order for "zero unstretched length" to have a physical meaning the masses would have to all coexist at the same spot without a rest separation distance. To avoid this we are required to having a static tension, to separate the masses from each other. As discussed in Section 2 above we would use an x solution given by (2.12).

With zero unstretched length we need to alter (2.13) which relates x to s and ξ. In all the above work we assumed that the s coordinate system was set up so that s and x were identical when the string was unstretched, without waves, and at rest. However with the zero unstretched length condition, all the masses have the same x (call it x = 0 )  location when the string is without tension. At the same time the s coordinate must vary through all its range for the various masses. We cannot accomodate this in (2.13). To fix this problem we change (2.13) into:

x = ξ     .        (3.12)

Using the solution (2.12) we see that the coordinate x will take on the form:

x = ξ = A cos(κs - ωt) + T0s/k1 ~ T0s/k1     .        (3.13)

So x will pretty much monotonically increase with s, plus a small oscillating component. The time varying part of the distortion will vary as a normal wave as  ξ(s,t) = A cos(κs - ωt) .  The component of the distortion η (distortion in the y direction, i.e. transverse waves) will have a similar wave function with amplitude, phase and frequency which are independent of those for the x component. These transverse waves will be strongly governed by the static tension (in the x direction) and weaky affected by the longitudinal waves. The longitudinal waves will weakly modulate their x coordinate unless the longitudinal waves have large amplitude creating oscillating tension comparable to the static tension.

As to the momentum content of these waves, the issue is fairly complicated as discussed at great length by Rowland and Pask (scroll down a page to see the text). As the text explains, momentum depends on a lot of factors and detailed constraints. Purely transverse waves have no x  momentum, but in some setups transverse waves will be accompanied by longitudinal waves which may bring momentum. Rowland and Pask maintain that there is a missing factor of 2 in the common assertion that the momentum equals the energy in the wave divided by its velocity. (See Elmore and Heald Section 1.11 as an example of the common formula.) Peskin derives the usual relationship that the momentum for a non-linear vibrating string equals the energy over the velocity (see Equation (100) in his reference) but surprisingly he concludes that the momentum is zero for the linear case which is consistent with my findings just above. He says in the linear case, wave motions in the three cases are easily separable, a fact that I find is true when the equilibrium length is zero.

4. Unstretchable string, no longitudinal mass effect, constant tension for transverse waves:

In this case the potential energy is:

    ,        (4.1)


is the amount of extra string needed because of the transverse wave per unit x.

The kinetic energy in terms of x is complicated by the change in density (per x) because of the transverse waves:

    .        (4.2)

The expression for in terms of x and y  is:

    .        (4.3)

Taking the derivative we will need for the Lagrangian analysis:

    .        (4.4)

The Lagrange equation:

    ,        (4.5)

with  L1 = K1 - V1  becomes:

    .        (4.6)

Taking the limit (small amplitude approximation) where   ,         and also         the above becomes:

    ,        (4.7)

which is a wave equation with the wave velocity given by    as is the usual result.

One interesting fact about a transverse wave traveling on an unstretchable string is that such wave transports mass. To see this consider a burst of such waves (similar to that in Fig. 2.8 above) containing perhaps 10 cycles of the wave with no waves before and after. Inside the burst, there is a longer length of string per Δx (as per Equation(4.3) above) than on either side of the burst. Presumably the source of the waves inserted this extra length (and string mass) during the creation of the burst. Since the burst is traveling at velocity v with extra mass only in the burst region where the waves are, the extra mass must also be traveling at the same velocity. The momentum of this wave is just this extra mass times its velocity. We could easy make the argument that a wave with average momentum along its direction of propagation is always accompanied with a corresponding mass transport. See Rowland and Pask Section VI for more on this.

5. Sound waves

Flash animation of model of masses and massless but spongy gas parcels for sound propagation. Mouse over the figure to activate it. Mouse off it to suspend it and click on it to restart it.
Fig. 5.1. Animation of sound propagation through model of masses (black rectangles) and massless spongy gas parcels (in red and blue). This model is ok as long as there are many masses in a wavelength of the wave being observed. Color has been added to indicate the amount of compression in the gas parcels and make the wave more visible. This animation assumes that there is an invisible source of waves on the left and an absorber of waves on the right.

Sound through a gas has the same kinetic energy term as above, but the potential energy is different. For adiabatic sound transmission we have the following equation:

Comparison of adiabatic ideal gas and linear spring
Fig. 5.2. The blue line shows the relative pressure of a gas parcel undergoing adiabatic compression versus x, the motion of a compressing piston. This is compared with the restoring force for an equivalent linear spring (in green), made to be equivalent around the point x = 1. We see that while the first derivatives (the slopes) of the two curves are the same around this point, their second derivatives (their curvatures) differ greatly.

    ,        (5.1)

where V = Ax is the volume of gas under consideration and γ = Cp/CV is the gas constant. Also, A
is the cross sectional area of the sound duct (assumed constant) and x0 is the length of a gas parcel.

At an ambient pressure P0 this becomes:

    ,        (5.2)

where V0 and x0 are the volume and length of a gas parcel under ambient pressure and temperature. Putting these equations together yields

    ,        (5.3)

that is to say, the pressure is inversely proportional to the volume to the gamma power or alternatively, inversely proportional to the length of the gas parcel to the gamma power.

To calculate the potential energy in a gas parcel, we consider a gas parcel in isolation where we have allowed its length x to become very long. Then we calculate the energy required to compress the gas to make its length equal l:

    .      (5.4)

We could define the potential energy we use here as shown in (5.4) minus a reference energy, such as the potential energy when the parcel's pressure equals P0 .  On the other hand this reference energy term would not survive the differentiation step in the Lagrange equation and so in the end makes no difference to the final equation. This is typical of potential energies: we can usually choose our reference, i.e. add an arbitrary constant, without affecting the outcome of our calculations.

We need the potential energy density per unit s, similar to m1 and k1 .  This potential energy density is given by:

    ,      (5.5)

where is the length of a gas parcel whose length at ambient pressure is x0 .  Also, ξ' is defined as /ds .  So our Lagranian becomes

    .        (5.6)

We substitute this into

    ,        (5.7)

to get

    ,      (5.8)


    .      (5.8aa)

Rearranged this becomes:

    .      (5.8a)

This equation has the feature that if ξ has small amplitude, the ξ' in the denominator in the right hand side can be ignored, i.e. the small signal linearization. This equation says that the small signal wave velocity is  ,   consistent with Elmore and Heald.

Next we investigate second order behavior of the waves predicted by (5.8a) using a power series solution of (5.8a).

The general binomial equation is:

    ,      (5.9)

or when x = 1  and y = e and e is much smaller than 1, this becomes:

    ,      (5.9a)

to the second power in e.

When we apply (5.9a) to the right hand side of Eqn. (5.8a), we get:

    .        (5.10)

This makes (5.8a) become:

    .      (5.10a)

A power series solution of sinusoids up through square terms is:

    ,      (5.11)


    ,      (5.12)

where X = κs - ωt . We only have the cosine term (to the first power) and not the sine because we are free to pick the phase of our solution, at least the phase of the first order sinusoids.

We now substitute (5.12) into (5.10a) one side at a time. On the left hand side of (5.10a), the a0 term disappears while the a1 term stays linear in cosX. All the other terms come through as second order in sinusoids of X and so are important (remember we are doing this to second order). The left side of (5.10a) now looks like:

    .      (5.13)

Next we do the right hand side of (5.10a) remembering that we are only keeping terms of power two or less. The right-most factor automatically gives us one power (or more), so we only need to use the linear term in . The right most factor (not the whole right hand side) is:

    ,      (5.14)

which is essentially the same as the previous formula (5.13).

The important part of the other factor on the right hand side is:

    .      (5.15)

After multiplying two factors together, (5.14) and (5.15), the right hand side of (5.10a) becomes:

    .      (5.16)

Equating the linear terms of the two sides (in (5.13) and (5.16) ) we get:

    ,      (5.17)

which gives the relationship between frequency and wavenumber:

    ,      (5.18)

which is the wave velocity.

Development of shocking
Fig. 5.3. We show the addition of higher order terms to the basic sinusoid (in blue) to create shock waves. The green line shows the result of adding a sinXcosX = 0.5sin2X term to the basic sinusoid as discussed in the text to the left. The red line shows the eventual fully developed shock with a great number of higher order terms. See an early blog for the complete infinite Fourier series of a saw tooth wave. One can also see this general reference on shocking of sound waves.

Equating the cos2X terms, we have:

    ,      (5.19)

which yields:

    ,      (5.20)

and is the same as before and doesn't tell us anything new.

The sin2X terms (in (5.13) and (5.16) ) give us:

    ,      (5.21)

again yielding the same relationship for the third time.

Equating the sinXcosX terms yields:

    ,      (5.22)

which can be solved:

    .      (5.23)

The addition of a sinXcosX term is the beginnings of the formation of shock waves. This behavior of sound is illustrated in Fig. 5.3. Shocking is a distortion of the sinusoid acoustic wave that results from the fluid parcels in the forward moving part of the sinusoid moving faster than the rest of the wave. This phenomenon takes place at high acoustic power levels.

Momentum is given by  :

    .      (5.24)

The cos2X factor averages to zero, implying that the average momentum in the s coordinate system is zero. As we saw with respect to (2.10) above, this means the average momentum is also zero in the Cartesian x coordinate system. This result disagrees with Peskin under "sound waves" Equation (40). Our results do agree with the findings of Rowland and Pask Section VI.   Rowland and Pask also discuss pseudomomentum which has applications in quantum mechanics. We need to remember that the above algebra is only been done to second order in the sinusoids. Higher order calculation, applicable to very intense sound, might yield a non-zero result. On the other hand, we have shown that momentum will not be proportional to the energy over the wave velocity.

We next discuss a computer simulation to confirm the zero result.

Octave simulation for sound … energy and momentum

To confirm some of the algebraic results of this section, we ran Octave simulations similar to those discussed at the end of Section 2 above. Because the circular or ring configuration was shown to be more accurate (see Section 2 above), we restricted the simulation to the circular configuration. Two of the figures and animations above in Section 2 (Fig. 2.5 and Fig. 2.9) are appropriate for this section. Compared with the simulations used for Section 2 we only changed the interaction force between adjacent masses to be appropriate for gas parcels, Eqn (5.3), instead of that appropriate for linear springs.

The results of the calculations of the energies, energy over waves velocity and momentum are shown in Table 5.1. We see that the momentum is orders of magnitude smaller than the energy over velocity, implying that the momentum is zero.

Table 5.1. Octave calculated values of energy and momentum at four different times during the wave propagation on a circular array. The wave velocity was set equal to 1.18 .  The amplitude was different than that in Section 2 above. These results show that the average momentum is much much less than the wave energy divided by the wave velocity.
Energy in the wave 7.6814e-006,  3.0705e-005,  6.8912e-005,  1.2211e-004
Energy over wave velocity 6.4920e-006,  2.5951e-005,  5.8241e-005,  1.0321e-004
Momentum in the wave -2.5324e-015,  -2.5307e-015,  -2.5278e-015.  -2.5322e-015

All postings by author previous:Spherical harmonics up: Contents of this set of postings next: Why only J’s in a standard Bessel expansion?