Solutions to the Wave Equation


General Overview

The Wave Equation applies to any uniform medium that provides a restoring force on every particle in the medium proportional to the displacement of that particle from its equilibrium position. In one spatial dimension the medium can be thought of as an ideal string, of uniform density and with no stiffness. In two dimensions the medium will be a membrane with similar properties -- an idealized drum head.

The equation describes a highly symmetrical relationship between the way such displacements vary with time and the way they vary with position inside the medium. The uniform speed of sound in the medium defines the relative scale of that symmetry.

Solutions of the wave equation must fulfill two types of constraints in order to describe any particular situation. Boundary conditions reflect the symmetrical constraints imposed by the shape of the medium and the nature of its boundaries (i.e. the nature of reflections that occur there). Initial conditions define the displacement and velocity of every particle of the medium at some specified instant in time.

The solutions take the form of expressions for the displacement of a particle from its equilibrium position; thus they must be functions of both position and time. Solutions to the wave equation, however, turn out to be products of two terms, one dependent only on time and the other dependent only on position within the medium. This allows wave equation problems to reduce to Eigenvalue problems characterized by families of discrete solutions [a harmonic series of solutions, for instance].

One Spatial Dimension: Normal Modes of an Ideal String

In one dimension, the wave equation may be written

where and are operators acting on the variable , which represents displacement from equilibrium at any point x along the string and any time t. These operators, respectively, perform second partial derivatives with respect to x and t . [In a partial derivative with respect to one variable, all other variables are held constant.] The speed of sound is the constant c relating the temporal and spatial variations in displacement. It is given by , where is the density of the one-dimensional medium (mass/length) and T is its tension.

We expect solutions of the form

where f(x) is independent of t and is an angular frequency (typically measured in radians/sec). Inserting this form into the wave equation gives us

since the factor involving the other variable is just a constant multiplier for the result of each partial derivative operator. Performing the operation on the time-dependent factor, we obtain

with two factors of and one sign change coming from the two differentiations with respect to time. Dividing both sides of the equation by the common cosine factor and by c2 then yields

,

and our task has reduced to solving a simpler equation in a single variable x, with no partial derivatives to worry about. If we define a wave number k such that , we can rewrite the remaining equation as

.

Notice that we have here a result very similar to the one we just described for two successive differentiations -- two factors of a constant and a change in sign -- only this time the derivatives are with respect to x rather than t and the constant is k rather than c. f(x) must involve either a sine or cosine of something proportional to x. If we assume our string has a finite length, i.e. , and impose boundary conditions appropriate to a string fixed at both ends, i.e. , it becomes clear that the trigonometric function in f(x) must be a sine rather than a cosine. [You may want to consider also the effects of different boundary conditions: both ends free to be displaced, or one end fixed and the other free. Both those cases are used in describing the longitudinal displacement of the air in organ pipes.]

This has become an Eigenvalue problem, with a whole familly of discrete solutions for f defining a set of normal vibrational modes for an ideal string:

where n can be any positive integer. Notice that each n defines its own wave number , and angular frequency . The frequencies of these solutions are harmonically related: . Notice how the fundamental frequency depends on the length, density, and tension of the string. The full set of coefficients An define the initial conditions: they can be chosen for any desired shape of the string at t = 0 and describe how much energy is stored in each of the normal modes of the string.

The full solution for our wave equation in one dimension, then -- combining the spatial and temporal factors and summing over all the normal modes -- becomes:

.

Two Spatial Dimensions: Normal Modes of an Ideal Rectangular Membrane

Although there have been square drums, our primary purpose in this section is to emphasize how directly solutions to the one-dimensional wave equation can be generalized to two spatial dimensions. A rectangular symmetry allows all references to position to be in terms of two orthogonal axes that can be dealt with identically. Note the close parallels between what follows and the discussion of the previous section. In rectangular symmetry, the wave equation can be written

.

is again the speed of sound in the medium. The tension T providing the restoring force proportional to the displacement is now surface tension, however -- force per unit length, measured across any line in the two-dimensional membrane's plane -- and the density now has units of mass per unit area. The "ideal" medium again is assumed to have no stiffness. The single spatial operator of the one-dimensional case has become a combination of two orthogonal operators

.

Again we expect solutions of the form

,

In that event, in close analogy to the one-dimensional case, the wave equation may be written

.

Performing the partial derivatives with respect to time we obtain

,

and dividing both sides by the common cosine factor and by c2 we have

,

where we again have defined a wave number k as . Let the dimensions of the membrane be given by and . Then the boundary conditions for fixed edges can be written

.

In a simple extension of the one dimensional case, this again yields an Eigenvalue prblem with solutions

.

The coefficients containing the initial conditions now carry two indices, one for each spatial dimension, as do the wave numbers

.

Each of the indices -- l and m -- can take on any positive integer value. The angular frequencies of each normal mode are given by

,

so the full solution may be written

.

Two Spatial Dimensions: Normal Modes of an Ideal Circular Membrane

In this case it is easiest to use a polar coordinate system with its origin in the center of the membrane. Any point then can be specified in terms of the orthogonal variables , distance from the center where R is the maximum radius of the membrane, and , an angle measured from an arbitarary reference radial line. In terms of these cylindrical coordinates, the wave equation becomes

where again is the speed of sound in the membrane, with surface tension T having dimensions of force/distance and density mass/surface. The spatial derivatives operator becomes, in polar coordinates,

and we again expect solutions of the form

where this time for the sake of convenience we have adopted a complex exponential notation for the time dependence, rather than a cosine function. In close analogy to our earlier cases, such a solution form causes the wave equation to take the form

or, taking the partial derivatives with respect to time,

which can be written

where wave number . The solutions to this Eigenvalue problem are called Bessel functions:

,

with separate factors for the radial and angular parts of the solution, and with the boundary conditions requiring . This limits the solutions to discrete values of kR and, hence, k; to wit.:

where j can be any positive integer. Normal mode frequencies correspond to values of knjR for which Jn(knjR) = 0. Since

,

the frequencies are simply

where c' is defined as the speed of sound in circumferences per unit time. To tabulate normal mode frequency ratios we need only compute

.

The full solution can be written

where the first of the two terms in the factor in square brackets describes the initial position of each particle in the membrane and the second term its initial velocity. One circular node, at radius R, is common to all the normal modes. For j > 1 there will be additional circular nodes, at certain distances from the center of the membrane. For j = 2, for example, the additional circular node is at , while for j = 3 there are two additional circular nodes, at and at .

Departures from Ideal Media

Two primary ways in which real media may differ from the ideal models we have been using are nonzero stiffness in a string or membrane and dissipation. of energy.

Stiffness in a string or membrane provides an additional restoring force, one not proportional to displacement from equilibrium position. This will increase the frequencies of normal modes, an effect that will be relatively stronger for higher order modes because of their shorter wavelengths and, thus, shorter bending radii within the medium. The effect of stiff strings, for instance, is most noticeable in the relative sharpness of higher harmonics.

Energy dissipation causes the amplitude of oscillation to decrease with time. Mechanisms that can cause such dissipation include radiation of sound into adjacent media, and conversion of vibrational energy into heat (due to friction within a real string or membrane, for instance, or viscous heating within a fluid medium).

Three Spatial Dimensions

The wave equation in the three dimensions x, y, and z is a straightforward generalization of the rectangular membrane case discussed above. As with two dimensions, different choices of coordinates can simplify the equations for problems with different symmetries. For sound radiating equally in all directions from a point source, for instance, the use of spherical coordinates allows simplification to a wave equation with the single spatial variable r.