【转】Courant–Friedrichs–Lewy condition(CFL条件,柯朗-弗里德里希斯-列维条件)
2012-11-26 17:30 H_OU 阅读(6939) 评论(0) 编辑 收藏 举报In mathematics, the Courant–Friedrichs–Lewy condition (CFL condition) isa necessary condition for convergence while solving certain partial differential equations (usually hyperbolic PDEs) numerically by the method of finite differences.[1]
虽然在Citcom解地幔对流中,每一个时间步是有限元法,但是在时间演化上,仍然是有限差分。我现在,还没见到过有限元用在时间演化(时间导数)上的。
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 incorrect results. The condition is named after Richard Courant, Kurt Friedrichs, and Hans Lewy who described it in their 1928 paper.[2]
Heuristic description
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,[3] then this length 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. (当空间网格的长度减小时,时间步的长度也应该减小,恩,怪不得当我Citcom中,把空间步长减小时,时间演化就会减慢。原来是在这里啊)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[4] and the time step: the CFL condition relates the length of the time step to a function interval lengths of each spatial variable.
Operatively, theCFLcondition is commonly prescribed for those terms of the finite-difference approximation of general partial differential equations which model theadvection phenomenon.[5]
The one-dimensional case
For one-dimensional case, the CFL has the following form:
where
u is the velocity (whose dimension is Length/Time)
\Delta t is the time step (whose dimension is Time)
\Delta x is the length interval (whose dimension is Length),
C is a dimensionless constant which depends only on the particular equation to be solved.[6]
The dimensionless number
is called the Courant number.
例如,对于波动方程而言,niu我一般取0.1,当然,其实niu=0.5就行了,但是,演化不了多少步就会发散。
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 n-dimensional case is the following one:
Note that the interval length it is not required to be the same for each spatial variable \Delta x_i, 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 conditionis 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 \Delta 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 \Delta x requires a fourth order decrease in the time step \Delta 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,[7] 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.
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 现代计算机视觉入门之:什么是图片特征编码
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· Spring AI + Ollama 实现 deepseek-r1 的API服务和调用
· 数据库服务器 SQL Server 版本升级公告
· C#/.NET/.NET Core技术前沿周刊 | 第 23 期(2025年1.20-1.26)
· 程序员常用高效实用工具推荐,办公效率提升利器!