To obtain higher accuracy and good stability, the semi-implicit or
CTCS method is very useful, which combines both explicit and implicit
schemes. We start by writing this method in a generic, abstract form:
![$\displaystyle \mbf{A}(t_{n+1}) = \mbf{A}(t_n) + \Delta t [\mbf{DA}(t_{n+1/2}) + \mbf{U}[\mbf{A}(t_{n+1/2})] + O(\Delta t^3) + O(\Delta t \Delta x^2).$](img317.png) |
(B.1) |
From the Taylor expansion, the error due to the time-discretisation,
for making each step in time is
. Hence, the
accumulated error in reaching a given time
is
approximately proportional to
in the CTCS methods. Thus, a much lower error can be achieved for the
same time-step, while at the same time there is usually an increase in
stability.
Paul Cochrane
2002-04-18