2.6.2 The CFL Condition
Measurable Outcome 2.2, Measurable Outcome 2.5, Measurable Outcome 2.6, Measurable Outcome 2.7
Recall from 2.2.3 that the domain of dependence for the convection equation at \((x,t)\) is the characteristic \(x(s<t)\). We can also consider the numerical domain of dependence of the solution at \((x_ i, t^ n)\). Figure 2.18 shows the numerical domain of dependence of the FTBS method. Note that the numerical domain of dependence at \(n-2\) includes the nodes \(i-2, i-1\), and \(i\). \((i-2,n-2)\) and (\(i-1,n-2\)) are in the domain of dependence of \((i,n)\) because these values are required in the calculation of \(U_{i-1}^{n-1}\), which in turn is required to calculate \(U_{i}^ n\). Similarly, \((i-1,n-2)\) and \((i,n-2)\) are in the domain of dependence of \((i,n)\) because they are required in the calculation of \(U_{i}^{n-1}\).
We now define two terms:
CFL condition: The discretization of a PDE satisfies the CFL condition if the numerical domain of dependence includes the physical domain of dependence.
CFL Theorem: The CFL condition is a necessary condition for the discretization of a time-dependent PDE to be convergent (i.e. in the limit as \(\Delta t \to 0\)). Courant, Friedrichs, and Levy are the authors who first described this requirement in 1928 (well before the first computer!)
In Figures 2.19 and 2.20, the CFL conditions are shown graphically for the one-dimensional convection equation discretized using the first-order upwind scheme.
For scheme, the CFL condition requires that
If this condition is not satisfied, then the characteristic will be outside of the numerical domain of dependence and the numerical solution cannot possibly converge to the exact solution if the same ration of \({\Delta t}/\Delta x\) is maintained as \({\Delta t}\to 0\), since the numerical method will not have the proper domain of dependence in that limit. The non-dimensional number \(\frac{|u|{\Delta t}}{\Delta x}\) is called the CFL Number or just the CFL. In general, the CFL condition for explicit finite difference methods for convection will require that the CFL number be bounded by a constant which will depend upon the particular numerical scheme (for FTBS, the constant is \(1\)).
Exercise
Suppose we wish to solve the 1-D convection equation with velocity \(u=2\) on a mesh with \(\Delta x = 0.1\) using the FTBS method. What is the largest possible time step for which the scheme satisfies the CFL condition?
We can solve for the maximum \(\Delta t\) as \(\Delta x/ |u| = 0.05.\)