Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random |  
Login
create new user
Username:
Password:
forget your password?
Main Menu
Sections

Meta

Talkback

Downloads

Information
[parent] example of Wave Mechanics: Initial Conditions and the d'Alembert Solution (Example)

Wave Mechanics Examples: Initial Conditions and the d’Alembert Solution

This companion article provides exercises for WM17, wave mechanics: Initial Conditions and the d’Alembert Solution. All exercises are stated first so they can be attempted without seeing the answers. Complete worked solutions follow in Part II.

For the one-dimensional constant-speed wave equation on the whole line,

|-----------|
|utt = c2uxx ,
------------
(1)

with initial data

u(x,0) = f(x )
(2)

and

ut(x,0) = g(x ),
(3)

the d’Alembert formula is

|--------------------------------------∫-x+ct--------|
|         1-                         1--             |
|u(x,t) = 2 [f (x − ct) + f(x + ct)] + 2c x− ct g (s) ds.|
------------------------------------------------------
(4)

The first term carries the initial displacement along the two characteristic directions. The integral term carries the effect of the initial velocity over the characteristic interval. This is the classical whole-line initial-value solution under the usual smoothness assumptions [436].

How to use this problem set

Attempt all exercises in Part I before consulting Part II. In each problem, identify the characteristic endpoints

x − ct    and     x + ct
(5)

before substituting into the formula. Keep the whole-line assumption in mind: a finite string also requires endpoint boundary conditions and generally involves reflections or mode expansions.

Part I: Exercises

Exercise 1: Read the d’Alembert formula

Consider

                                      ∫ x+ct
         1-                         1--
u(x,t) = 2 [f (x − ct) + f(x + ct)] + 2c x− ct g (s) ds.
(6)

Answer the following.

  1. Which function represents the initial displacement?
  2. Which function represents the initial velocity?
  3. What are the two characteristic footpoints on the line t = 0?
  4. Which part of the formula survives if g(x) = 0?
  5. Which part survives if f(x) = 0?

Exercise 2: Released displacement with zero initial velocity

Suppose

g(x) = 0.
(7)

  1. Simplify the d’Alembert formula.
  2. Explain why the initial profile splits into two copies.
  3. What is the amplitude of each copy relative to the original amplitude?
  4. What direction does each copy move?
  5. Show explicitly that the two copies add back to f(x) at t = 0.

The geometry is illustrated below.

PIC

Figure. With zero initial velocity, the initial displacement separates into equal right- and left-moving half-amplitude copies.

Exercise 3: Initial velocity required for one-way motion

An initial shape f(x) is prescribed.

  1. What initial velocity g(x) produces a pure right-moving solution?
  2. What initial velocity produces a pure left-moving solution?
  3. For
    f(x) = A cos(kx ),
    (8)

    find g(x) for pure right-moving motion.

  4. For the same f(x), find g(x) for pure left-moving motion.

PIC

Figure. The initial displacement alone does not determine propagation direction. The initial velocity carries the directional information.

Exercise 4: Recover the two traveling components

Define

        ∫  x
H (x) =     g(s)ds,
         x∗
(9)

where x is a fixed reference point.

Starting from

u(x,t) = F (x − ct) + G (x + ct),
(10)

show that one convenient choice of the component functions is

        1        1
F(x) =  -f(x) − ---H (x)
        2       2c
(11)

and

        1        1
G (x) = --f(x) + --H (x).
        2        2c
(12)

Explain why adding a constant to F and subtracting the same constant from G does not change u(x,t).

Exercise 5: Direct numerical evaluation

Let

c = 2.5 m/s,
(13)

and let the initial data be

f(x ) = αx2,     α = 0.020m − 1,
(14)

g(x) = βx,     β = 0.30 s−1.
(15)

Find

u(1.0 m,0.20 s).
(16)

Show separately the displacement contribution and the initial-velocity contribution.

Exercise 6: Velocity-only compact pulse

Suppose the initial displacement is zero,

f(x ) = 0,
(17)

and the initial velocity is

       {
         V0,  |x | ≤ a,
g(x) =
         0,   |x | > a.
(18)

At the observation point x = 0:

  1. Find u(0,t) for 0 t a∕c.
  2. Find u(0,t) for t a∕c.
  3. Evaluate the result for
    V0 = 0.40m/s,      a = 0.60m,     c = 3.0 m/s
    (19)

    at t = 0.10 s and t = 0.30 s.

The integral is controlled by the overlap between the characteristic interval and the support of g.

PIC

Figure. Only the portion of the initial-velocity profile inside [x ct,x + ct] contributes to the velocity integral.

Exercise 7: Sinusoidal initial displacement

Let

f(x ) = A cos(kx ),    g(x ) = 0.
(20)

  1. Use the d’Alembert formula to find u(x,t).
  2. Use a trigonometric identity to simplify the result.
  3. Show that the angular frequency is ω = ck.
  4. Explain why the final expression is a standing-wave form even though d’Alembert’s solution is built from traveling waves.

Exercise 8: Sinusoidal initial velocity

Let

f(x) = 0,     g(x) = V0 sin(kx ).
(21)

  1. Evaluate the velocity integral.
  2. Show that
             V0
u(x,t) = ck-sin (kx)sin(kct).
    (22)

  3. Rewrite the coefficient using ω = ck.
  4. Differentiate the result with respect to time and verify ut(x, 0) = g(x).

Exercise 9: Domain of dependence

At the observation point

x0 =  4.0 m,     t0 = 0.30s,
(23)

suppose the wave speed is

c = 5.0 m/s.
(24)

  1. Find the two characteristic footpoints on the initial line.
  2. State the interval of initial data that can influence (x0,t0).
  3. Can a change in the initial data at x = 2.2 m influence the observation point?
  4. Can a change at x = 5.4 m influence it?
  5. Can a change at x = 6.0 m influence it?

PIC

Figure. The point (x0,t0) depends only on initial data inside the interval [x0 ct0,x0 +ct0].

Exercise 10: Finite propagation from compactly supported data

Suppose both initial functions vanish outside

− 1.0 m ≤ x ≤ 1.0m.
(25)

The wave speed is

c = 3.0 m/s.
(26)

At time

t = 0.40 s,
(27)

find the largest interval in which the solution can possibly be nonzero. Explain how this demonstrates finite propagation speed.

Exercise 11: Check the initial conditions directly

Starting from

         1                          1 ∫ x+ct
u(x,t) = --[f(x − ct) + f(x + ct)] +---      g(s) ds,
         2                         2c  x−ct
(28)

show directly that

u(x,0) = f(x )
(29)

and

ut(x,0) = g(x ).
(30)

For the velocity check, use the Leibniz rule for differentiating an integral whose limits depend on time.

Exercise 12: Whole-line formula versus a finite string

A string occupies only

0 ≤ x ≤ L
(31)

and is fixed at both ends:

u (0, t) = 0,     u(L, t) = 0.
(32)

Answer the following.

  1. Does the basic whole-line d’Alembert formula automatically enforce these endpoint conditions?
  2. What additional physical phenomenon appears when waves reach the endpoints?
  3. Name two methods that can be used to solve the finite-string problem after boundary conditions are imposed.
  4. Explain why applying the whole-line formula blindly to a finite string can give the wrong physical answer.

Exercise 13: Connect the initial-value solution to string mechanics

An ideal string has

T =  100N,      μ = 0.025kg/m.
(33)

The initial displacement is

                [   (  x )2 ]
f(x) = 0.008 exp   −  ----     m
                     0.20
(34)

and the initial velocity is zero.

  1. Find the wave speed.
  2. Write the complete d’Alembert solution.
  3. At t = 0.010 s, where are the centers of the two Gaussian pulses?
  4. What is the peak amplitude of each separated pulse?

Exercise 14: Synthesis – design a pure right-moving sinusoid

An ideal string has

T =  160N,      μ = 0.100kg/m.
(35)

The desired initial displacement is

f (x ) = 0.006 cos(5x) m.
(36)

The goal is for this entire shape to travel purely to the right.

  1. Find the wave speed c.
  2. Find the required initial velocity g(x).
  3. Write the resulting traveling-wave solution u(x,t).
  4. Identify k and ω.
  5. Find the wavelength and ordinary frequency.
  6. Explain why choosing g(x) = 0 instead would not produce the same motion.

Part II: Complete Worked Solutions

Solution 1: Read the d’Alembert formula

The formula is

                                      ∫ x+ct
u(x,t) = 1-[f (x − ct) + f(x + ct)] + 1--     g (s) ds.
         2                          2c x− ct
(37)

  1. The initial displacement is
    |--------------|
f (x) = u(x,0).|
----------------
    (38)

  2. The initial velocity is
    |---------------|
g(x) = u (x,0 ). |
---------t-------
    (39)

  3. The characteristic footpoints are
    |--------------------|
x-−--ct---and--x-+-ct.-
    (40)

  4. If g = 0, only the displaced copies of f remain:
    |----------------------------------|
|u(x,t) = 1-[f (x − ct) + f(x + ct)].|
----------2------------------------|
    (41)

  5. If f = 0, only the velocity integral remains:
    |--------------------------|
|          1 ∫ x+ct        |
|u(x,t) = ---      g(s)ds. |
----------2c--x−ct---------
    (42)

Solution 2: Released displacement with zero initial velocity

Setting

g (x ) = 0
(43)

removes the integral term, so

|----------------------------------|
|         1-           1-          |
|u(x,t) = 2 f(x − ct) + 2 f(x + ct).
-----------------------------------
(44)

The first term is a right-moving copy because its argument is x ct. The second is a left-moving copy because its argument is x + ct.

Each copy has one-half the initial amplitude. At t = 0,

u(x, 0) = 1
--
2f(x) + 1
--
2f(x) (45)
= f(x) . (46)

Thus the two half-amplitude copies exactly reconstruct the original profile at the initial instant.

Solution 3: Initial velocity required for one-way motion

For a pure right-moving wave,

u (x, t) = f (x − ct).
(47)

Differentiate with respect to time:

              ′
ut(x,t) = − cf (x − ct).
(48)

At t = 0,

|----------------|
|g(x) = − cf′(x).|
-----------------
(49)

Similarly, for a pure left-moving wave,

u (x, t) = f (x + ct),
(50)

so

|----------------|
|g(x) = +cf ′(x).|
-----------------
(51)

Now let

f(x) = A cos(kx ).
(52)

Then

f ′(x) = − Ak sin(kx).
(53)

Therefore the required right-moving initial velocity is

------------------------
|                      |
gright(x)-=-cAk-sin(kx),-
(54)

whereas the required left-moving initial velocity is

|------------------------|
-gleft(x-) =-− cAk-sin(kx-).|
(55)

The sign of the initial velocity therefore selects the propagation direction.

Solution 4: Recover the two traveling components

WM17 gives

  ′      1- ′     -1-
F  (x) = 2f (x) − 2c g(x)
(56)

and

  ′      1-′      1--
G (x ) = 2f (x) + 2cg(x ).
(57)

Define

        ∫  x
H (x) =     g(s)ds.
         x∗
(58)

Integrating gives one convenient choice,

|------------------------|
|        1-      -1-     |
-F(x)-=--2f(x)-−-2c-H-(x)-
(59)

and

|-------1--------1--------|
G (x) = --f(x) + --H (x). |
--------2--------2c-------|
(60)

Their sum is

F (x) + G (x) = f(x),
(61)

as required by the initial displacement.

If a constant C is added to F and subtracted from G, then

[F  + C ] + [G − C ] = F + G.
(62)

Therefore the physical field u(x,t) is unchanged by this constant redistribution between the two component functions.

Solution 5: Direct numerical evaluation

The observation event is

x =  1.0 m,     t = 0.20s,
(63)

with

c = 2.5 m/s.
(64)

The characteristic endpoints are

x ct = 1.0 (2.5)(0.20) = 0.50 m, (65)
x + ct = 1.0 + (2.5)(0.20) = 1.50 m. (66)

The displacement contribution is

1
--
2[f(0.50) + f(1.50)] = 1
--
2α[(0.50)2 + (1.50)2] (67)
= 1-
2(0.020)(0.25 + 2.25) m (68)
= 0.025 m . (69)

The velocity contribution is

1
---
2c 0.501.50βsds =  β
---
2c[s2 ]
 --
  20.501.50 (70)
= 0.30
-5.0[           ]
 2.25 − 0.25
 -----2----- m (71)
= 0.060 m . (72)

Therefore

|--------------------------|
-u(1.0m,-0.20-s)-=-0.085-m.-|
(73)

Solution 6: Velocity-only compact pulse

Because f = 0,

            ∫ ct
u(0,t) = -1-    g(s)ds.
         2c  −ct
(74)

For early times with

ct ≤ a,
(75)

the entire interval [ct,ct] lies inside the region where g = V 0. Thus

u(0,t) = 1--
2cV 0(2ct) (76)
= V 0t , 0 t                                                                          a-
                                                                         c. (77)

For later times,

ct ≥ a,
(78)

the integration interval already covers the full support [a,a]. Hence

u(0,t) = -1-
2cV 0(2a) (79)
= V0a-
 c , t                                                                         a-
                                                                        c. (80)

With

a-=  0.60-=  0.20 s,
c    3.0
(81)

we have at t = 0.10 s,

                                |--------|
u(0,0.10) = V0t = (0.40 )(0.10) = -0.040m--.
(82)

At t = 0.30 s,

             (0.40-)(0.60)-  |--------|
u (0,0.30 ) =     3.0     = -0.080-m-.
(83)

The plateau occurs because, after t = a∕c, expanding the characteristic interval does not include any new region where g is nonzero.

Solution 7: Sinusoidal initial displacement

With

f(x ) = A cos(kx ),    g(x ) = 0,
(84)

the d’Alembert formula becomes

         A-                A-
u(x,t) = 2 cos[k(x − ct)] + 2 cos[k (x +  ct)].
(85)

Use

cos(α − β ) + cos(α + β) = 2cos αcos β.
(86)

With

α =  kx,     β = kct,
(87)

we obtain

|----------------------------|
|u(x,t) = A cos(kx )cos(kct).|
-----------------------------
(88)

The temporal factor has angular frequency

|--------|
-ω-=-ck.-|
(89)

The result is a standing-wave form because equal right- and left-moving waves of the same amplitude superpose to make a nontranslating pattern with fixed nodes and antinodes.

Solution 8: Sinusoidal initial velocity

Here

f(x) = 0,     g(x) = V0 sin(kx ).
(90)

Therefore

             ∫
          V0-  x+ct
u (x,t) = 2c       sin (ks )ds.
              x−ct
(91)

Integrate:

u(x,t) = V0
----
2ck[cos(kx − kct) − cos(kx + kct)] . (92)

Using

cos(α −  β) − cos(α + β) = 2 sin α sin β,
(93)

we obtain

|----------------------------|
|         V0-                |
|u(x,t) = ck sin (kx)sin(kct).|
-----------------------------
(94)

Since

ω = ck,
(95)

this can be written as

|----------------------------|
|u(x, t) = V0-sin(kx) sin(ωt ). |
-----------ω-----------------|
(96)

Differentiate with respect to time:

ut(x,t) = V0sin(kx) cos(ωt).
(97)

At t = 0,

|----------------------------|
-ut(x,0)-=-V0-sin(kx-) =-g(x).|
(98)

Solution 9: Domain of dependence

The characteristic endpoints are

x0 ct0 = 4.0 (5.0)(0.30) = 2.5 m, (99)
x0 + ct0 = 4.0 + (5.0)(0.30) = 5.5 m. (100)

Therefore the domain of dependence on the initial line is

|------------------|
|2.5 m ≤  x ≤ 5.5m. |
--------------------
(101)

A change at x = 2.2 m lies outside this interval, so it cannot influence (x0,t0).

A change at x = 5.4 m lies inside the interval, so it can influence the observation point.

A change at x = 6.0 m lies outside the interval, so it cannot influence the observation point.

Thus

|------------------------------------|
|2.2 : no,     5.4 : yes,    6.0 : no.|
--------------------------------------
(102)

Solution 10: Finite propagation from compactly supported data

Initially, both f and g vanish outside

[− 1.0,1.0] m.
(103)

In time t, information can propagate outward by at most the distance

ct = (3.0)(0.40) = 1.2m.
(104)

Therefore the largest possible support at that time is

[− 1.0 − 1.2, 1.0 + 1.2] m.
(105)

Hence

|--------------------|
-− 2.2-m-≤-x ≤-2.2m.--
(106)

Outside this interval, neither characteristic reaching (x,t) can intersect nonzero initial data. The disturbance therefore cannot appear arbitrarily far away instantly. This is finite propagation at speed c.

Solution 11: Check the initial conditions directly

Start with

                                      ∫
         1-                        -1-  x+ct
u(x,t) = 2 [f(x − ct) + f(x + ct)] + 2c      g(s) ds.
                                       x−ct
(107)

Set t = 0:

u(x, 0) = 1-
2[f(x) + f(x)] + 1--
2c xxg(s) ds (108)
= f(x). (109)

Thus

|--------------|
u (x,0) = f(x).|
----------------
(110)

Now differentiate with respect to time. The first two terms give

− cf ′(x − ct) + cf ′(x + ct).
  2             2
(111)

For the integral term, the Leibniz rule gives

∂
---
∂t[ 1 ∫ x+ct       ]
 ---      g(s)ds
 2c  x−ct = 1
---
2c[cg(x + ct) + cg (x − ct)] (112)
= 1-
2[g(x + ct) + g(x ct)]. (113)

Therefore

ut(x,t) = c
2-f(x ct) + c
2-f(x + ct) (114)
+ 1
--
2[g(x + ct) + g(x ct)]. (115)

At t = 0, the two fterms cancel:

u (x,0) = 1[g(x) + g(x)] = g(x).
 t        2
(116)

Hence

-----------------
|               |
ut(x,0)-=-g(x-).--
(117)

Solution 12: Whole-line formula versus a finite string

  1. No. The basic whole-line d’Alembert formula is constructed for an unbounded spatial domain. It does not automatically impose
    u (0, t) = 0,     u(L, t) = 0.
    (118)

  2. When a wave reaches a finite endpoint, it reflects. The boundary condition determines the reflected behavior.
  3. Two standard methods are reflected or extended initial data and normal-mode/Fourier-series expansion.
  4. Applying the whole-line formula blindly would allow the wave to pass beyond the physical endpoints instead of reflecting and satisfying the endpoint constraints. The resulting field could therefore violate the actual boundary conditions.

Thus the whole-line formula solves an initial-value problem, while a finite string requires an initial-boundary-value problem.

Solution 13: Connect the initial-value solution to string mechanics

For an ideal string,

    ∘ ---
       T-
c =    μ.
(119)

With

T =  100N,      μ = 0.025kg/m,
(120)

we obtain

c = ∘ ------
  -100--
  0.025 (121)
= √-----
 4000 (122)
63.2 m/s . (123)

Because g = 0, the solution is

u(x,t) = 0.004 exp [   (       )2]
      x −-ct
 −    0.20 (124)
+ 0.004 exp [             ]
    ( x + ct)2
 −    ------
      0.20 m. (125)

At

t = 0.010 s,
(126)

the pulse centers are at

x = ±ct = ± (63.2)(0.010)m.
(127)

Therefore

|--------------|
|x ≈ ±0.632 m. |
----------------
(128)

Each separated pulse has one-half the initial peak amplitude:

|--------------------------|
Apulse = 0.004m  = 4.0mm.  |
----------------------------
(129)

Solution 14: Synthesis – design a pure right-moving sinusoid

The wave speed is

c = ∘  ---
   T
   --
   μ (130)
= ∘  ------
    160
   ------
   0.100 (131)
= √ -----
  1600 (132)
= 40 m/s . (133)

The initial displacement is

f (x ) = 0.006 cos(5x) m.
(134)

Differentiate:

f ′(x) = − 0.030sin(5x).
(135)

For pure right-moving motion,

g(x) = − cf′(x).
(136)

Therefore

g(x) = 40[0.030 sin(5x)] (137)
= 1.20 sin(5x) m/s . (138)

The complete right-moving solution is

|--------------------------------|
|u(x,t) = 0.006cos[5(x − 40t)]m. |
---------------------------------
(139)

Expanding the phase,

5(x − 40t) = 5x − 200t.
(140)

Hence

|------------|    |--------------|
k = 5 rad/m  ,    |ω = 200 rad/s .
--------------    ---------------
(141)

The wavelength is

λ = 2π
---
k (142)
= 2π-
5 (143)
1.26 m . (144)

The ordinary frequency is

fwave = ω
---
2π (145)
= 200-
2π (146)
31.8 Hz . (147)

If instead g(x) = 0, the d’Alembert formula would produce

1f(x − ct) + 1f (x +  ct),
2            2
(148)

which is an equal two-way split, not a single right-moving waveform. The nonzero initial velocity is what removes the unwanted left-moving component.

Common mistakes

  • Mistake: forgetting that both f and g are required for a second-order-in-time wave equation.
  • Mistake: omitting the factor 12 on the two displaced copies of f.
  • Mistake: reversing the limits of the velocity integral.
  • Mistake: omitting the factor 1(2c) in front of the velocity integral.
  • Mistake: assuming g = 0 creates a single traveling copy. It creates two equal half-amplitude copies moving in opposite directions.
  • Mistake: confusing the domain of dependence with the entire initial line. Only [x ct,x + ct] can affect (x,t).
  • Mistake: using the whole-line d’Alembert formula on a finite string without enforcing boundary conditions.
  • Mistake: using the wrong sign in g = cffor one-way motion.

What WM17E1 reinforces

The d’Alembert formula converts physical initial data directly into a wave field:

|--------------------------------------∫-------------|
|         1                          1   x+ct        |
|u(x,t) = --[f(x − ct) + f(x + ct)] +---      g(s) ds.|
----------2-------------------------2c--x−ct---------
(149)

The exercises reinforce four structural ideas:

  • initial displacement and initial velocity play different roles,
  • the wave equation contains independent right- and left-moving information,
  • information propagates only within characteristic cones at speed c, and
  • pure one-way motion requires compatible initial displacement and velocity.

These ideas complete the basic whole-line initial-value problem for the one-dimensional ideal wave equation and prepare the way for more advanced topics such as energy transport, Fourier decompositions, modal expansions, interfaces, and dispersive waves.

References

References

[1]   A. P. French, Vibrations and Waves, M.I.T. Introductory Physics Series, W. W. Norton & Company, 1971.

[2]   Frank S. Crawford, Jr., Waves, Berkeley Physics Course, Volume 3, McGraw-Hill, 1968.

[3]   Richard P. Feynman, Robert B. Leighton, and Matthew Sands, The Feynman Lectures on Physics, Volume I, Chapter 47, “Sound. The Wave Equation.”

[4]   Walter A. Strauss, Partial Differential Equations: An Introduction, Second Edition, John Wiley & Sons, 2008.

[5]   Massachusetts Institute of Technology, 8.03SC Physics III: Vibrations and Waves, Lecture 10, “Traveling Waves,” MIT OpenCourseWare.

[6]   Gilbert Strang and Cleve Moler, Learn Differential Equations: Up Close, “Wave Equation,” MIT OpenCourseWare, 2015.


"example of Wave Mechanics: Initial Conditions and the d'Alembert Solution" is owned by bloftin.
(view preamble)
View style:
Other names:  WM17E1
Keywords:  wave mechanics, one-dimensional wave equation, d'Alembert solution, initial displacement, initial velocity, characteristics, domain of dependence, finite propagation speed, right-moving waves, left-moving waves, exercises, worked solutions

This object's parent.

Cross-references: energy, domain, field, speed, identity, motion, function, boundary, velocity, formula, wave equation, mechanics, wave, WM17

This is version 1 of example of Wave Mechanics: Initial Conditions and the d'Alembert Solution, born on 2026-09-12.
Object id is 1182, canonical name is ExampleOfWaveMechanicsInitialConditionsAndTheDAlembertSolution.
Accessed 5 times total.

Classification:
Physics Classification46.40.Cd (Mechanical wave propagation (including diffraction, scattering, and)
 46.40.-f (Vibrations and mechanical waves )
 02.30.Jr (Partial differential equations)
Pending Errata and Addenda
None.
Discussion
Style: Expand: Order:

No messages.

Interact
rate | post | correct | update request | add example | add (any)