3第二章-有限差分方法基础解读课件

计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput第二章有限差分方法基础§2.1有限差分方法概述§2.2导数的数值逼近方法§2.3差分格式的性质§2.4发展方程的稳定性分析第二章有限差分方法基础§2.1有限差分方法概述§2.1有限差分方法概述以一维非定常热传导方程为例，介绍有限差分方法的概念、简单构造方法和求解过程。2.1.1基本方程和定解问题方程(2.1.1)和初边条件(2.1.2)构成了一个适定的定解问题。有限差分方法：对于一个偏微分方程，如果把方程中的所有偏导数近似地用代数差商(AlgebraicDifferenceQuotient)代替，则可以用一组代数方程近似地替代这个偏微分方程，进而得到数值解，这种方法称为有限差分方法(FiniteDifferenceMethod)。§2.1有限差分方法概述以一维非定常热传导方程2.1.2求解域及偏导数的离散化为了用有限差分方法求解式(2.1.1)，需要把其中的偏导数表示为代数形式，为此，首先要把自变量从连续的分布变为离散形式。这个过程称为求解域的离散化。

1.空间求解域的离散化把空间求解域分为M段（均匀剖分）

2.时间变量的离散化把感兴趣的时间段(t=T之前)分为N段（均匀剖分），则时间方向的求解域可以划分为2.1.2求解域及偏导数的离散化为了用有限差分

求解域被划分为一系列离散的时空网格点图2.1求解域的离散化

3.解的离散表示目标：求出所有网格点上物理量u的近似解。求解域被划分为一系列离散的时空网格点图2.1

4.导数的数值逼近把方程中的偏导数项近似表示为代数形式。4.导数的数值逼近把方程中的偏导数项近似表示为代3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.1.3差分格式同一偏导数可以有不同的近似方法，不同的导数近似方法导致方程的不同的有限差分近似。FTCS(ForwarddifferenceinTime,CentraldifferenceinSpace)格式时间方向用前差近似，空间二阶导数用中心差分近似。对初始条件和边界条件的离散化式(2.1.9)~(2.1.12)称为方程(2.1.1)的一个有限差分方程或有限差分格式(finitedifferencescheme)。2.1.3差分格式同一偏导数可以有不同的近似方法，不同的2.BTCS(BackwarddifferenceinTime,CentraldifferenceinSpace)格式时间方向用后差近似，空间二阶导数用中心差分近似。在研究数值方法时，通常把tn时刻的物理量视为已知量，而把tn+1时刻的物理量作为待求的未知量。因此，式(2.1.13)可以改写成2.BTCS(Backwarddifference2.1.4差分方程的求解FTCS格式可以改写为可见，在FTCS格式中，某一点的数值解只依赖于前一时间步的三个点，如图2.2所示。图2.2：FTCS格式的模板点2.1.4差分方程的求解FTCS格式可以改写为可见，在FTCS格式的求解过程FTCS格式的求解过程2.BTCS格式可以改写为跟FTCS格式不同，BTCS格式中同时涉及到n+1时刻的多个未知量，不能递推求解，称为隐式格式(implicitscheme)。图2.3：BTCS格式的模板点BTCS格式的求解过程2.BTCS格式可以改写为跟FTCS格式不同，BTC3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.1.5用时间相关方法求解定常问题考虑非定常热传导方程和定解条件2.1.5用时间相关方法求解定常问题考虑非定常热传导方程3第二章-有限差分方法基础解读课件BTCS格式的求解过程FTCS格式的求解过程BTCS格式的求解过程FTCS格式的求解过程§2.2导数的数值逼近方法2.2.1精度分析在上一节，我们得到了一阶偏导数的前差、后差和中心差分近似，以及二阶导数的中心差分近似。这些近似方法逼近偏导数的程度如何呢？可以用Taylor展开式进行分析。§2.2导数的数值逼近方法2.2.1精度分析3第二章-有限差分方法基础解读课件一般来讲，对偏导数的近似精度越高，差分格式的精度越高。一般来讲，对偏导数的近似精度越高，差分格式的精度例：一维非定常热传导方程的FTCS格式中涉及的导数差分近似的精度。例：一维非定常热传导方程的FTCS格式中涉及的导数差分近似的2.2.2导数差分近似的待定系数法2.2.2导数差分近似的待定系数法3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.2.3导数差分近似方法的差分算子法1.差分算子的定义算子，一种前置运算符。算子和它后面的作用量一起代表一种确定的运算过程。在引入差分算子的定义之前，先介绍一种特殊的算子——移位算子。移位算子的运算规则为移位算子的下标表示移位的方向，上标表示移位的步数。2.2.3导数差分近似方法的差分算子法1.差分算子的定差分算子：移位算子和可以表示为移位算子函数的算子。差分方法中常用的算子：差分算子：移位算子和可以表示为移位算子函数的算子。差分方法中2.差分算子之间的关系2.差分算子之间的关系所有的差分算子均可用Taylor展开式来估算截断误差项的量级。所有的差分算子均可用Taylor展开式来估算截断误差项的量级3.微分算子与差分算子的关系3.微分算子与差分算子的关系4.导数的近似根据差分算子之间的转化关系，可以建立微分算子与其它差分算子之间的联系，从而得到导数的数值近似公式。即：4.导数的近似根据差分算子之间的转化关系，可以即：与待定系数法得到的结果一致。即：与待定系数法得到的结果一致。即：即：5.紧致格式从上面的推导可以看出，导数的有限差分近似精度越高，所需要的模板点越多。对于一阶导数，一般需要5个点才能得到四阶精度的差分近似。模板点数太多不仅使数值方法变得复杂，也给边界附近的处理带来一定困难。紧致格式：用较少的模板点构造导数的高阶近似。5.紧致格式从上面的推导可以看出，导数的有限差3第二章-有限差分方法基础解读课件基于Pade近似的导数近似方法，称为紧致格式(compactscheme)。基于Pade近似的导数近似方法，称为紧致格式(compac3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件§2.3差分格式的性质2.3.1范数的定义及性质1.向量范数§2.3差分格式的性质2.3.1范数的定义及性质1.2.算子范数2.算子范数3第二章-有限差分方法基础解读课件2.3.2差分格式的精度差分格式是微分方程的近似，通常用局部截断误差(localtruncationerror)衡量差分格式逼近微分方程的程度。2.3.2差分格式的精度差分格式是微分方程的近似，通常用3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件如果时间步长和空间步长之间满足一定的关系，FTCS格式时间方向可达到二阶精度，空间方向可达到四阶精度。根据差分格式精度的定义，按照上面的分析，FTCS格式时间方向是一阶精度，空间方向是二阶精度。如果时间步长和空间步长之间满足一定的关系，FTCS格式时间方2.3.3差分格式的相容性截断误差是在网格点上逐点定义的。定义中每个网格点上的数值解构成一个解向量，每一个网格点上的截断误差也构成一个向量。因此，可以用向量范数来刻画差分格式的局部截断误差。2.3.3差分格式的相容性截断误差是在网格点上逐点定义的2.3.4差分格式的收敛性和稳定性1.差分方程的矩阵形式考虑线性的发展方程（双曲型方程和抛物型方程）的差分格式。发展型方程的一般形式：以非定常热传导方程的FTCS格式为例，将差分格式写成矩阵形式：FTCS格式：解向量记为：考虑到边界条件，则差分格式可以写为：2.3.4差分格式的收敛性和稳定性1.差分方程的矩阵形2.整体截断误差局部截断误差：差分方程逼近微分方程的程度整体截断误差：差分方程的解逼近微分方程的精确解的程度2.整体截断误差局部截断误差：差分方程逼近微分方程的程度整3第二章-有限差分方法基础解读课件3.差分格式的收敛性和稳定性差分格式的收敛性对于保证数值解的有效性是非常重要的。如果差分格式是收敛的，那么，当计算网格足够密时，数值解将相当接近精确解。差分格式的稳定性等价于差分方程数值解的一致有界性。3.差分格式的收敛性和稳定性差分格式的收敛性对于保证数值解3第二章-有限差分方法基础解读课件上述定理建立了算子范数的一致有界性与稳定性之间的关系。当差分格式稳定时，整体截断误差和局部截断误差量级相同。上述定理建立了算子范数的一致有界性与稳定性之间的关系。当差分3第二章-有限差分方法基础解读课件Lax等价性定理是计算流体力学中的一个重要定理。直接分析差分格式的收敛性比较困难，而稳定性分析则比较简单。Lax定理告诉我们，在一定条件下，收敛性和稳定性是等价的；通过稳定性分析，即可确定差分格式的收敛条件。4.稳定性的意义Lax等价性定理是计算流体力学中的一个重要定理。直接分析差分§2.4发展方程的稳定性分析2.4.1矩阵方法§2.4发展方程的稳定性分析2.4.1矩阵方法3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.4.2VonNeumann稳定性理论2.4.2VonNeumann稳定性理论3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.4.3稳定性分析实例2.4.3稳定性分析实例3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput第二章有限差分方法基础§2.1有限差分方法概述§2.2导数的数值逼近方法§2.3差分格式的性质§2.4发展方程的稳定性分析第二章有限差分方法基础§2.1有限差分方法概述§2.1有限差分方法概述以一维非定常热传导方程为例，介绍有限差分方法的概念、简单构造方法和求解过程。2.1.1基本方程和定解问题方程(2.1.1)和初边条件(2.1.2)构成了一个适定的定解问题。有限差分方法：对于一个偏微分方程，如果把方程中的所有偏导数近似地用代数差商(AlgebraicDifferenceQuotient)代替，则可以用一组代数方程近似地替代这个偏微分方程，进而得到数值解，这种方法称为有限差分方法(FiniteDifferenceMethod)。§2.1有限差分方法概述以一维非定常热传导方程2.1.2求解域及偏导数的离散化为了用有限差分方法求解式(2.1.1)，需要把其中的偏导数表示为代数形式，为此，首先要把自变量从连续的分布变为离散形式。这个过程称为求解域的离散化。

1.空间求解域的离散化把空间求解域分为M段（均匀剖分）

2.时间变量的离散化把感兴趣的时间段(t=T之前)分为N段（均匀剖分），则时间方向的求解域可以划分为2.1.2求解域及偏导数的离散化为了用有限差分

求解域被划分为一系列离散的时空网格点图2.1求解域的离散化

3.解的离散表示目标：求出所有网格点上物理量u的近似解。求解域被划分为一系列离散的时空网格点图2.1

4.导数的数值逼近把方程中的偏导数项近似表示为代数形式。4.导数的数值逼近把方程中的偏导数项近似表示为代3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.1.3差分格式同一偏导数可以有不同的近似方法，不同的导数近似方法导致方程的不同的有限差分近似。FTCS(ForwarddifferenceinTime,CentraldifferenceinSpace)格式时间方向用前差近似，空间二阶导数用中心差分近似。对初始条件和边界条件的离散化式(2.1.9)~(2.1.12)称为方程(2.1.1)的一个有限差分方程或有限差分格式(finitedifferencescheme)。2.1.3差分格式同一偏导数可以有不同的近似方法，不同的2.BTCS(BackwarddifferenceinTime,CentraldifferenceinSpace)格式时间方向用后差近似，空间二阶导数用中心差分近似。在研究数值方法时，通常把tn时刻的物理量视为已知量，而把tn+1时刻的物理量作为待求的未知量。因此，式(2.1.13)可以改写成2.BTCS(Backwarddifference2.1.4差分方程的求解FTCS格式可以改写为可见，在FTCS格式中，某一点的数值解只依赖于前一时间步的三个点，如图2.2所示。图2.2：FTCS格式的模板点2.1.4差分方程的求解FTCS格式可以改写为可见，在FTCS格式的求解过程FTCS格式的求解过程2.BTCS格式可以改写为跟FTCS格式不同，BTCS格式中同时涉及到n+1时刻的多个未知量，不能递推求解，称为隐式格式(implicitscheme)。图2.3：BTCS格式的模板点BTCS格式的求解过程2.BTCS格式可以改写为跟FTCS格式不同，BTC3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.1.5用时间相关方法求解定常问题考虑非定常热传导方程和定解条件2.1.5用时间相关方法求解定常问题考虑非定常热传导方程3第二章-有限差分方法基础解读课件BTCS格式的求解过程FTCS格式的求解过程BTCS格式的求解过程FTCS格式的求解过程§2.2导数的数值逼近方法2.2.1精度分析在上一节，我们得到了一阶偏导数的前差、后差和中心差分近似，以及二阶导数的中心差分近似。这些近似方法逼近偏导数的程度如何呢？可以用Taylor展开式进行分析。§2.2导数的数值逼近方法2.2.1精度分析3第二章-有限差分方法基础解读课件一般来讲，对偏导数的近似精度越高，差分格式的精度越高。一般来讲，对偏导数的近似精度越高，差分格式的精度例：一维非定常热传导方程的FTCS格式中涉及的导数差分近似的精度。例：一维非定常热传导方程的FTCS格式中涉及的导数差分近似的2.2.2导数差分近似的待定系数法2.2.2导数差分近似的待定系数法3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件2.2.3导数差分近似方法的差分算子法1.差分算子的定义算子，一种前置运算符。算子和它后面的作用量一起代表一种确定的运算过程。在引入差分算子的定义之前，先介绍一种特殊的算子——移位算子。移位算子的运算规则为移位算子的下标表示移位的方向，上标表示移位的步数。2.2.3导数差分近似方法的差分算子法1.差分算子的定差分算子：移位算子和可以表示为移位算子函数的算子。差分方法中常用的算子：差分算子：移位算子和可以表示为移位算子函数的算子。差分方法中2.差分算子之间的关系2.差分算子之间的关系所有的差分算子均可用Taylor展开式来估算截断误差项的量级。所有的差分算子均可用Taylor展开式来估算截断误差项的量级3.微分算子与差分算子的关系3.微分算子与差分算子的关系4.导数的近似根据差分算子之间的转化关系，可以建立微分算子与其它差分算子之间的联系，从而得到导数的数值近似公式。即：4.导数的近似根据差分算子之间的转化关系，可以即：与待定系数法得到的结果一致。即：与待定系数法得到的结果一致。即：即：5.紧致格式从上面的推导可以看出，导数的有限差分近似精度越高，所需要的模板点越多。对于一阶导数，一般需要5个点才能得到四阶精度的差分近似。模板点数太多不仅使数值方法变得复杂，也给边界附近的处理带来一定困难。紧致格式：用较少的模板点构造导数的高阶近似。5.紧致格式从上面的推导可以看出，导数的有限差3第二章-有限差分方法基础解读课件基于Pade近似的导数近似方法，称为紧致格式(compactscheme)。基于Pade近似的导数近似方法，称为紧致格式(compac3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件§2.3差分格式的性质2.3.1范数的定义及性质1.向量范数§2.3差分格式的性质2.3.1范数的定义及性质1.2.算子范数2.算子范数3第二章-有限差分方法基础解读课件2.3.2差分格式的精度差分格式是微分方程的近似，通常用局部截断误差(localtruncationerror)衡量差分格式逼近微分方程的程度。2.3.2差分格式的精度差分格式是微分方程的近似，通常用3第二章-有限差分方法基础解读课件3第二章-有限差分方法基础解读课件如果时间步长和空间步长之间满足一定的关系，FTCS格式时间方向可达到二阶精度，空间方向可达到四阶精度。根据差分格式精度的定义，按照上面的分析，FTCS格式时间方向是一阶精度，空间方向是二阶精度。如果时间步长和空间步长之间满足一定的关系，FTCS格式时间方2.3.3差分格式的相容性截断误差是在网格点上逐点定义的。定义中每个网格点上的数值解构成一个解向量，每一个网格点上的截断误差也构成一个向量。因此，可以用向量范数来刻画差分格式的局部截断误差。2.3.3差分格式的相容性截断误差是在网格点上逐点定义的2.3.4差分格式的收敛性和稳定性1.差分方程的矩阵形式考虑线性的发展方程（双曲型方程和抛物型方程）的差分格式。发展型方程的一般形式：以非定常热传导方程的FTCS格式为例，将差分格式写成矩阵形式：FTCS格式：解向量记为：考虑到边界条件，则差分格式可以写为：2.3.4差分格式的收敛性和稳定性1.差分方程的矩阵形2.整体截断误差局部截断误差：差分方程逼近微分方程的程度整体截断误差：差分方程的解逼近微分方程的精确解的程度2.整体截断误差局部截