In mathematics, the Courant–Friedrichs–Lewy (CFL) condition is a necessary condition for convergence while solving certain partial differential equations (usually hyperbolic PDEs) numerically by the method of finite differences. It arises in the numerical analysis of explicit time integration schemes, when these are used for the numerical solution. As a consequence, the time step must be less than a certain time in many explicit time-marching computer simulations, otherwise the simulation will produce incorrect results. The condition is named after Richard Courant, Kurt Friedrichs, and Hans Lewy who described it in their 1928 paper.
The principle behind the condition is that, for example, if a wave is moving across a discrete spatial grid and we want to compute its amplitude at discrete time steps of equal duration, then this duration must be less than the time for the wave to travel to adjacent grid points. As a corollary, when the grid point separation is reduced, the upper limit for the time step also decreases. In essence, the numerical domain of dependence of any point in space and time (as determined by initial conditions and the parameters of the approximation scheme) must include the analytical domain of dependence (wherein the initial conditions have an effect on the exact value of the solution at that point) in order to assure that the scheme can access the information required to form the solution.
The CFL condition
In order to make a reasonably formally precise statement of the condition, it is necessary to define the following quantities
- Spatial coordinate: it is one of the coordinates of the physical space in which the problem is posed.
- Spatial dimension of the problem: it is the number of spatial dimensions i.e. the number of spatial coordinates of the physical space where the problem is posed. Typical values are , and .
- Time: it is the coordinate, acting as a parameter, which describes the evolution of the system, distinct from the spatial coordinates.
The spatial coordinates and the time are supposed to be discrete-valued independent variables, which are placed at regular distances called the interval length and the time step, respectively. Using these names, the CFL condition relates the length of the time step to a function of the interval lengths of each spatial coordinate and of the maximum speed with which information can travel in the physical space.
The one-dimensional case
For one-dimensional case, the CFL has the following form:
where the dimensionless number is called the Courant number,
- is the magnitude of the velocity (whose dimension is length/time)
- is the time step (whose dimension is time)
- is the length interval (whose dimension is length).
The value of changes with the method used to solve the discretised equation, especially depending on whether the method is explicit or implicit. If an explicit (time-marching) solver is used then typically . Implicit (matrix) solvers are usually less sensitive to numerical instability and so larger values of may be tolerated.
The two and general n-dimensional case
In the two-dimensional case, the CFL condition becomes
with obvious meaning of the symbols involved. By analogy with the two-dimensional case, the general CFL condition for the -dimensional case is the following one:
The interval length is not required to be the same for each spatial variable . This "degree of freedom" can be used in order to somewhat optimize the value of the time step for a particular problem, by varying the values of the different interval in order to keep it not too small.
- See reference Courant, Friedrichs & Lewy 1928. There exists also an English translation of the 1928 German original: see references Courant, Friedrichs & Lewy 1956 and Courant, Friedrichs & Lewy 1967.
- This situation commonly occurs when a hyperbolic partial differential operator has been approximated by a finite difference equation, which is then solved by numerical linear algebra methods.
- This quantity is not necessarily the same for each spatial variable, as it is shown in "The two and general n–dimensional case" section of this entry : it can be chosen in order to somewhat relax the condition.
- Precisely, this is the hyperbolic part of the PDE under analysis.
- Courant, R.; Friedrichs, K.; Lewy, H. (1928), "Über die partiellen Differenzengleichungen der mathematischen Physik", Mathematische Annalen (in German), 100 (1): 32–74, Bibcode:1928MatAn.100...32C, JFM 54.0486.01, MR 1512478, doi:10.1007/BF01448839.
- Courant, R.; Friedrichs, K.; Lewy, H. (September 1956) , On the partial difference equations of mathematical physics, AEC Research and Development Report, NYO-7689, New York: AEC Computing and Applied Mathematics Centre – Courant Institute of Mathematical Sciences, pp. V + 76, archived from the original on October 23, 2008.: translated from the German by Phyllis Fox. This is an earlier version of the paper Courant, Friedrichs & Lewy 1967, circulated as a research report.
- Courant, R.; Friedrichs, K.; Lewy, H. (March 1967) , "On the partial difference equations of mathematical physics", IBM Journal of Research and Development, 11 (2): 215–234, Bibcode:1967IBMJ...11..215C, MR 0213764, Zbl 0145.40402, doi:10.1147/rd.112.0215. A freely downloadable copy can be found here.