4.2 Case 1

\(\tau =0.0001\) sec, \(h=0.1\ ft\)

Speed Method CPU time (sec) RMSE Animation (2D) plots
U=2 Explicit FTCS 20 0.0546 HTML HTML
Explicit LAX 31 0.0543 HTML
Implicit FTCS 45 0.0548 HTML
C-R 49 0.0544 HTML
U=t/20 Explicit FTCS 21 0.003 HTML HTML
Explicit LAX 31 0.0031 HTML
Implicit FTCS 67 0.0031 HTML
C-R 69 0.0032 HTML

Note the following: The explicit FTCS remained stable throughout the run due to the small time step. All other methods were stable as well during the run. For the CPU for the varying \(u\) case, notice that for the implicit methods this value is larger than the CPU for the same methods but when \(u\) is fixed. This is due to the fact that the matrix \(A\) is no longer constant, and must be recomputed at each time step before calling Thomas algorithm to solve \(Ax=b\) system.

Also notice that the CPU time for the implicit methods is larger than the explicit methods. This is due to the extra computational cost in solving \(Ax=b\). Even when using Thomas algorithm, this is still more expensive than the explicit methods when number of time steps is large.