4 Stability

There is a stability constraint that is common to many methods for obtaining numerical solutions of the wave equation. Issues relating to stability of numerical methods can be extremely complicated, but the following Key Point is enough for our purposes.

Key Point 25

The numerical method seen in this Section requires that

[maths rendering]

for solutions not to grow unrealistically with [maths rendering] .

This is called the CFL condition (named after an acronym of three mathematicians Courant, Friedrichs and Lewy).

Exercises
  1. Suppose that [maths rendering] satisfies the wave equation [maths rendering] in [maths rendering] and [maths rendering] . It is given that [maths rendering] satisfies boundary conditions [maths rendering] [maths rendering] and initial conditions that need not be stated for the purposes of this question. The application is such that the wave speed [maths rendering] .

    The numerical method

    [maths rendering]

    where [maths rendering] , is implemented using [maths rendering] and [maths rendering] .

    Suppose that, after 7 time-steps, the following data forms part of the numerical solution:

    [maths rendering]

    Carry out the next time-step so as to find an approximation to [maths rendering] at [maths rendering] .

  2. Suppose that [maths rendering] satisfies the wave equation [maths rendering] in [maths rendering] and [maths rendering] . It is given that [maths rendering] satisfies boundary conditions [maths rendering] [maths rendering] . The initial elevation may be summarised as

    [maths rendering]

    and the string is initially at rest (that is, [maths rendering] ). The application is such that the wave speed [maths rendering] .

    Carry out the first two time-steps of the numerical method

    [maths rendering]

    where [maths rendering] in which [maths rendering] and [maths rendering] .

  1. In this case [maths rendering] and the required time-step is carried out as follows: [maths rendering]

    to 4 decimal places and these are the approximations to [maths rendering] , [maths rendering] , [maths rendering] , [maths rendering] and [maths rendering] , respectively.

  2. In this case [maths rendering] and the first time-step is carried out as follows: [maths rendering]

    to 4 decimal places.

    The second time-step is as follows:

    [maths rendering]

    to 4 decimal places.