Tables 6.4.5, 6.4.6 and 6.4.7 are the numerical solutions respectively with r=0.1, 0.5 and 1.0. Unlike FTCS scheme the solution doesn't oscillate even with r=1. The solutions with r=1 at various t and x is compared in the table 6.4.8 with analytical solution. It is very clear from this table that the numerical and analytical solutions are in good agreement. However the only disadvantage with Crank-Nicolson scheme over FTCS scheme is the computational effort involved at each time level which is relatively high.
Weighted average scheme
The weighted average of schemes(6.4.2) and (6.4.5) given by
................. (6.4.8)
for i = 1, 2, ..............N-1 and n = 0, 1,2, 3.............
is called the weighted average scheme for one dimensional heat equation. Here, q =0, 1/2 and 1 gives FTCS, Crank - Nicolson and fully implicit schemes respectively.
Three
level schemes
All the schemes discussed until now are two level since only two time levels that is nth and (n+1)th are involed. If we use central differences to approximate both time and space derivative terms we get
................. (6.4.9)
for i = 1, 2, 3 ...........N-1 and n = 1, 2............
For n = 1 the scheme needs the solution at level which is outside the domain of interest and we may not have any information about u at this time level. Hence at n=1 it
is necessary to use a different scheme and from n = 2 onwards the schems(6.4.9) can be used. If we replace term in (6.4.9) with
Then we get
................. (6.4.10)
This difference scheme is called Dufort and Frankel explicit
method.