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

§5.5关于对流项离散格式假扩散特性的讨论一、假扩散的含义本来含义：对流-扩散方程中一阶导数项的离散格式的截断误差小于二阶而引起较大数值计算误差。分析：纯对流方程，显式，FUD*第30页，共59页，星期日，2025年，2月5日一、假扩散的含义（续）假扩散系数*第31页，共59页，星期日，2025年，2月5日一、假扩散的含义（续）拓宽含义：由以下三种原因引起的数值计算误差。1、非稳态项或对流项采用一阶截差的格式2、流动方向与网格线呈倾斜交叉（多维问题）。3、建立差分格式时没有考虑到非常数源项的影响。*第32页，共59页，星期日，2025年，2月5日二、由于一阶导数截差阶数低而引起的假扩散一维无源项稳态模型方程*第33页，共59页，星期日，2025年，2月5日二、由于一阶导数截差阶数低而引起的假扩散（续）一维非稳态对流问题流向扩散*第34页，共59页，星期日，2025年，2月5日三、流速与网格线倾斜交叉引起的假扩散物理问题不考虑扩散作用的结果交叉扩散*第35页，共59页，星期日，2025年，2月5日四、由非常数源项引起的假扩散带源项的模型方程两点边值问题*第36页，共59页，星期日，2025年，2月5日五、低阶格式引起显著数值计算误差的例子Smith-Hutton问题出口截面结果乘方格式QUICK格式*第37页，共59页，星期日，2025年，2月5日五、低阶格式引起显著数值计算误差的例子(续)FUDHDPLDQUICK方腔自然对流宽高比：33Grashof数：9500*第38页，共59页，星期日，2025年，2月5日§5.6可以克服或减轻假扩散的格式或方法一、采用高阶格式以有效地克服流向扩散1、二阶迎风格式（SUD）2、三阶迎风格式（TUD）3、QUICK格式4、SGSD格式二、减轻或克服交叉扩散的方法1、对一阶迎风格式采用有效扩散系数2、采用自适应网格技术3、采用斜中心差分*第39页，共59页，星期日，2025年，2月5日1、二阶迎风格式（SUD）二阶精度，绝对稳定*第40页，共59页，星期日，2025年，2月5日2、三阶迎风格式（TUD）三阶精度，条件稳定*第41页，共59页，星期日，2025年，2月5日3、QUICK格式三阶精度，条件稳定三阶迎风格式修正系数为1/6*第42页，共59页，星期日，2025年，2月5日传热与流体流动的数值计算对流扩散方程的离散格式*第1页，共59页，星期日，2025年，2月5日§5.1对流项离散格式的重要性

及两种离散方式一、对流项离散格式的重要性1、数值解的准确性（假扩散）2、数值解的稳定性3、数值解的经济性二、构造离散格式的两种方式1、Taylor展开法2、控制容积积分法两种定义截差阶数一致，但截差首项系数有所不同。*第2页，共59页，星期日，2025年，2月5日§5.2对流项的中心差分与迎风格式一、一维对流-扩散问题模型方程的精确解边界条件：*第3页，共59页，星期日，2025年，2月5日一、一维对流-扩散问题模型方程的精确解（续）Peclet数：Pe表示对流与扩散作用的相对大小。00*第4页，共59页，星期日，2025年，2月5日二、对流项的中心差分对方程采用控制容积积分法记：F=ru通过界面的流量。界面上单位面积扩散阻力的倒数(扩导)。*第5页，共59页，星期日，2025年，2月5日二、对流项的中心差分（续）在数值计算过程中，如果连续性方程始终得到满足，则：在求解过程中，始终保持连续性方程满足非常重要。常物性条件下均分网格：*第6页，共59页，星期日，2025年，2月5日二、对流项的中心差分（续）例：在一维模型方程离散求解的均分网格中，已知fW=100，fE=200。试对PD=0，1，2及4四种情况按中心差分格式计算fP之值。负系数会导致物理上不真实的解。*第7页，共59页，星期日，2025年，2月5日三、对流项的迎风格式Taylor展开法控制容积积分法e界面w界面*第8页，共59页，星期日，2025年，2月5日三、对流项的迎风格式（续）