有限差分法.docx

《有限差分法.docx》由会员分享，可在线阅读，更多相关《有限差分法.docx（10页珍藏版）》请在三一办公上搜索。

1、有限差分法有限差分法 有限差分法 finite difference method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似， 积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组 ， 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原微分

2、方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n1层的近似值时要用到第n层的近似值 ，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差

3、的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以T

4、aylor级 数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 从差分的空间形式来考虑，可分为中心格式和逆风格式。 考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

5、 构造差分的方法有多种形式， 目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等， 其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几 种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式 ，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数 形

6、式，便构成不同的有限元方法 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟 。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单 元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的 基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元 上的近似解构成。 在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同 ，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二 乘法和伽辽金法，从计算单元

7、网格的形状来划分，有三角形网格、四边形网格和多边形 网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合 同样构成不同的有限元计算格式。 对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数 ； 最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小； 在配置法中，先在计算域 内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘 积表示，但最常用的多项式插值函数。 有限元插值函数分为两大类， 一类只要求插值多项式本身在

8、插值点取已知值，称为拉格朗日(Lagrange)多项式插值； 一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。 单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。 在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。 对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线 性插值函数、二阶或更高阶插值函数等。 对于

9、有限元方法，其基本思路和解题步骤可归纳为 (1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。 (2)区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干 相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工 作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节 点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。 (3)确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条 件的插值函数作为单元基函数。有限元方法中

10、的基函数是在单元中选取的，由于各单元 具有规则的几何形状，在选取基函数时可遵循一定的法则。 (4)单元分析：将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近；再将近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数(即单元中各节点 的参数值)的代数方程组，称为单元有限元方程。 (5)总体合成：在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。 (6)边界条件的处理：一般边界条件有三种形式，分为本质边界条件(狄里克雷边界条件 )、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。对于自然边界条件， 一般在积分表达式中可自动得

11、到满足。对于本质边界条件和混合边界条件，需按一定法则对总体有限元方程进行修正满足。 (7)解有限元方程：根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。 有限体积法又称为控制体积法。 其基本思路是：将计算区域 划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。其中的未知数是网格点上的因变量 的数值。为了求出控制体积的积分，必须假定 值在网格点之间的变化规律，即假设 值 的分段的分布的分布剖面。从积分区域的选取方法看来，有限体积法属于加权剩余法中 的子区域

12、法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。简 言之，子区域法属于有限体积发的基本方法。 有限体积法的基本思路易于理解，并能得出直接的物理解释。离散方程的物理意义，就 是因变量 在有限大小的控制体积中的守恒原理，如同微分方程表示因变量在无限小的控 制体积中的守恒原理一样。 限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足， 对整个计算区域，自然也得到满足。这是有限体积法吸引人的优点。 有一些离散方法，例如有限差分法，仅当网格极其细密时，离散方程才满足积分守恒；而有限体积法即使在粗网格情况下，也显示出准确的积分守恒。 就离散方法而言，有限体积法可视作

13、有限单元法和有限差分法的中间物。有限单元法必 须假定值在网格点之间的变化规律，并将其作为近似解。有限差分法只 考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求的结点值，这与有限差分法相类似；但有限体积法在寻求控制体积的积分时，必须假定值在网格 点之间的分布，这又与有限单元法相类似。在有限体积法中，插值函数只用于计算控制 体积的积分，得出离散方程之后，便可忘掉插值函数；如果需要的话，可以对微分方程 中不同的项采取不同的插值函数。 第二节 有限差分法(FDM) 一、基本思想 按时间步长和空间步长将时间和空间区域剖分成若干网格，用未知函数在网格结(节)点上的值所构成的差分近似代替所

14、用偏微分方程中出现的各阶导数，从而把表示变量连续变化关系的偏微分方程离散为有限个代数方程，然后解此线性代数方程组，以求出溶质在各网格结(节)点上不同时刻的浓度。 二、基本步骤 剖分渗流区，确定离散点。将所研究的水动力弥散区域按某种几何形状(如矩形、任意多边形等)剖分成网络系统。 建立水动力弥散问题的差分方程组。 求解差分方程组。采用各种迭代法，如点逐次超松驰方法(SOR)、线逐次超松驰方法(LSOR)、迭代的交替方向隐式方法(IADI)及强隐式方法(SID)等。 三、几种导数的差分近似 1、一阶向前差分 设f(x)为任一足够光滑的函数,将其沿x的正向展开为Tayler级数: 于是 若略去O(x

15、), 则有 2、一阶向后差分 3、一阶中心差分 4、二阶导数的差分 四、差分方程的相容性、收敛性、稳定性 1、相容性 导数与其差分近似式之间存在截断误差。因此，差分方程的解并不是严格的，而是近似地满足原来的偏微分方程。但是，当时间步长t和空间步长x都趋近于零时，差分方程的截差(截断误差)也趋近于零，差分方程的极限形式就是原偏微分方程。这时，认为差分方程与偏微分方程是相容的，这种相容性表示差分方程“收敛”于原偏微分方程。 2、收敛性 指差分方程的解，即当步长t、x0时收敛于原偏微分方程的解。 3、稳定性 由于差分方程的求解是以步进方式进行的，在逐步推进的过程中，误差也逐步积累。若这种误差积累保持

16、有界，则差分方程是稳定的，若这种误差积累无界则差分方程是不稳定的。 五、差分格式 1、显式差分格式 将及差分中的浓度取为tn，便可得到其显式差分格式：2、隐式差分格式 将及差分中的浓度取为tn+1，便可得到其隐式差分格式： 3、Crank-Nicolson差分格式 取显式与隐式差分格式的平均，便得到该格式： 上述三种差分格式截断误差不同，稳定性条件也不同。 有限差分法 finite difference method 用差分代替微分，是有限差分法的基本出发点。是一种微分方程和积分微分方程数值解的方法。 目录 隐藏 1 基本思想 2 主要内容 3 构造方法 有限差分法-基本思想 把连续的定解区域

17、用有限个离散点构成的网格来代替， 这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似， 积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组 ， 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 有限差分法-主要内容 如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差

18、分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n1层的近似值时要用到第n层的近似值 ，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。 有限差分法-构造方法 最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式