对流扩散方程的离散格式.ppt

《对流扩散方程的离散格式.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、§5.1对流项离散格式的重要性及两种离散方式一、对流项离散格式的重要性1、数值解的准确性（假扩散）2、数值解的稳定性3、数值解的经济性二、构造离散格式的两种方式1、Taylor展开法2、控制容积积分法两种定义截差阶数一致，但截差首项系数有所不同。1§5.2对流项的中心差分与迎风格式一、一维对流-扩散问题模型方程的精确解边界条件：2一、一维对流-扩散问题模型方程的精确解（续）Peclet数：Pe表示对流与扩散作用的相对大小。003二、对流项的中心差分对方程采用控制容积积分法记：F=ru通过界面的流量。界面上单位面积扩散阻力的倒数(扩导)。4二、对流项的中心差分（续）在

2、数值计算过程中，如果连续性方程始终得到满足，则：在求解过程中，始终保持连续性方程满足非常重要。常物性条件下均分网格：5二、对流项的中心差分（续）例：在一维模型方程离散求解的均分网格中，已知fW=100，fE=200。试对PD=0，1，2及4四种情况按中心差分格式计算fP之值。负系数会导致物理上不真实的解。6三、对流项的迎风格式Taylor展开法控制容积积分法e界面w界面7三、对流项的迎风格式（续）e界面w界面8三、对流项的迎风格式（续）迎风格式离散形式：9四、中心差分与一阶迎风格式的讨论1、对流项中心差分在不发生振荡的参数范围内，比一阶迎风格式的误差更小。2、一阶迎风

3、格式离散方程系数永远大于零，不会引起解的振荡，得到物理上看似合理的解。3、一阶迎风格式截差阶数低，除非采用相当密的网格，否则计算结果的误差较大。4、一阶迎风格式的启示：应当在迎风方向取更多的信息构造格式，更好地反映对流过程的物理本质。5、在调试程序或计算的中间过程仍可以采用一阶迎风格式。10§5.3对流-扩散方程的混合格式及乘方格式一、系数aE与aW之间的内在联系aE(i)与aW(i+1)共享同一个界面。对流项中心差分：对流项一阶迎风：11二、混合格式(Spalding,1971)12三、指数格式利用精确解得到相邻节点间符合精确解的关系式。13三、指数格式（续）14四

4、、乘方格式(Patankar,1979)15五、5种3点格式系数汇总只需给出定义式16§5.4对流-扩散方程5种3点格式系数特性的分析总通量密度J：单位时间内、单位面积上由扩散及对流作用而引起的某一物理量的总转移量。一、通量密度及其离散表达式17一、通量密度及其离散表达式（续）J*的离散表达式：BehindAhead界面后的项界面前的项以坐标轴正方向为依据的“前”、“后”。18二、系数A、B间的关系1、和差特性当时，界面上的扩散通量为零，于是：19二、系数A、B间的关系（续）2、对称特性坐标系I：坐标系II：因为：于是：20二、系数A、B间的关系（续）指数格式系数A、

5、B间的关系21三、系数特性的重要推论和差特性：对称特性：重要推论：对5种3点格式的任何一种，若在PD>0时，A(PD)的计算式为已知，则在的范围内A(PD)、B(PD)的计算式均可得出。22三、系数特性的重要推论（续）证明：23四、aE、aW的通用表达式J通量密度守恒方程24四、aE、aW的通用表达式（续）不同格式的区别仅在于的计算式不同。25五、5种3点格式的26五、5种3点格式的（续）271、从一维到多维的推广在每一个坐标方向上都按一维问题处理。2、所得出的系数表达式便于编制通用性程序，由专用模块处理。3、利用aE(i)与aW(i+1)间的关系，可以节省计算系数的

6、工作量。六、关于格式定义与系数特性的说明28§5.5关于对流项离散格式假扩散特性的讨论一、假扩散的含义本来含义：对流-扩散方程中一阶导数项的离散格式的截断误差小于二阶而引起较大数值计算误差。分析：纯对流方程，显式，FUD29一、假扩散的含义（续）假扩散系数30一、假扩散的含义（续）拓宽含义：由以下三种原因引起的数值计算误差。1、非稳态项或对流项采用一阶截差的格式2、流动方向与网格线呈倾斜交叉（多维问题）。3、建立差分格式时没有考虑到非常数源项的影响。31二、由于一阶导数截差阶数低而引起的假扩散一维无源项稳态模型方程32二、由于一阶导数截差阶数低而引起的假扩散（续）一维

7、非稳态对流问题流向扩散33三、流速与网格线倾斜交叉引起的假扩散物理问题不考虑扩散作用的结果交叉扩散34四、由非常数源项引起的假扩散带源项的模型方程两点边值问题35五、低阶格式引起显著数值计算误差的例子Smith-Hutton问题出口截面结果乘方格式QUICK格式36五、低阶格式引起显著数值计算误差的例子(续)FUDHDPLDQUICK方腔自然对流宽高比：33Grashof数：950037§5.6可以克服或减轻假扩散的格式或方法一、采用高阶格式以有效地克服流向扩散1、二阶迎风格式（SUD）2、三阶迎风格式（TUD）3、QUICK格式4、SGSD格式二、