《计算水力学第三章教程文件.ppt》由会员分享，可在线阅读，更多相关《计算水力学第三章教程文件.ppt（102页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、计算水力学第三章平面、计算网格平面、计算网格 网格节点网格节点 节点函数值节点函数值 n网格剖分使得每一空间步长、时间步长网格剖分使得每一空间步长、时间步长均相等，则称该网格为一均相等，则称该网格为一均匀网格均匀网格，否，否则称之为则称之为非均匀网格非均匀网格n数值解主要是求解节点上的末知变量的数值解主要是求解节点上的末知变量的数值，利用数值，利用有限的节点有限的节点上的值来代替整上的值来代替整个求解域内的个求解域内的连续函数连续函数值。值。n概念：离散、插值、误差概念：离散、插值、误差 构造差分方程、分析数值误差构造差分方程、分析数值误差第二节第二节 偏导数的差商近似偏导数的差商近似一、差分
2、、差商的基本概念一、差分、差商的基本概念n解析函数 n导数定义 n差分 n差商 n向前差分 n向后差分 n中心差分 n一阶导数，对应的差分称为一阶差分。一阶导数，对应的差分称为一阶差分。对一阶差分再作一阶差分，所得到的称对一阶差分再作一阶差分，所得到的称之为二阶差分。二阶向前差分：之为二阶差分。二阶向前差分：任何阶差分都可以由其低一阶的差分得到：任何阶差分都可以由其低一阶的差分得到：函数的差分与自变量的差分之比，即为函函数的差分与自变量的差分之比，即为函数对自变量的差商数对自变量的差商 n一阶向前差商一阶向前差商 n一阶向后差商一阶向后差商 n一阶中心差商一阶中心差商 n二阶中心差商二阶中心差
3、商二、偏导数的差商近似二、偏导数的差商近似展开法展开法n通过对差商近似点（，）的通过对差商近似点（，）的展开，可以分析差商对偏导近似展开，可以分析差商对偏导近似的精度的精度 n一阶向前差商一阶向前差商n一阶向后差商一阶向后差商n二阶中心差商二阶中心差商边界处偏导数的差商近似n对点（，）进行展开 n构造一阶偏导数的二阶精度的差商近似必须有 解得 n构造二阶偏导数的差商近似必须有 解得 n构造二阶偏导数具有二阶精度的差商近似必须有n解得：多项式插值法多项式插值法用多项式插值法把待求函数表示成含待定系数的解析函数，由节点函数值确定该系数，然后对此函数求偏导数，得到逼近偏导数的差商表达式。设函数可用抛
4、物插值公式来近似：设原点在点的位置上，则有 解出待定系数 n用高阶多项式插值可得到高阶差商表达式。用高阶多项式插值可得到高阶差商表达式。n高阶多项式插值具有龙格不稳定性，使得插值高阶多项式插值具有龙格不稳定性，使得插值对计算误差十分敏感。对计算误差十分敏感。n多项式插值法在计算流体力学中多用于处理边多项式插值法在计算流体力学中多用于处理边界处的差商近似。界处的差商近似。n偏导数的差商近似还有其它多种方法，但最终偏导数的差商近似还有其它多种方法，但最终均需用展开来计算其近似的误差，均需用展开来计算其近似的误差，因此在实际计算中通常均用展开因此在实际计算中通常均用展开法来构造，因为此法在构造差商近
5、似的同时还法来构造，因为此法在构造差商近似的同时还得出了其近似的误差精度。得出了其近似的误差精度。第三节第三节 差分方程差分方程n偏导数用其差商近似来代替偏导数用其差商近似来代替n偏微分方程转变为相应的代数方程称之偏微分方程转变为相应的代数方程称之为差分方程。为差分方程。对流方程对流方程在点（，）成立在点（，）成立 在点（，）的对流方程可以近似 数 差分方程 n设，求解域（），n定解条件n离散表达定解问题定解问题 格式差分方程格式差分方程 n例定解问题 n采用格式 取，0.5，C=1.0 则 b取，2，C=1.0则 n对同一定解问题的同一差分格式（）其不同的空间与时间步长，将得到不同的结果，如
6、果作为原始定解问题的近似解，那一个解精度高呢？。n不稳定的解是不能作为原定解问题的近似解的。n偏导数的差商近似并非一种，同一偏微分方程的差分方程也并非一个，可以有若干个，对原始定解问题也相应有若干种差分格式。格式格式 格式格式 蛙跳格式蛙跳格式 n显式格式：由第时间层上的值，可直接算出显式格式：由第时间层上的值，可直接算出第时间层上的值的格式。第时间层上的值的格式。n隐式格式：不能直接从时间层上值直接解出，隐式格式：不能直接从时间层上值直接解出，需联立求解层上的值的格式。需联立求解层上的值的格式。对同一个定解问题，可以有多种差分格式，对同一个定解问题，可以有多种差分格式，多种步长参数来近似，从
7、而也得到若干个差分多种步长参数来近似，从而也得到若干个差分近似解。那么这些解是否可以都作为原定解问近似解。那么这些解是否可以都作为原定解问题的近似解？那些解精度高？为什么？题的近似解？那些解精度高？为什么？相容性、稳定性及收敛性分析相容性、稳定性及收敛性分析第四节第四节 截断误差和相容性截断误差和相容性n以格式为例 n等价方程等价方程 n截断误差截断误差 n格式的截断误差 n格式的截断误差 n蛙跳格式的截断误差 n对流方程的差分方程等价形式 n差分方程和相应的微分方程相容 n定解条件 n差分算子 n截断误差 n定解问题相容 第五节第五节 收敛性收敛性n相容性：是指当自变量的步长趋于零时相容性：
8、是指当自变量的步长趋于零时,差分格式与微分问题的截误差的范数是差分格式与微分问题的截误差的范数是否趋于零否趋于零,从而可看出是否能用此差分格从而可看出是否能用此差分格式来逼近微分问题。式来逼近微分问题。n收敛性：是指当自变量步长趋于零时收敛性：是指当自变量步长趋于零时,要要求差分格式的解趋于微分方程定解问题求差分格式的解趋于微分方程定解问题的解。要求差分格式的解的解。要求差分格式的解(数值解数值解)与微与微分方程定解问题的解分方程定解问题的解(精确解精确解)是一致的。是一致的。n差分格式的解差分格式的解 n微分问题的解微分问题的解 n离散误差离散误差 n差分格式收敛差分格式收敛相容性是收敛性的
9、必要条件相容性是收敛性的必要条件,相容性相容性是形式上的逼近，收敛性是解的逼近，是形式上的逼近，收敛性是解的逼近，相容性不一定能保证收敛性相容性不一定能保证收敛性 n例微分方程定解问题n解析解为 将，等分为段则步长n差分解为微分问题 格式 n离散误差 n由截断误差分析有 n当和,即,则 n说明当时,本问题的格式收敛。这种离散化误差的最大绝对值趋于零的收敛性情况称为一致性收敛。第六节第六节 稳定性稳定性n定解问题 依赖区间AB 和决定域pAB影响区域 n同一微分问题同一微分问题,当采用不同差分格式时当采用不同差分格式时,其依赖区间、决定区域和影响区域可以其依赖区间、决定区域和影响区域可以是不一致
10、的。是不一致的。n依赖区间、决定区域和影响区域是由差依赖区间、决定区域和影响区域是由差分格式本身的构造所决定的并与步长比分格式本身的构造所决定的并与步长比 有关有关 定解问题定解问题FTBSFTBS格式计算格式计算n假设在第假设在第j层上的第层上的第i点点,由于计算误差，由于计算误差，得到。得到。设设i,，相应于，相应于FTBS格式格式 n算例表明了当 值不同时计算误差所产生的影响。n 误差逐渐衰减传播n 误差无衰减传播n 误差震荡放大传播n数值误差有不同的传播方式,格式使误差逐渐衰减传播称为差分格式稳定，否则称为不稳定。单增长型的不稳定称为静力不稳定性单增长型的不稳定称为静力不稳定性过冲型振
11、荡的不稳定称为动力不稳定过冲型振荡的不稳定称为动力不稳定 von Neumann稳定性分析方法稳定性分析方法n定解问题nFTBS格式n初值误差 n误差传播方程 n误差展开成傅氏级数 n代入误差传播方程n对任意的有 为放大因子 nFTBS格式稳定条件格式稳定条件 nFTBS格式稳定条件n格式为一不稳定格式 n格式稳定条件n蛙跳格式稳定条件 nVon Neumann稳定性分析法主要用于稳定性分析法主要用于线性初值问题的稳定性分析。对于非线线性初值问题的稳定性分析。对于非线性问题用局部线性化的方法加以推广。性问题用局部线性化的方法加以推广。局部线性化方法假定非线性系数变化得局部线性化方法假定非线性系
12、数变化得很缓慢，因而可用局部网格结点上的函很缓慢，因而可用局部网格结点上的函数值代入后作为常数处理，并认为每一数值代入后作为常数处理，并认为每一网格结点上的计算稳定性与相邻结点无网格结点上的计算稳定性与相邻结点无关，以网格结点上最小的局部稳定极限关，以网格结点上最小的局部稳定极限值作为整个差分问题的稳定极限值。值作为整个差分问题的稳定极限值。第七节第七节 Lax等价定理等价定理 相容性是收敛性的必要条件，稳定性与收敛性相容性是收敛性的必要条件，稳定性与收敛性有一定的联系。有一定的联系。Lax等价定理就是阐述相容性、等价定理就是阐述相容性、收敛性和稳定性三者之间的关系的。收敛性和稳定性三者之间的
13、关系的。nLax等价定理等价定理:对一个适定的线性微分问题及一对一个适定的线性微分问题及一个与其相容的差分格式，如果该格式稳定则必个与其相容的差分格式，如果该格式稳定则必收敛，不稳定必不收敛。换言之，若线性微分收敛，不稳定必不收敛。换言之，若线性微分问题适定，差分格式相容，则稳定性是收敛性问题适定，差分格式相容，则稳定性是收敛性的必要和充分的条件。的必要和充分的条件。n根据此定理,在线性适定和格式相容的条件下，只要证明了格式是稳定的，则一定收敛；若不稳定，则不收敛。由于收敛性的证明往往比稳定性更难，故人们就可以把注意力集中在稳定性的研究上。第八节第八节 差分方程数值效应差分方程数值效应n微分方
14、程是描述物理量在时间和空间上的连续变化的规律n差分方程来描述离散化后物理量的变化规律n离散误差使原系统的物理性质和规律遭到歪曲和破坏的作用称为数值效应或离散近似的伪物理效应。n必须对这些效应有明确的概念，从物理上来考虑数值格式的合理性，减少数值效应的影响。一、一、“逆风逆风”效应效应n物质的对流输运出现了与波速相反方向传播的不合理现象，称为“逆风”效应，是一伪物理现象的数值效应。对流方程FTCS格式 n假定在某瞬时在某一断面处引入某一物理量，由表可见，物理量向上、向下游两个方向传播，出现由表可见，物理量向上、向下游两个方向传播，出现了与波速相反方向传播的不合理现象，称之为了与波速相反方向传播的
15、不合理现象，称之为“逆风逆风”效应。效应。FTCS格式所描述的物理量的运动规律与它所近似的原格式所描述的物理量的运动规律与它所近似的原问题固有的规律相差甚大，不仅计算结果误差很大，问题固有的规律相差甚大，不仅计算结果误差很大，而且也往往是引起差分格式不稳定的一个因素，前面而且也往往是引起差分格式不稳定的一个因素，前面已证明了该差分格式是无条件不稳定的。已证明了该差分格式是无条件不稳定的。对流方程采用一阶精度对流方程采用一阶精度FTBS格式或格式或FTFS格式能近似格式能近似原问题的物理现象。采用何种格式还与波速的方向有原问题的物理现象。采用何种格式还与波速的方向有关，当为正时采用格式，当为负时
16、采用关，当为正时采用格式，当为负时采用格式。若的符号在计算过程中会改变，例格式。若的符号在计算过程中会改变，例如潮水河道，则可以采用如潮水河道，则可以采用“逆风逆风”格式：格式：“逆风逆风”格式格式 二、物理耗散与弥散二、物理耗散与弥散n一维行波的波高 n在一固定时刻，空间上相差一个波长其波高相等 n所以称为波数 n在一固定位置，时间相差一个周期，其波高相等 n为频率。n表示单位时间传播的距离，称为相速度。n物理耗散物理耗散是指波幅因阻尼作用而衰减的现象n弥散弥散是指波的相速度随波数发生变化的现象。n深水波相速度 为弥散波或称色散波 n浅水波相速度为非弥散波、无耗散和色散的模型、无耗散和色散的
17、模型 n放大因子的模，说明波在传播过程中，经过时刻后振幅没有衰减。放大因子的幅角,表明传播时刻后，同一位置处波的相位迟后为，它是相函数的差；而相速度为与波数无关，没有色散现象。方程描述了无衰减、无弥散的物理现象。、有耗散的模型、有耗散的模型 分别为二阶、四阶耗散系数。n方程描述了有物理耗散而无色散的波运动。耗散系数引起波幅的衰减，但相速度不发生改变。、有弥散的模型、有弥散的模型 、称为三阶、五阶弥散系数 方程描述了有弥散的波动，波分量的振幅值不随时间变化，而相速度是波数的函数。三、数值耗散和弥散三、数值耗散和弥散n用差分方程逼近微分方程时引入了误差，用差分方程逼近微分方程时引入了误差，有时这些
18、误差项使计算结果的幅值衰减有时这些误差项使计算结果的幅值衰减和相速度发生变化，其作用相当于流动和相速度发生变化，其作用相当于流动中的物理耗散和弥散，这种虚假的物理中的物理耗散和弥散，这种虚假的物理效应称作数值耗散和数值弥散。效应称作数值耗散和数值弥散。n例如：对流方程描述的是既无耗散又无例如：对流方程描述的是既无耗散又无弥散的流体运动弥散的流体运动 n格式格式的等价方程 n截断面误差 n所以n代入等价方程n忽略高阶小量得 n二阶数值耗散系数 n方程适定性的要求方程适定性的要求n差分格式的适定性的必要条件为：差分格式的适定性的必要条件为：n与该格式的稳定性条件相同，可见通过上述分析出差与该格式的
19、稳定性条件相同，可见通过上述分析出差分格式稳定性的一些必要条件，这种方法称之为分格式稳定性的一些必要条件，这种方法称之为方法，可以用于方程组的非线性问题的分析。方法，可以用于方程组的非线性问题的分析。n上述是分析了截断误差的主要部分，如果将后面的高上述是分析了截断误差的主要部分，如果将后面的高阶小量展开后可知，该格式同样存在弥散，但与耗散阶小量展开后可知，该格式同样存在弥散，但与耗散相比要小。一般分析主要部分的数值效应而略去高阶相比要小。一般分析主要部分的数值效应而略去高阶小量。小量。同理有同理有n格式等价方程n二阶数值耗散系数 n格式等价方程n二阶数值耗散系数 n蛙跳格式等价方程n三阶弥散系
20、数 四、混淆误差四、混淆误差 n空间离散无法辨认小于空间离散无法辨认小于的波长，这的波长，这时小于时小于的短波分量会补充到长波分的短波分量会补充到长波分量中，从而使长波分量发生畸变，甚至量中，从而使长波分量发生畸变，甚至可能引起计算的不稳定。这种使短波分可能引起计算的不稳定。这种使短波分量补充到长波分量中，从而使长波分量量补充到长波分量中，从而使长波分量发生畸变而导致的误差称为混淆误差。发生畸变而导致的误差称为混淆误差。因此在具体问题中选取因此在具体问题中选取时，应保证波时，应保证波的主要分量变得充分的主要分量变得充分“长长”，即充分大。，即充分大。n函数函数 n导数导数n差商差商n差商与导数的比值为衰减比：n长波有，衰减系数：n因为，所以长波幅值很小，差商是微商好的近似，且当时 n考虑可辨认的短波。如考虑可辨认的短波。如，则，则，衰减比为，衰减比为，这时差商带来了很大的误差。，这时差商带来了很大的误差。n对于和空间步长对于和空间步长接近的短波，差商无接近的短波，差商无法近似导数。法近似导数。此课件下载可自行编辑修改，仅供参考！此课件下载可自行编辑修改，仅供参考！感谢您的支持，我们努力做得更好！谢谢感谢您的支持，我们努力做得更好！谢谢