Wave equation (classical physics): Difference between revisions
imported>Paul Wormer |
imported>Paul Wormer |
||
(2 intermediate revisions by the same user not shown) | |||
Line 34: | Line 34: | ||
It is to be shown that the product function | It is to be shown that the product function | ||
:<math> | :<math> | ||
h | h \equiv f\circ g\quad \Longleftrightarrow\quad h(x,t) \equiv f\big(g(x,t)\big) \equiv f\big(x-vt\big) | ||
</math> | </math> | ||
is a solution of the wave equation for ''any'' real twice-differentiable function ''f'' of a single variable. Indeed, by the [[chain rule]] one gets, | is a solution of the wave equation for ''any'' real twice-differentiable function ''f'' of a single variable. Indeed, by the [[chain rule]] one gets, | ||
Line 101: | Line 101: | ||
==Normal modes of a string== | ==Normal modes of a string== | ||
The archetypal example of a wave is the transverse motion of a string constrained to vibrate in a plane. Let ''y''(''x'',''t'') be the transverse displacement of the string from its equilibrium. The wave equation is | The archetypal example of a wave is the transverse motion of a string constrained to vibrate in a plane. Let ''y''(''x'',''t'') be the transverse displacement of the string from its equilibrium. The wave equation is | ||
:<math> | :<math> | ||
Line 133: | Line 130: | ||
</math> | </math> | ||
General solutions are linear combinations of these functions. | General solutions are linear combinations of these functions. | ||
{{Image|SW1.gif|right||<small>Fig. 2. Fundamental mode, frequency ν<sub>0</sub>.</small>}} | |||
{{Image|SW2.gif|right||<small>Fig. 3. First overtone, frequency 2ν<sub>0</sub>.</small>}} | |||
{{Image|SW3.gif|right||<small>Fig. 4. Second overtone, frequency 3ν<sub>0</sub>.</small>}} | |||
{{Image|Moving wave.gif|right||<small>Fig. 5. Linear combination of lowest 3 modes </small> }} | |||
Let us next assume that the string is clamped for ''x'' = 0, that is, the displacement is zero for ''x'' = 0. This boundary condition excludes a cosine type solution, so that the spatial equation is of the form sin''kx''. If the string has length ''L'' and its end is also clamped, it follows that | Let us next assume that the string is clamped for ''x'' = 0, that is, the displacement is zero for ''x'' = 0. This boundary condition excludes a cosine type solution, so that the spatial equation is of the form sin''kx''. If the string has length ''L'' and its end is also clamped, it follows that | ||
:<math> | :<math> | ||
Line 145: | Line 145: | ||
Here λ, a quantity with a dimension of length, is called [[wave length]]. In figures 2, 3. and 4: λ= 2''L''/1, 2''L''/2, and 2''L''/3, for the first vibration (the "fundamental"), the second (the "first overtone"), and the third (the "second overtone"), respectively. These sinusoidal modes of the string are its ''normal modes''. | Here λ, a quantity with a dimension of length, is called [[wave length]]. In figures 2, 3. and 4: λ= 2''L''/1, 2''L''/2, and 2''L''/3, for the first vibration (the "fundamental"), the second (the "first overtone"), and the third (the "second overtone"), respectively. These sinusoidal modes of the string are its ''normal modes''. | ||
With regard to the temporal equation we remark the following. The string is pulled to | With regard to the temporal equation we remark the following. The string is pulled to certain shape ''y''(''x''), held in position for a brief instant, and then let go at ''t'' = 0, so that ''dy/dt'' = 0 at ''t'' = 0 and ''dy/dt'' ≠ 0 for ''t'' > 0, | ||
:<math> | :<math> | ||
\left(\frac{\partial y(x,t)}{\partial t}\right)_{t=0} = 0. | \left(\frac{\partial y(x,t)}{\partial t}\right)_{t=0} = 0. | ||
Line 162: | Line 162: | ||
\nu = \frac{1}{t_0} = \frac{\omega}{2\pi} = \frac{n}{2L} \sqrt{\frac{T}{\rho}} \equiv n \nu_0 , | \nu = \frac{1}{t_0} = \frac{\omega}{2\pi} = \frac{n}{2L} \sqrt{\frac{T}{\rho}} \equiv n \nu_0 , | ||
</math> | </math> | ||
where ν<sub>0</sub> ≡ '' | where ν<sub>0</sub> ≡ ''1''/(2L) √T/ρ is the fundamental frequency. | ||
One may wonder whether the solution has the general form found by d'Alembert and discussed in the previous section. Note to that end that | |||
:<math> | |||
\sin(k x)\cos(\omega t) = \frac{1}{2}\left[ \sin k(x+ v t) + \sin k(x-vt) \right], | |||
</math> | |||
and that the sine function is odd, sin(−''u'') = −sin(''u''). | |||
===Examples=== | |||
In figures 2, 3, and 4 are shown the functions | In figures 2, 3, and 4 are shown the functions | ||
:<math> | :<math> | ||
Line 169: | Line 176: | ||
</math> | </math> | ||
which implies the following choices: (i) the initial displacement (amplitude) ''A'' = 0.3, (ii) the length of the string ''L'' = π, (iii) the propagation velocity ''v'' = 1, and (iv) the radial velocity ω = ''kv'' = ''n''. These choices are compatible: the choice of ''L'' fixes a unit of length and the choice of ''v'' fixes a unit of time; these units can be fixed independently. | which implies the following choices: (i) the initial displacement (amplitude) ''A'' = 0.3, (ii) the length of the string ''L'' = π, (iii) the propagation velocity ''v'' = 1, and (iv) the radial velocity ω = ''kv'' = ''n''. These choices are compatible: the choice of ''L'' fixes a unit of length and the choice of ''v'' fixes a unit of time; these units can be fixed independently. | ||
The <i>d'Alembertian operator</i> defined by | |||
:<math> | |||
\square^2 \quad \stackrel{\mathrm{def}}{=}\quad \frac{\partial^2}{\partial x^2} - \frac{1}{v^2}\frac{\partial^2}{\partial t^2} | |||
</math> | |||
is linear. Hence | |||
:<math> | :<math> | ||
\sin | \square^2 \Big(\sin x\cos t + 0.7 \sin 2x\cos 2 t + 0.5 \sin 3x\cos 3 t \Big) = 0. | ||
</math> | </math> | ||
This solution of the wave equation is shown in figure 5. | |||
==Acoustical waves== | ==Acoustical waves== |
Latest revision as of 02:27, 24 March 2010
![](http://s9.addthis.com/button1-share.gif)
In classical physics, the wave equation is the name given to a certain real partial differential equation in which the second derivative with respect to position x is proportional to the second derivative with respect to time t. Hence, the wave equation has the general form
The real proportionality constant v2 had dimension (length)2 over (time)2, i.e., v has the dimension of speed. In the simplest case v is independent of x and t (homogeneous space) and independent of f (no dispersion).
The wave equation can be generalized to three dimensions,
The quantity v has again dimension of speed (length over time); ∇2 is the Laplace operator; t is the time; r represents the space coordinates, r ≡ (x, y, z). In classical physics it is commonly the case that the function f(r,t) is real and observable. Examples of the meaning of f and v will be discussed below. In an isotropic space v does not depend on direction.
History
In 1727 John (Johann) Bernoulli studied the vibrations of a discrete number of masses connected by springs that satisfy Hooke's law. He derived a difference equation for the displacement ψk(t) of the kth mass as a function of time. In 1746 the work was extended by Jean Le Rond d'Alembert in his attack on the vibrating string problem typified by the violin string.[1] To treat the continuous string d'Alembert increased the number of masses, letting their number go to infinity, while the size and mass of each individual "bead" went to zero, such that the mass density ρ of the string (mass per unit length) remained finite and constant over the string. In this manner d'Alembert derived the equation
where T is the tension (force along the string, also taken constant over the string). He derived a solution of the form
and emphasized that infinitely many curves f and g other than sinusoidal curves are modes of vibration. The 1747 article contains the first appearance of the wave equation in print. Note that the dimension of ρ/T is: M/L (M L/t2)−1 = (t/L)2, where the force T (mass times acceleration) has dimension ML/t2. In total, ρ/T has indeed the dimension of one over speed squared.
Solution of the one-dimensional equation
Any arbitrary differentiable function of the form f(x − vt) is a solution of the one-dimensional wave equation for constant (independent of x and t), real, and positive v. Hence any function f of x − vt may be called a wave. In contrast to what is often assumed, a wave (solution of the wave equation) does not have to be periodic, so that concepts like wave length and period are not necessarily defined. An example is given in Figure 1, where the red curve depicts an arbitrary wave propagating along the x-axis in positive direction. Snapshots at the times t=0, t=t0, and t=2t0 are depicted. At t=0 the red curve is the function f(x) that has non-zero values only in the neighborhood of x=0. At t=t0 the red curve has the very same shape, but has non-zero values only in the neighborhood of x=vt0. Note that the argument of f in the latter point on the x-axis is: x−vt0 = 0.
In order to prove that f(x − vt) satisfies the wave equation, we write
It is to be shown that the product function
is a solution of the wave equation for any real twice-differentiable function f of a single variable. Indeed, by the chain rule one gets,
and
and also
so that
and it follows that h(x,t) ≡ f(x−vt) satisfies the wave equation.
The same steps can be followed for f(x + vt). Hence f(x±vt) satisfies the wave equation. In a notation less strict than just used—but notwithstanding obvious—this can be written as
Note that f(x+vt) for x=−vt has the same value as f(0) at t=0. In other words, the function moves along the negative x-axis.
The initial condition, the form f(x) of f at t = 0, specifies (partly) the wave. [2] Because f(x ± vt) differs from f(x) only by a translation of magnitude vt along the x-axis, (−vt in positive direction, +vt in negative direction), it is apparent that the function moves along the x-axis with velocity v, the velocity of wave propagation.
A general solution of the one-dimensional wave equation (first found by d'Alembert) is
with
where f and g are arbitrary twice differentiable functions of one variable.
In order to specify physical solutions the wave equation must be supplemented by boundary conditions. Let us impose the conditions:
When physically ψ(x,t) stands for the motion of a string of length L, this condition means that both ends of the string are clamped, which is the case for musical string instruments. In order to satisfy the boundary conditions, the functions f and g must be related by the equations
By putting f(u) = g(u) for all u ≡ vt and requiring f to be an odd function,
the first boundary condition is satisfied:
The second boundary condition enforces periodicity of f with wave length 2L
where the substitution vt = −(L+u) was made to obtain the expression after the arrow.
Notice that the part of f(x) lying between x = 0 and x = L is arbitrary, which corresponds to the physical fact that a string can be started vibrating from any shape whatsoever (subject to its ends being fixed). In particular, its shape need not be sinusoidal, proving that solutions other than the normal modes (to be discussed in the next section) can occur.
Normal modes of a string
The archetypal example of a wave is the transverse motion of a string constrained to vibrate in a plane. Let y(x,t) be the transverse displacement of the string from its equilibrium. The wave equation is
where, as first shown by d'Alembert, the velocity v is √T/ρ. The tension T as well as the linear mass density ρ (mass per length) are constant over the string. The wave equation is derived under the assumption that the slope of the spring is small: . This holds if the maximum displacement of the spring is much less than its length.
The equation may be solved by a method known as separation of variables. One makes the Ansatz
where F depends on x only and G on t only. Substitution of the Ansatz into the wave equation gives
The equation on the right has a left-hand side solely depending on x, while its right-hand side depends solely on t. This can only be correct if both equations are equal to a constant, independent of x and t. Let us write −k2 for this constant, where k may be complex in principle. We will, however, assume k to be real and hence will retrieve a restricted set of solutions, the normal modes. The wave equation separates into two equations,
Clearly k has a dimension of reciprocal length and vk has a dimension of reciprocal time. It is common to write
and to call ω the angular velocity. Special solutions of the spatial and temporal equations are:
General solutions are linear combinations of these functions.
Let us next assume that the string is clamped for x = 0, that is, the displacement is zero for x = 0. This boundary condition excludes a cosine type solution, so that the spatial equation is of the form sinkx. If the string has length L and its end is also clamped, it follows that
It is common to write
Here λ, a quantity with a dimension of length, is called wave length. In figures 2, 3. and 4: λ= 2L/1, 2L/2, and 2L/3, for the first vibration (the "fundamental"), the second (the "first overtone"), and the third (the "second overtone"), respectively. These sinusoidal modes of the string are its normal modes.
With regard to the temporal equation we remark the following. The string is pulled to certain shape y(x), held in position for a brief instant, and then let go at t = 0, so that dy/dt = 0 at t = 0 and dy/dt ≠ 0 for t > 0,
This condition excludes the function sinωt, because the first derivative of the sine is the cosine which is not zero at t = 0. Hence the time function is of the form cosωt. In total, the solution with given boundary and initial conditions is
A function of this form is a standing wave. Clearly, this function repeats itself after a time period t0 with
where T is the tension and ρ is the mass density. The reciprocal period is the frequency,
where ν0 ≡ 1/(2L) √T/ρ is the fundamental frequency.
One may wonder whether the solution has the general form found by d'Alembert and discussed in the previous section. Note to that end that
and that the sine function is odd, sin(−u) = −sin(u).
Examples
In figures 2, 3, and 4 are shown the functions
which implies the following choices: (i) the initial displacement (amplitude) A = 0.3, (ii) the length of the string L = π, (iii) the propagation velocity v = 1, and (iv) the radial velocity ω = kv = n. These choices are compatible: the choice of L fixes a unit of length and the choice of v fixes a unit of time; these units can be fixed independently. The d'Alembertian operator defined by
is linear. Hence
This solution of the wave equation is shown in figure 5.
Acoustical waves
To a very good approximation the propagation of pressure waves in a fluid (liquid or gaseous) medium is described by the wave equation. Such waves are called acoustical waves, because the propagation of sound in a medium is an example.
Consider an undisturbed fluid particle at position x and let its displacement, caused by a disturbance at time t, be in the z direction, so that after being disturbed the particle is at position: x + z(x,t). If the slope satisfies (the "acoustic approximation") the displacement z satisfies
where
Here ρ is the density of the fluid, ρ0 the undisturbed density, and P(ρ) the pressure in the fluid.
If the fluid is air, then is the speed of sound. There may be 3-dimensional motion in the general case, satisfying the 3-dimensional wave equation, with an appropriate modification of the parameters. However, if is the fluid is homogeneous and isotropic, the extension is without any difficulties.
Electromagnetic waves
In a region of space where the charge and current densities are zero, Maxwell's equations for the electromagnetic field in vacuum lead to
where E is the electric field and B is magnetic induction (also known magnetic flux density). The quantity c is the speed of light in vacuum. The vector potential A and the scalar potential φ give rise to both E and B through
where in SI units k = 1 and in Gaussian units k = c. Both A and φ satisfy wave equations (that hold in absence of charge and current densities); they are
Notes
- ↑ J. d'Alembert, Recherches sur la courbe que forme une corde tenduë mise en vibration [Researches on the curve formed by a stretched string set into vibration], Histoire de l’Académie des sciences et belles-lettres (Berlin) vol. 3, pp. 214–249 (1747, published 1749).
- ↑ Since the differential equation is of second order, a second initial condition, for instance the time derivative of f at t = 0, is required for the general solution of the wave equation.