2.5.5 Linear stability analysis of naive FTCS scheme

Applying ansatz (123) to the naive scheme Eq. (92) gives

$$A e^{- ji k \Delta x} = e^{- ji k \Delta x} - \frac{\D... ... e^{- j\left ( i - 1 \right) k \Delta x} \right) \enspace .$$ (124)

Multiplying with $e^{ji k \Delta x}$ and using the Courant number

$$\alpha := \frac{\Delta t}{\Delta x} v_{}$$ (125)

we get

$$A = 1 - \alpha \frac{1}{2} \left( e^{- jk \Delta x} ... ...Delta x} \right) = 1 + j\alpha \sin k \Delta x \enspace ,$$ (126)

$$\mbox{abs} \! \left( A \right) = \left( 1 + \alpha^2 \sin... ...( 1 + \alpha^2 \sin^2 \pi/N \right)^{-\frac{1}{2}} \enspace ,$$ (127)

$${\color{HIGH1color} \mbox{abs} \! \left( A \right) > 1 ... ...for} \enspace N > 1 \enspace , \enspace \alpha > 0 \enspace .$$ (128)

All waves escept the ones with smallest wavenumber grow exponentially in time:

The scheme is unconditionally unstable, independent of the time-step.