For a schemethe Fourier ansatz gives the amplification polynomialIts roots are the amplification factors of that Fourier mode.
Centered space and leapfrog time discretization of givesIts amplification polynomial is . The two roots have unit modulus for , but a repeated unit root at equality violates the uniform root condition.
Applying a centered leapfrog step in time and the centered second difference in space to givesFor every , each nonconstant Fourier mode has an amplification root of modulus greater than one, so the scheme is unconditionally unstable.
Articles by others on the same topic
There are currently no matching articles.