Module 1: Introduction to Finite Difference Method and Fundamentals of CFD
  Lecture 4:
 

Consistency

A finite difference representation of a partial differential equation (PDE) is said to be consistent if we can show that the difference between the PDE and its finite difference (FDE) representation vanishes as the mesh is refined, i.e,

  lim (PDE-FDE)  =    lim (TE)=0

A questionable scheme would be one for which the truncation error is and not explicitly or or higher orders.

In such cases the scheme would not be formally consistent unless the mesh were refined in a manner such that . Let us take Eq. (3.3) and use the Dufort-Frankel (1953) differencing scheme. The FDE is

(4.1)

Now the leading terms of truncated series form the truncation error for the complete equation:

 

The above expression for truncation for error is meaningful if together with and .However, and may individually approach zero in such a way that . Then if we reconstitute the PDE from FDE and TE, we shall obtain

     lim (PDE-FDE)   =     lim (TE)  =

And finally PDE becomes

 

We started with a parabolic one and ended with a hyperbolic one!

So, DuFort-Frankel scheme is not consistent for the 1D unsteady state heat conduction equation unless together with and .