传热学第五章2
- 格式:ppt
- 大小:1.11 MB
- 文档页数:58
主讲陶文铨西安交通大学能源与动力工程学院热流中心CFD-NHT-EHT CENTER2010年10月18日, 西安数值传热学第五章对流扩散方程的离散格式(2)对流项离散格式的重要性及两种离散方式5.5.1假扩散的含义与成因5.5.2一阶截差格式引起严重假扩散举例1.本来的含义2.扩充的含义3.Taylor 展开法的分析5.5关于假扩散的讨论5.5.3网格倾斜交叉引起的计算误差5.5.4 非常数源项引起的假扩散5.5.5 两个名例以一维非稳态纯对流过程为例俩分析,其中有两n nφφ2(,O x φΔΔ其中关于时间的二阶导数项可做如下变化:时才没有这部分的计算误差。
2. 扩充的含义现有文献中常常将较大的计算误差都称为假扩散,大致有以下几项原因:(1) 一阶导数的一阶截差格式;(2) 流动方向与网格线呈倾斜交叉;(3) 离散格式未计及非常数源项的影响。
5.5.2一阶截差格式引起严重假扩散举例1.一维稳态对流扩散问题对流项用FUD,扩散项用CD,当Pe较大时,数值计算结果严重偏离精确解。
Physically plausible solution纯对流传递纯对流传递由离散方程:1n−1此时只有对流,没有扩散!时则有严重假扩散!0.8C =0.8C =当时,产生了严重的扩散作此种误差称为流向假扩散Γ≠Γ气流01. 设UE对P 控制容积,有2. 设控制容积,此时:计算误差纯对流传递三个对流问题的归纳这就是假扩散纯对流传递3)网格倾斜交叉引起的计算误差E冷热流体之间产生了温度均匀化的过程,即交叉5.5.5 已知流场计算温度场232(1),2(1)u y x v x y =−=−−参考解xT严重假扩散2) Leonard细高方腔中的自然对流换热5.6.1采用高阶格式克服流向假扩散5.6可以克服或减轻假扩散的格式与方法5.2.2 克服、减轻交叉假扩散的方法1. 采用二阶迎风2.采用三阶迎风3. 采用QUICK 格式1. 采用有效扩散系数2.采用自适应网格4. 采用SGSD 格式可以克服或减轻假扩散的格式与方法相当于界面上的中心差分)W WWxφ+Δ如型线上凹,则(2) FVM向上游取两点定义界面插值2.采用三阶迎风展开定义-一阶导数的三阶偏差分格式3. 采用定义-界面的插值在中心差分基础上考虑曲中心差分插值率修正?需要满足两个条件:插值的正确修正:相邻(2)0W PE φφφ−+<型线下凹8Cur −对e-界面u e 小于零时,取,,W P φφφu e 大于零时,取怎样相邻的三点?QUICK(2)e φφ=1/2w i φφ−=有:4. 采用CD条件稳定,但没有二阶假扩散;二阶迎风绝对稳定,组合起来,但是:如何确定值,特别是如何由计算结果来5. 高阶格式实施中的问题f u f计算边界:固o2) 代数方程的求解:等时,5.6.2用减小扩散系采用自适应网格(以减轻流5.7 对流-扩散方程离散形式稳定性分析5.7.1 数值计算中常见的三种不稳定性5.7.2 分析对流项格式不稳定性的“符号不变原则”5.7.3 稳定性分析结果讨论5.7.4 对流项格式问题讨论小结2.“符号不变”原则的基本思想3. “符号不变”原则的实施步骤4. “符号不变”原则的实施例子1. 研究背景扩散方程离散形式稳定性分析也会产生振荡的解,称为对流项离散格式的不稳态定性,研究目的是,找出产生振荡的临界Peclet 数。