计算水力学第三章.ppt

计算水力学第三章.ppt

《计算水力学第三章.ppt》由会员分享，可在线阅读，更多相关《计算水力学第三章.ppt（102页珍藏版）》请在冰豆网上搜索。

计算水力学计算水力学第三章第三章有限差分的基本理论有限差分的基本理论第一节第一节基本概念基本概念n一维对流方程一维对流方程n计算平面计算平面为的上半平面。

在平面上画出为的上半平面。

在平面上画出两族平行于坐标轴的直线两族平行于坐标轴的直线,把求解域分成矩形把求解域分成矩形的的计算计算网格网格。

网格线的交点称为。

网格线的交点称为节点节点,方向方向上网格线之间的距离上网格线之间的距离称为称为空间步长空间步长,轴方轴方向上网格线之间的距离向上网格线之间的距离t称为称为时间步长时间步长平面、计算网格平面、计算网格网格节点网格节点节点函数值节点函数值n网格剖分使得每一空间步长、时间步长网格剖分使得每一空间步长、时间步长均相等，则称该网格为一均相等，则称该网格为一均匀网格均匀网格，否，否则称之为则称之为非均匀网格非均匀网格n数值解主要是求解节点上的末知变量的数值解主要是求解节点上的末知变量的数值，利用数值，利用有限的节点有限的节点上的值来代替整上的值来代替整个求解域内的个求解域内的连续函数连续函数值。

值。

n概念：

离散、插值、误差概念：

离散、插值、误差构造差分方程、分析数值误差构造差分方程、分析数值误差第二节第二节偏导数的差商近似偏导数的差商近似一、差分、差商的基本概念一、差分、差商的基本概念n解析函数n导数定义n差分n差商n向前差分n向后差分n中心差分n一阶导数，对应的差分称为一阶差分。

一阶导数，对应的差分称为一阶差分。

对一阶差分再作一阶差分，所得到的称对一阶差分再作一阶差分，所得到的称之为二阶差分。

二阶向前差分：

之为二阶差分。

二阶向前差分：

任何阶差分都可以由其低一阶的差分得到：

任何阶差分都可以由其低一阶的差分得到：

函数的差分与自变量的差分之比，即为函函数的差分与自变量的差分之比，即为函数对自变量的差商数对自变量的差商n一阶向前差商一阶向前差商n一阶向后差商一阶向后差商n一阶中心差商一阶中心差商n二阶中心差商二阶中心差商二、偏导数的差商近似二、偏导数的差商近似展开法展开法n通过对差商近似点（，）的通过对差商近似点（，）的展开，可以分析差商对偏导近似展开，可以分析差商对偏导近似的精度的精度n一阶向前差商一阶向前差商n一阶向后差商一阶向后差商n二阶中心差商二阶中心差商边界处偏导数的差商近似n对点（，）进行展开n构造一阶偏导数的二阶精度的差商近似必须有解得n构造二阶偏导数的差商近似必须有解得n构造二阶偏导数具有二阶精度的差商近似必须有n解得：

多项式插值法多项式插值法用多项式插值法把待求函数表示成含待定系数的解析函数，由节点函数值确定该系数，然后对此函数求偏导数，得到逼近偏导数的差商表达式。

设函数可用抛物插值公式来近似：

设原点在点的位置上，则有解出待定系数n用高阶多项式插值可得到高阶差商表达式。

用高阶多项式插值可得到高阶差商表达式。

n高阶多项式插值具有龙格不稳定性，使得插值高阶多项式插值具有龙格不稳定性，使得插值对计算误差十分敏感。

对计算误差十分敏感。

n多项式插值法在计算流体力学中多用于处理边多项式插值法在计算流体力学中多用于处理边界处的差商近似。

界处的差商近似。

n偏导数的差商近似还有其它多种方法，但最终偏导数的差商近似还有其它多种方法，但最终均需用展开来计算其近似的误差，均需用展开来计算其近似的误差，因此在实际计算中通常均用展开因此在实际计算中通常均用展开法来构造，因为此法在构造差商近似的同时还法来构造，因为此法在构造差商近似的同时还得出了其近似的误差精度。

得出了其近似的误差精度。

第三节第三节差分方程差分方程n偏导数用其差商近似来代替偏导数用其差商近似来代替n偏微分方程转变为相应的代数方程称之偏微分方程转变为相应的代数方程称之为差分方程。

为差分方程。

对流方程对流方程在点（，）成立在点（，）成立在点（，）的对流方程可以近似数差分方程n设，求解域（），n定解条件n离散表达定解问题定解问题格式差分方程格式差分方程n例定解问题n采用格式取，0.5，C=1.0则b取，2，C=1.0则n对同一定解问题的同一差分格式（）其不同的空间与时间步长，将得到不同的结果，如果作为原始定解问题的近似解，那一个解精度高呢？

。

n不稳定的解是不能作为原定解问题的近似解的。

n偏导数的差商近似并非一种，同一偏微分方程的差分方程也并非一个，可以有若干个，对原始定解问题也相应有若干种差分格式。

格式格式格式格式蛙跳格式蛙跳格式n显式格式：

由第时间层上的值，可直接算出显式格式：

由第时间层上的值，可直接算出第时间层上的值的格式。

第时间层上的值的格式。

n隐式格式：

不能直接从时间层上值直接解出，隐式格式：

不能直接从时间层上值直接解出，需联立求解层上的值的格式。

需联立求解层上的值的格式。

对同一个定解问题，可以有多种差分格式，对同一个定解问题，可以有多种差分格式，多种步长参数来近似，从而也得到若干个差分多种步长参数来近似，从而也得到若干个差分近似解。

那么这些解是否可以都作为原定解问近似解。

那么这些解是否可以都作为原定解问题的近似解？

那些解精度高？

为什么？

题的近似解？

那些解精度高？

为什么？

相容性、稳定性及收敛性分析相容性、稳定性及收敛性分析第四节第四节截断误差和相容性截断误差和相容性n以格式为例n等价方程等价方程n截断误差截断误差n格式的截断误差n格式的截断误差n蛙跳格式的截断误差n对流方程的差分方程等价形式n差分方程和相应的微分方程相容n定解条件n差分算子n截断误差n定解问题相容第五节第五节收敛性收敛性n相容性：

是指当自变量的步长趋于零时相容性：

是指当自变量的步长趋于零时,差分格式与微分问题的截误差的范数是差分格式与微分问题的截误差的范数是否趋于零否趋于零,从而可看出是否能用此差分格从而可看出是否能用此差分格式来逼近微分问题。

式来逼近微分问题。

n收敛性：

是指当自变量步长趋于零时收敛性：

是指当自变量步长趋于零时,要要求差分格式的解趋于微分方程定解问题求差分格式的解趋于微分方程定解问题的解。

要求差分格式的解的解。

要求差分格式的解（数值解数值解）与微与微分方程定解问题的解分方程定解问题的解（精确解精确解）是一致的。

是一致的。

n差分格式的解差分格式的解n微分问题的解微分问题的解n离散误差离散误差n差分格式收敛差分格式收敛相容性是收敛性的必要条件相容性是收敛性的必要条件,相容性相容性是形式上的逼近，收敛性是解的逼近，是形式上的逼近，收敛性是解的逼近，相容性不一定能保证收敛性相容性不一定能保证收敛性n例微分方程定解问题n解析解为将，等分为段则步长n差分解为微分问题格式n离散误差n由截断误差分析有n当和,即,则n说明当时,本问题的格式收敛。

这种离散化误差的最大绝对值趋于零的收敛性情况称为一致性收敛。

第六节第六节稳定性稳定性n定解问题依赖区间AB和决定域pAB影响区域n同一微分问题同一微分问题,当采用不同差分格式时当采用不同差分格式时,其依赖区间、决定区域和影响区域可以其依赖区间、决定区域和影响区域可以是不一致的。

是不一致的。

n依赖区间、决定区域和影响区域是由差依赖区间、决定区域和影响区域是由差分格式本身的构造所决定的并与步长比分格式本身的构造所决定的并与步长比有关有关定解问题定解问题FTBSFTBS格式计算格式计算n假设在第假设在第j层上的第层上的第i点点,由于计算误差，由于计算误差，得到。

得到。

设设i,，相应于，相应于FTBS格式格式n算例表明了当值不同时计算误差所产生的影响。

n误差逐渐衰减传播n误差无衰减传播n误差震荡放大传播n数值误差有不同的传播方式,格式使误差逐渐衰减传播称为差分格式稳定，否则称为不稳定。

单增长型的不稳定称为静力不稳定性单增长型的不稳定称为静力不稳定性过冲型振荡的不稳定称为动力不稳定过冲型振荡的不稳定称为动力不稳定vonNeumann稳定性分析方法稳定性分析方法n定解问题nFTBS格式n初值误差n误差传播方程n误差展开成傅氏级数n代入误差传播方程n对任意的有为放大因子nFTBS格式稳定条件格式稳定条件nFTBS格式稳定条件n格式为一不稳定格式n格式稳定条件n蛙跳格式稳定条件nVonNeumann稳定性分析法主要用于稳定性分析法主要用于线性初值问题的稳定性分析。

对于非线线性初值问题的稳定性分析。

对于非线性问题用局部线性化的方法加以推广。

性问题用局部线性化的方法加以推广。

局部线性化方法假定非线性系数变化得局部线性化方法假定非线性系数变化得很缓慢，因而可用局部网格结点上的函很缓慢，因而可用局部网格结点上的函数值代入后作为常数处理，并认为每一数值代入后作为常数处理，并认为每一网格结点上的计算稳定性与相邻结点无网格结点上的计算稳定性与相邻结点无关，以网格结点上最小的局部稳定极限关，以网格结点上最小的局部稳定极限值作为整个差分问题的稳定极限值。

值作为整个差分问题的稳定极限值。

第七节第七节Lax等价定理等价定理相容性是收敛性的必要条件，稳定性与收敛性相容性是收敛性的必要条件，稳定性与收敛性有一定的联系。

有一定的联系。

Lax等价定理就是阐述相容性、等价定理就是阐述相容性、收敛性和稳定性三者之间的关系的。

收敛性和稳定性三者之间的关系的。

nLax等价定理等价定理:

对一个适定的线性微分问题及一对一个适定的线性微分问题及一个与其相容的差分格式，如果该格式稳定则必个与其相容的差分格式，如果该格式稳定则必收敛，不稳定必不收敛。

换言之，若线性微分收敛，不稳定必不收敛。

换言之，若线性微分问题适定，差分格式相容，则稳定性是收敛性问题适定，差分格式相容，则稳定性是收敛性的必要和充分的条件。

的必要和充分的条件。

n根据此定理,在线性适定和格式相容的条件下，只要证明了格式是稳定的，则一定收敛；若不稳定，则不收敛。

由于收敛性的证明往往比稳定性更难，故人们就可以把注意力集中在稳定性的研究上。

第八节第八节差分方程数值效应差分方程数值效应n微分方程是描述物理量在时间和空间上的连续变化的规律n差分方程来描述离散化后物理量的变化规律n离散误差使原系统的物理性质和规律遭到歪曲和破坏的作用称为数值效应或离散近似的伪物理效应。

n必须对这些效应有明确的概念，从物理上来考虑数值格式的合理性，减少数值效应的影响。

一、一、“逆风逆风”效应效应n物质的对流输运出现了与波速相反方向传播的不合理现象，称为“逆风”效应，是一伪物理现象的数值效应。

对流方程FTCS格式n假定在某瞬时在某一断面处引入某一物理量，由表可见，物理量向上、向下游两个方向传播，出现由表可见，物理量向上、向下游两个方向传播，出现了与波速相反方向传播的不合理现象，称之为了与波速相反方向传播的不合理现象，称之为“逆风逆风”效应。

效应。

FTCS格式所描述的物理量的运动规律与它所近似的原格式所描述的物理量的运动规律与它所近似的原问题固有的规律相差甚大，不仅计算结果误差很大，问题固有的规律相差甚大，不仅计算结果误差很大，而且也往往是引起差分格式不稳定的一个因素，前面而且也往往是引起差分格式不稳定的一个因素，前面已证明了该差分格式是无条件不稳定的。

已证明了该差分格式是无条件不稳定的。

对流方程采用一阶精度对流方程采用一阶精度FTBS格式或格式或FTFS格式能近似格式能近似原问题的物理现象。

采用何种格式还与波速的方向有原问题的物理现象。

采用何种格式还与波速的方向有关，当为正时采用格式，当为负时采用关，当为正时采用格式，当为负时采用格式。

若的符号在计算过程中会改变，例格式。

若的符号在计算过程中会改变，例如潮水河道，则可以采用如潮水河道，则可以采用“逆风逆风”格式：

格式：

“逆风逆风”格式格式二、物理耗散与弥散二、物理耗散与弥散n一维行波的波高n在一固定时刻，空间上相差一个波长其波高相等n所以称为波数n在一固定位置，时间相差一个周期，其波高相等n为频率。

n表示单位时间传播的距离，称为相速度。

n物理耗散物理耗散是指波幅因阻尼作用而衰减的现象n弥散弥散是指波的相速度随波数发生变化的现象。

n深水波相速度为弥散波