In mathematics, the Courant–Friedrichs–Lewy condition (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 when explicit time-marching schemes 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 wildly incorrect results. The condition is named after Richard Courant, Kurt Friedrichs, and Hans Lewy who described it in their 1928 paper.
- 1 Heuristic description
- 2 The CFL condition
- 3 Implications of the CFL condition
- 4 See also
- 5 Notes
- 6 References
- 7 External links
The information 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 length, then this length must be less than the time for the wave to travel 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 (which data values in the initial conditions affect the numerical computed value at that point) must include the analytical domain of dependence (where in the initial conditions has 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 n of spatial dimensions i.e. the number of spatial coordinates of the physical space where the problem is posed. Typical values are n = 1, n = 2 and n = 3.
- 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, whose minimal steps are called respectively the interval length and the time step: the CFL condition relates the length of the time step to a function interval lengths of each spatial variable.
The one–dimensional case
For one-dimensional case, the CFL has the following form:
- u is the velocity (whose dimension is Length/Time)
- Δt is the time step (whose dimension is Time)
- Δx is the length interval (whose dimension is Length),
- C is a dimensionless constant which depends only on the particular equation to be solved.
The dimensionless number
is called the Courant number.
The two and general n–dimensional case
In the two–dimensional case, the CFL condition becomes
with obvious meaning of the symbols involved. In analogy with the two–dimensional case, the general CFL condition for the n–dimensional case is the following one
Note that the interval length it is not required to be the same for each spatial variable Δxi, i = 1 ,..., n . 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.
Implications of the CFL condition
The CFL condition is only a necessary one
As already remarked, the CFL condition is a necessary condition, but may not be sufficient for the convergence of the Finite-difference approximation of a given numerical problem. Thus, in order to establish the convergence of the finite-difference approximation, it is necessary to use other methods, which in turn could imply further limitations on the length of the time step and/or the lengths of the spatial intervals.
The CFL condition can be a very strong requirement
The CFL condition can be a very limiting constraint on the time step Δt: for example, in the finite-difference approximation of certain fourth-order nonlinear partial differential equations, it can have the following form
meaning that a decrease in the length interval Δx requires a fourth order decrease in the time step Δt for the condition to be fulfilled. Therefore, when solving particularly stiff problems, efforts are often made to avoid the CFL condition, for example by using implicit methods. However, in a recent work, a modern dynamical systems approaches to modeling, based upon center manifold theory, is demonstrated to provide theoretical support for the construction of non-traditional discretisations that automatically overcome the CFL restriction: see the article by Roberts (2003) for further information.
- Richard Courant
- Finite difference method
- Kurt Otto Friedrichs
- Implicit method
- Hans Lewy
- Numerical analysis
- ^ In general, it is not a sufficient condition: also, it can be also a demanding condition for some problems. See the "Implications of the CFL condition" section of this article for a brief survey of this issues.
- ^ 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 for each spatial variable, as it is shown in the "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.
- ^ Precisely, it does not depend on the time step Δt and on the length interval Δx.
- ^ See section 3 in the article by Roberts (2003)
- Courant, R.; Friedrichs, K.; Lewy, H. (1928), "Über die partiellen Differenzengleichungen der mathematischen Physik" (in German), Mathematische Annalen 100 (1): 32–74, Bibcode 1928MatAn.100...32C, doi:10.1007/BF01448839, JFM 54.0486.01, MR1512478, http://resolver.sub.uni-goettingen.de/purl?GDZPPN002272636 .
- Courant, R.; Friedrichs, K.; Lewy, H. (September 1956) , On the partial difference equations of mathematical physics, AEC Researh 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, http://www.archive.org/stream/onpartialdiffere00cour#page/n0/mode/2up .: 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, MR0213764, Zbl 0145.40402, http://domino.research.ibm.com/tchjr/journalindex.nsf/a3807c5b4823c53f85256561006324be/769774a3c9f3685f85256bfa00683f8a!OpenDocument . A freely downlodable copy can be found here.
- Roberts, A. J. (2003), "A holistic finite difference approach models linear dynamics consistently", Mathematics of Computation 72: 247–262, doi:10.1090/S0025-5718-02-01448-5, MR1933820, Zbl 1008.37047, http://www.ams.org/mcom/2003-72-241/S0025-5718-02-01448-5 .
Wikimedia Foundation. 2010.
Look at other dictionaries:
Richard Courant — Infobox Scientist box width = 300px name = Richard Courant image size = 300px caption = birth date = January 8, 1888 birth place = Lublinitz, Kingdom of Prussia death date = January 27 1972 death place = residence = citizenship = nationality =… … Wikipedia
Nombre de Courant — Le nombre de Courant Co est un nombre sans dimension utilisé en informatique et en mathématiques et plus particulièrement en calcul par différences finies. Ce nombre porte le nom de Richard Courant, un mathématicien allemand. Définition On le… … Wikipédia en Français
List of mathematics articles (C) — NOTOC C C closed subgroup C minimal theory C normal subgroup C number C semiring C space C symmetry C* algebra C0 semigroup CA group Cabal (set theory) Cabibbo Kobayashi Maskawa matrix Cabinet projection Cable knot Cabri Geometry Cabtaxi number… … Wikipedia
Критерий Куранта — Фридрихса — Леви — (критерий КФЛ) необходимое условие устойчивости явного численного решения некоторых дифференциальных уравнений в частных производных. Как следствие, во многих компьютерных симуляциях временной шаг должен быть меньше определённого значения, иначе… … Википедия
Критерий Куранта — Фридрихса Леви (критерий КФЛ) необходимое условие устойчивости явного численного решения некоторых дифференциальных уравнений в частных производных. Как следствие, во многих компьютерных симуляциях временной шаг должен быть меньше определённого… … Википедия
List of numerical analysis topics — This is a list of numerical analysis topics, by Wikipedia page. Contents 1 General 2 Error 3 Elementary and special functions 4 Numerical linear algebra … Wikipedia
Upwind scheme — In computational fluid dynamics, the upwind schemes are any of a class of discretization methods to solve hyperbolic partial differential equations numerically. The wave equation, the advection equation, the Euler equations in fluid dynamics, etc … Wikipedia
Numerical ordinary differential equations — Illustration of numerical integration for the differential equation y = y,y(0) = 1. Blue: the Euler method, green: the midpoint method, red: the exact solution, y = et. The step size is h = 1.0 … Wikipedia
Explicit and implicit methods — In applied mathematics, explicit and implicit methods are approaches used in computer simulations of physical processes, or in other words, they are numerical methods for solving time variable ordinary and partial differential equations.Explicit… … Wikipedia
Godunov's theorem — Godunov s theorem, also known as Godunov s order barrier theorem, is a mathematical theorem important in the development of the theory of high resolution schemes for the numerical solution of partial differential equations.Professor Sergei K.… … Wikipedia