The artifacts can be explained by some theory

Exact solution of the scheme:

$$ u^n = A^n,\quad A = \frac{1 - (1-\theta) a\Delta t}{1 + \theta a\Delta t}\thinspace .$$
Concluding remarks:

Only the Backward Euler scheme is guaranteed to always give qualitatively correct results.

\
« Previous