Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Mathematics LibreTexts

2.7: d’Alembert’s Solution of the Wave Equation

( \newcommand{\kernel}{\mathrm{null}\,}\)

A general solution of the one-dimensional wave equation can be found. This solution was first Jean-Baptiste le Rond d’Alembert (1717- 1783) and is referred to as d’Alembert’s formula. In this section we will derive d’Alembert’s formula and then use it to arrive at solutions to the wave equation on infinite, semi-infinite, and finite intervals.

We consider the wave equation in the form utt=c2uxx and introduce the transformation

u(x,t)=U(ξ,η),whereξ=x+ctandη=xct.

Note that ξ, and η are the characteristics of the wave equation.

We also need to note how derivatives transform. For example

ux=U(ξ,η)x=U(ξ,η)ξξx+U(ξ,η)ηηx=U(ξ,η)ξ+U(ξ,η)η.

Therefore, as an operator, we have

x=ξ+η.

Similarly, one can show that

t=cξcη.

Using these results, the wave equation becomes

0=uttc2uxx=(2t2c22x2)u=(t+cx)(tcx)u=(cξcη+cξ+cη)(cξcηcξcη)U=4c2ξηU.

Therefore, the wave equation has transformed into the simpler equation,

Uηξ=0.

A further integration gives

U(ξ,η)=ηΓ(η)dη+F(ξ)=G(η)+F(η).

Therefore, we have as the general solution of the wave equation,

u(x,t)=F(x+ct)+G(xct),

where F and G are two arbitrary, twice differentiable functions. As t is increased, we see that F(x+ct) gets horizontally shifted to the left and G(xct) gets horizontally shifted to the right. As a result, we conclude that the solution of the wave equation can be seen as the sum of left and right traveling waves.

Let’s use initial conditions to solve for the unknown functions. We let

u(x,0)=f(x),ut(x,0)=g(x),|x|<.

Applying this to the general solution, we have

f(x)=F(x)+G(x)g(x)=c[F(x)Gx)].

We need to solve for F(x) and G(x) in terms of f(x) and g(x). Integrating Equation (2.7.13), we have

1cx0g(s)dx=F(x)G(x)F(0)+G(0).

Adding this result to Equation (2.7.13), gives

F(x)=12f(x)+12cx0g(s)ds+12[F(0)G(0)].

Subtracting from Equation (2.7.13), gives

G(x)=12f(x)12cx0g(s)ds12[F(0)G(0)].

Now we can write out the solution u(x,t)=F(x+ct)+G(xct), yielding d’Alembert’s solution

u(x,t)=12[f(x+ct)+f(xct)]+12cx+ctxctg(s)ds.

When f(x) and g(x) are defined for all xR, the solution is well-defined. However, there are problems on more restricted domains. In the next examples we will consider the semi-infinite and finite length string problems.In each case we will need to consider the domain of dependence and the domain of influence of specific points. These concepts are shown in Figure 2.7.1. The domain of dependence of point P is red region. The point P depends on the values of u and ut at points inside the domain. The domain of influence of P is the blue region. The points in the region are influenced by the values of u and ut at P.

clipboard_e1f1c237d80de5831ebe73d7955bd73df.png
Figure 2.7.1: The domain of dependence of point P is red region. The point P depends on the values of u and ut at points inside the domain. The domain of influence of P is the blue region. The points in the region are influenced by the values of u and ut at P.

Example 2.7.1

Use d’Alembert’s solution to solve

utt=c2uxx,u(x,0)=f(x),ut(x,0)=g(x),0x<.

Solution

The d’Alembert solution is not well-defined for this problem because f(xct) is not defined for xct<0 for c, t>0. There are similar problems for g(x). This can be seen by looking at the characteristics in the xt-plane. In Figure 2.7.2 there are characteristics emanating from the points marked by η0 and ξ0 that intersect in the domain x>0. The point of intersection of the blue lines have a domain of dependence entirely in the region x, t>0, however the domain of dependence of point P reaches outside this region. Only characteristics ξ=x+ct reach point P, but characteristics η=xct do not. But, we need f(η) and g(x) for x<ct to form a solution.

clipboard_ecceb88db28a993cc1eea9c6225fc5f4c.png
Figure 2.7.2: The characteristics for the semi-infinite string.

This can be remedied if we specified boundary conditions at x=0. For example, Fixed end boundary condition we will assume the end x=0 is fixed,

u(0,t)=0,t0.

Imagine an infinite string with one end (at x=0) tied to a pole.

Since u(x,t)=F(x+ct)+G(xct), we have

u(0,t)=F(ct)+G(ct)=0.

Letting ζ=ct, this gives G(ζ)=F(ζ), ζ0.

Note that

G(ζ)=12f(ζ)12cζ0g(s)dsF(ζ)=12f(ζ)12cζ0g(s)ds=12f(ζ)+12cζ0g(σ)dσ

Comparing the expressions for G(ζ) and F(ζ), we see that

f(ζ)=f(ζ),g(ζ)=g(ζ).

These relations imply that we can extend the functions into the region x<0 if we make them odd functions, or what are called odd extensions. An example is shown in Figure 2.7.3.

Another type of boundary condition is if the end x=0 is free,

ux(0,t)=0,t0.

In this case we could have an infinite string tied to a ring and that ring is allowed to slide freely up and down a pole.

One can prove that this leads to

f(ξ)=f(ξ),g(ξ)=g(ξ).

Thus, we can use an even extension of these function to produce solutions.

Example 2.7.2

Solve the initial-boundary value problem

utt=c2uxx,0x<,t>0.u(x,0)={x,0x1,2x,1x2,0x<0,x>2,ut(x,0)=0,0x<.u(0,t)=0,t>0.

Solution

This is a semi-infinite string with a fixed end. Initially it is plucked to produce a nonzero triangular profile for 0x2. Since the initial velocity is zero, the general solution is found from d’Alembert’s solution,

u(x,t)=12[f0(x+ct)+f0(xct)],

where f0(x) is the odd extension of f(x)=u(x,0). In Figure 2.7.3 we show the initial condition and its odd extension. The odd extension is obtained through reflection of f(x) about the origin.

clipboard_ec3aaf0e322f77f0b216096a66b895844.png
Figure 2.7.3: The initial condition and its odd extension. The odd extension is obtained through reflection of f(x) about the origin.

The next step is to look at the horizontal shifts of f0(x). Several examples are shown in Figure 2.7.4.These show the left and right traveling waves.

clipboard_e666a6f55da4b721be3b87a91a814e465.png
Figure 2.7.4: Examples of f0(x+ct) and f0(xct).

In Figure 2.7.5 we show superimposed plots of f0(x+ct) and f0(xct) for given times. The initial profile in at the bottom. By the time ct=2 the full traveling wave has emerged. The solution to the problem emerges on the right side of the figure by averaging each plot.

clipboard_e496c99d4bb855adedbba35e74ca4e82f.png
Figure 2.7.5: Superimposed plots of f0(x+ct) and f0(xct) for given times. The initial profile in at the bottom. By the time ct=2 the full traveling wave has emerged.
clipboard_ea4a1b0f235c00720b8d36d1b13e9b818.png
Figure 2.7.6: On the left is a plot of f(x+ct), f(xct) from Figure 2.7.5 and the average, u(x,t). On the right the solution alone is shown for ct=0 at bottom to ct=1 at top for the semi-infinite string problem

Example 2.7.3

Use d’Alembert’s solution to solve

utt=c2uxx,u(x,0)=f(x),ut(x,0)=g(x),0x.

Solution

The general solution of the wave equation was found in the form

u(x,t)=F(x+ct)+G(xct).

However, for this problem we can only obtain information for values of x and t such that 0x+ct and 0xct. In Figure 2.7.7 the characteristics x=ξ+ct and x=ηct for 0ξ, η. The main (gray) triangle, which is the domain of dependence of the point (,2,/2c), is the only region in which the solution can be found based solely on the initial conditions. As with the previous problem, boundary conditions will need to be given in order to extend the domain of the solution.

clipboard_ebc496012dfe0fc5fb76bacf4746216e2.png
Figure 2.7.7: The characteristics emanating from the interval 0x for the finite string problem.

In the last example we saw that a fixed boundary at x=0 could be satisfied when f(x) and g(x) are extended as odd functions. In Figure 2.7.8 we indicate how the characteristics are affected by drawing in the new one as red dashed lines. This allows us to now construct solutions based on the initial conditions under the line x=ct for 0x. The new region for which we can construct solutions from the initial conditions is indicated in gray in Figure 2.7.8.

clipboard_e4d58a1a25a2b0f72ffdf0f3ee354ff01.png
Figure 2.7.8: The red dashed lines are the characteristics from the interval [,0] from using the odd extension about x=0.

We can add characteristics on the right by adding a boundary condition at x=. Again, we could use fixed u(,t)=0, or free, ux(,t)=0, boundary conditions. This allows us to now construct solutions based on the initial conditions for x2.

Let’s consider a fixed boundary condition at x=. Then, the solution must satisfy

u(,t)=F(+ct)+G(ct)=0.

To see what this means, let ζ=+ct. Then, this condition gives (since ct=ζ)

F(ξ)=G(2ξ),ξ2.

Note that G(2ζ) is defined for 02ζ. Therefore, this is a well-defined extension of the domain of F(x).

Note that

F(ξ)=12f(ξ)+12c0g(s)ds.G(2ξ)=12f(2ξ)+12c2ξ0g(s)ds=12f(2ξ)12cξ0g(2σ)dσ

Comparing the expressions for G(ζ) and G(2ζ), we see that

f(ξ)=f(2ξ),g(ξ)=g(2ξ).

These relations imply that we can extend the functions into the region x> if we consider an odd extension of f(x) and g(x) about x=. This will give the blue dashed characteristics in Figure 2.7.9 and a larger gray region to construct the solution.

clipboard_e9b622fe3cd5a60260f2b5d79e07e1b33.png
Figure 2.7.9: The red dashed lines are the characteristics from the interval [,0] from using the odd extension about x=0 and the blue dashed lines are the characteristics from the interval [,2] from using the odd extension about x=.

So far we have extended f(x) and g(x) to the interval x2 in order to determine the solution over a larger xt-domain. For example, the function f(x) has been extended to

fext(x)={f(x),<x<0,f(x),0<x<,f(2x),<x<2.

A similar extension is needed for g(x). Inserting these extended functions into d’Alembert’s solution, we can determine u(x,t) in the region indicated in Figure 2.7.9.

Even though the original region has been expanded, we have not determined how to find the solution throughout the entire strip, [0,]×[0,). This is accomplished by periodically repeating these extended functions with period 2. This can be shown from the two conditions

f(x)=f(x),x0,

f(x)=f(2x),x2.

Now, consider

f(x+2)=f(2(x2))=f(x)=f(x).

This shows that f(x) is periodic with period 2. Since g(x) satisfies the same conditions, then it is as well.

In Figure 2.7.10 we show how the characteristics are extended throughout the domain strip using the periodicity of the extended initial conditions. The characteristics from the interval endpoints zig zag throughout the domain, filling it up. In the next example we show how to construct the odd periodic extension of a specific function.

clipboard_ef5022e1ce02caa92cc328f6a2208bda2.png
Figure 2.7.10: Extending the characteristics throughout the domain strip.

Example 2.7.4

Construct the periodic extension of the plucked string initial profile given by

f(x)={x,0x2,x,2x,

satisfying fixed boundary conditions at x=0 and x=.

Solution

We first take the solution and add the odd extension about x=0. Then we add an extension beyond x=. This process is shown in Figure 2.7.11.

clipboard_ec069df4d294f3e481255cd196fa6ec09.png
Figure 2.7.11: Construction of odd periodic extension for (a) The initial profile, f(x). (b) Make f(x) an odd function on [,]. (c) Make the odd function periodic with period 2.

We can use the odd periodic function to construct solutions. In this case we use the result from the last example for obtaining the solution of the problem in which the initial velocity is zero, u(x,t)=12[f(x+ct)+f(xct)]. Translations of the odd periodic extension are shown in Figure 2.7.12.

clipboard_e29b0d195dc8edea174a458925010de0f.png
Figure 2.7.12: Translations of the odd periodic extension.

In Figure 2.7.13 we show superimposed plots of f(x+ct) and f(xct) for different values of ct. A box is shown inside which the physical wave can be constructed. The solution is an average of these odd periodic extensions within this box. This is displayed in Figure 2.7.14.

clipboard_e1a8e0b6966b746dc69898d23b4a4ea21.png
Figure 2.7.13: Superimposed translations of the odd periodic extension.
clipboard_e3621b77251a0b32e30d8aba774f6dd32.png
Figure 2.7.14: On the left is a plot of f(x+ct), f(xct) from Figure 2.7.13 and the average, u(x,t). On the right the solution alone is shown for ct=0 to ct=1.

This page titled 2.7: d’Alembert’s Solution of the Wave Equation is shared under a CC BY-NC-SA 3.0 license and was authored, remixed, and/or curated by Russell Herman via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?