第五章数值传热学

合集下载

数值传热学(课件)

数值传热学(课件)

02 数值传热学的基本原理
控制方程
控制方程
数值传热学的核心是求解控制方 程,这些方程描述了热量传递过 程中的物理规律。
偏微分方程
控制方程通常以偏微分方程的形 式给出,包含了温度、时间、空 间等变量的变化关系。
初始条件和边界条

为了求解控制方程,需要给出初 始条件和边界条件,这些条件限 定了问题的解的范围。
详细描述
传热过程模拟是数值传热学的另一重要应用,通过建立传热过程的数学模型,可以模拟物体内部的温 度分布和热量传递过程。这对于能源、化工、电子等领域中的热工设备设计和优化具有重要意义。
04 数值传热学面临的挑战与 解决方案
计算精度与稳定性问题
总结词
计算精度和稳定性是数值传热学中的核心问题,直接关系到模拟结果的准确性和可靠性。
详细描述
多尺度问题要求数值方法能够捕捉到不同尺度的物理现象,并准确地将它们联系起来。 这需要发展具有多尺度分辨率的数值方法,如多重网格法、谱方法和自适应网格法等。
非线性问题
总结词
非线性问题在传热过程中广泛存在,如 流动、相变和化学反应等,给数值模拟 带来很大难度。
VS
详细描述
非线性问题需要数值方法能够处理高度非 线性的物理方程,并能够准确地捕捉到非 线性现象。这需要发展高效的数值算法, 如有限元法和有限体积法等,同时还需要 考虑非线性问题的特殊性质,如初始条件 和边界条件等。
02
它涉及传热学的基本原理、数学 建模、数值计算和计算机技术等 多个领域,是计算流体动力学和 计算传热学的重要组成部分。
数值传热学的重要性
随着科技的发展,传热问题在能源、 环境、航空航天、化工等领域越来越 突出,数值传热学的应用也越来越广 泛。

数值传热学 -回复

数值传热学 -回复

数值传热学 -回复
数值传热学(Numerical Heat Transfer)是一门研究热传递现象的学科,通过数值模拟和计算方法来分析热传导、对流和辐射等传热过程。

本文将介绍数值传热学的基本原理、方法和应用。

1. 基本原理
数值传热学基于传热学原理和计算数学方法,将传热过程建模为数学方程,并通过数
值方法求解这些方程,从而得到热传递的数值解。

主要的传热模型包括热传导、对流和辐
射传热。

2. 数值方法
数值传热学常用的方法包括有限差分法、有限元法和边界元法等。

有限差分法是最常
用的方法之一,将传热区域离散化为网格,通过差分近似计算网格点上的温度或热流量。

有限元法则是另一种常用的方法,将传热区域划分为元素,通过建立元素之间的关系来计
算温度场或热流场。

边界元法则是将问题转化为边界上的积分方程,通过求解积分方程得
到温度场或热流场。

3. 应用领域
数值传热学在各个领域都有广泛的应用。

在工程领域,数值传热学用于优化热交换器
的设计、预测电子器件温度分布、模拟流体在管道内的传热过程等。

在材料科学领域,数
值传热学用于研究材料的导热性能、相变过程以及焊接和烧结等工艺。

在能源领域,数值
传热学用于分析太阳能热收集器的性能、燃烧过程中的传热机制等。

通过数值传热学的研究,我们可以更加深入地了解热传递过程,并可以通过数值模拟
方法来预测和优化热传递的效果。

数值传热学也为各个领域的工程和科学研究提供了重要
的工具和方法。

通过不断的发展和创新,数值传热学将进一步推动热传递理论和应用的发展。

数值传热学第五章-数值计算

数值传热学第五章-数值计算
即, F 2D时,有可能使 aE 或aW 为负 产生不切实际的结果
2016-1-30
太 原 理 工 大 学
9 /70
Thermal
这就是中心差分格式求解对流换热问题时仅限于低
Fw Fe Re(低的F/D)的原因 . aP Dw De aW aE Fe Fw 2 2
若E 200, W 100 P 50
两个值均不 符合实际
若E 100, W 200 P 250
aW Dw Fw 2 1 2 3
aE De Fe 2 1 2 1 违背了正系数规则
aP aE aW 1 3 2, 而 anb 1 3 4 这样,aP anb , 违反了斯卡巴勒准则( 主对角占优)
ui 0
i-1 W
w
二类迎风格式) 控制容积界面上值的规定: 界面上的值等于界面上风侧 网格节点上的值。
e P e E
Fe 0 Fe 0
i
P
e
i+1
E
ui 0
类似地,w界面上
w W Fw 0 w P Fw 0
上述条件语句紧凑格式的写法:
Γe、Γw可以用算术平均法或调和平均法求得。
定义:
D 扩散传导性. x
Thermal
F u 对流或流动强度,可正 、可负,由流动方向定
整理后的离散化方程 其中:
2016-1-30
a p P aE E aWW
Fw 2 太 原 理 工 大 学 aW Dw
7 /70
太 原 理 工 大 学
3 /70
Thermal

数值传热学的通用方程

数值传热学的通用方程

数值传热学的通用方程数值传热学的通用方程引言:传热学是研究热量在物体内传递的学科,它在实际生活中具有广泛的应用。

数值传热学是传热学的一个重要分支,借助数值计算方法和计算机模拟,能够更准确地预测和模拟热量的传递过程。

在数值传热学中,通用方程是一种重要的工具,它能够描述和计算物体内热量的传递方式。

本文将以数值传热学的通用方程为主题,通过分析其深度和广度,以全面评估和解释这一概念。

一、数值传热学的基础概念1.1 热量传递的三种方式热量传递有三种方式:传导、对流和辐射。

传导是指热量通过物质的直接接触和振动传递,对流是指热量通过流体的传输,辐射是指热量通过电磁波辐射传递。

这三种方式在不同的情况下起着不同的作用,同时它们也相互影响和耦合。

1.2 数值计算方法在传热学中的应用数值计算方法是数值传热学的核心工具,它可以通过数学模型和离散计算,模拟和预测物体内热量的传递过程。

常用的数值计算方法有有限元法、有限差分法和有限体积法等。

通过这些方法,我们可以更准确地计算和研究热量的传递规律。

二、数值传热学的通用方程2.1 传热方程的基本形式传热方程是描述热量传递过程的数学方程,它以物体内部的温度分布、热流和热导率等参数为基础,通过各种数学方法和推导,得到不同传热方式下的通用方程。

2.2 热传导方程热传导方程是描述热量通过传导方式传递的方程。

在传热过程中,热量会从高温处传向低温处,而传热率又与温度梯度和材料的热导率成正比。

热传导方程能够计算和描述热量在物体内部的传递过程,为热传导问题的分析和计算提供了基础。

2.3 流体传热方程流体传热方程是描述热量通过对流方式传递的方程。

流体传热过程中,流体的流动状态和温度梯度会影响热量的传递速率。

流体传热方程能够计算和描述流体内部的热量传递过程,对于流体传热问题的研究和分析具有重要意义。

2.4 辐射传热方程辐射传热方程是描述热量通过辐射方式传递的方程。

辐射传热过程中,热量通过电磁波的辐射传输,与物体的温度和辐射特性有关。

【免费下载】数值传热学第五章作业

【免费下载】数值传热学第五章作业

5-2解:根据课本p158式(5—1a )得一维稳态无源项的对流-扩散方程如下所示: (取常物性)22x x u ∂∂Γ=∂∂φφρ边界条件如下:L L x x φφφφ====,;,00由(5—2)得方程的精确解为: 11)/(00--=--⋅Pe L x Pe L e e φφφφΓ=/uL Pe ρ将分成15等份,有:L ∆=P Pe 15对于中心差分、一阶迎风、混合格式和QUICK 格式分别分析如下:1)(CD)中心差分节点离散方程: 2)5.01()5.01(11-∆+∆++-=i i i P P φφφ10,2 =i 2)一阶迎风节点离散方程: ∆-∆++++=P P i i i 2)1(11φφφ10,2 =i 3)混合格式当时,节点离散方程:,1=∆P 2)5.01()5.01(11-∆+∆++-=i i i P P φφφ10,2 =i 当时,节点离散方程: , 10,5=∆P 1-=i i φφ10,2 =i 4)QUICK 格式,节点离散方程: , ⎥⎦⎤⎢⎣⎡--++++++=+-∆∆-∆∆+∆)336(81221211111i i i i i i P P P P P φφφφφφ2=i , ⎥⎦⎤⎢⎣⎡---++++++=+--∆∆-∆∆+∆)35(812212112111i i i i i i i P P P P P φφφφφφφ2≠i用matlab 编程如下:(本程序在x/L=0-1范围内取16个节点进行离散计算,假设y(1)= =0,y(16)==1,程序中Pa 为,x 为题中所提的x/L 。

由于本程序假设y(1)=0φL φ∆P =0,y(16)==1,所以)0φL φy y y y y y L =--=--=--010)1()16()1(00φφφφPa=input('请输入Pa=')x=0:1/15:1Pe=15*Pa;y=(exp(Pe*x)-1)/(exp(Pe)-1)plot(x,y,'-*k') %精确解hold ony(1)=0,y(16)=1;for i=2:15y(i)=((1+0.5*Pa)*y(i-1)+(1-0.5*Pa)*y(i+1))/2;endplot(x,y(1:16),'-or') %中心差分hold onfor i=2:15y(i)=((1+Pa)*y(i-1)+y(i+1))/(2+Pa);endplot(x,y(1:16),'-.>g') %一阶迎风hold onfor i=2:15if Pa==1y(i)=((1+0.5*Pa)*y(i-1)+(1-0.5*Pa)*y(i+1))/2;elsey(i)=y(i-1)endendplot(x,y(1:16),'-+y') %混合格式hold onfor i=2:15if i==2y(i)=y(i+1)/(2+Pa)+(1+Pa)*y(i-1)/(2+Pa)+(Pa/(2+Pa))*(6*y(i)-3*y(i-1)-3*y(i+1))/8 elsey(i)=y(i+1)/(2+Pa)+(1+Pa)*y(i-1)/(2+Pa)+(Pa/(2+Pa))*(5*y(i)-y(i-1)-y(i-2)-3*y(i+1))/8 endendplot(x, y(1:16),'-<b') %QUICK 格式hold onlegend('精确解','中心差分','一阶迎风','混合格式','QUICK 格式')运行结果如下图所示:当 :1=∆P当:5=∆P当:10=∆P5-3 解:根据课本式(5-19)得:乘方格式:⎪⎪⎩⎪⎪⎨⎧<-≤≤--+≤≤->=∆∆∆∆∆∆∆∆10,010,)1.01(100,)1.01(10,055P P P P P P P P D a e E 当时有:1.0=∆P 951.0)1.01.01()1.01(55=⨯-=-=∆P D a e E 301.0/3)()()()()()(===Γ=Γ=∆ee e e e e e e e P u x u u x D ρδρρδ5297.2830951.0951.0=⨯==e E D a 由系数关系可得:∆=-P D a D a e E w W 53.3130)951.01.0((=⨯+=⨯+=∆w e E W D D a P a根据式(5-51g )得: 205.01.010=⨯=∆∆=tx a P p ρ根据式(4-12)得: (本题方程中无源项)0P W E P a fa fa a ++=当采用隐式时,则得到:1=f 0597.62253.315297.280=++=++=P W E P a fa fa a 即:时,,,,1.0=∆P 5297.28=E a 53.31=W a 20=p a 0597.62=P a 当时,按照以上算法得出:10=∆P ,, , 0=E a 3=W a 20=p a 5=P a。

数值传热学

数值传热学

数值传热学
t为了更好地理解热学中的非稳态传热现象,需要对其进行数值模拟,在数值传热学方法中,有一种方法叫做有限元方法,它是一种基于网格方法的非线性有限元方法。

ttt在研究和处理复杂工程问题时,为简化计算机求解代价高的无限大规模的实际物理问题,常采用网格技术,对复杂的多相流动或物体的运动状态进行模拟,并将该计算过程和成果称之为“数值模拟”。

ttt在应用数值传热学方法的过程中要注意这样几点:一是网格划分、初始条件及边界条件的选取要适当二是系统初始化要合理三是尽可能使所有的网格之间相互独立四是保证结果的重现性五是不要
忽视分辨率的概念六是分析与综合要紧密联系起来七是数值计算过
程要符合数学规,使输出的数据便于人们分析比较八是在数值计算过程中若发现新的或难以理解的情况或事件,应记录下来,待分析完后再去验证九是对所得到的结果要进行认真检查。

t有限单元法是在有限空间或无限体积中把某些大块区域作为节点,其他区域为单元,用有限个节点(单元)组成有限个相互连接的单元链。

这种方法将无限的区域离散化成有限个单元,在每个单元内假定一定的约束条件和单元本身的物理属性。

网格在三维空间中的布置形式,可以由连续函数来描述。

有限元法通过把物理问题分解成许多微小的单元,然后按照一定的节点连接关系进行组合,并假定这些单元遵循各自的约束条件。

当计算机通过网络将数据存入存储器中时,有限元法就得到了充分发挥,可以利用计算机快速运算获得高精度的解。

但由于有限元法是一
种离散化方法,因此如果计算时出现局部收敛性差的问题,很可能导致整个求解过程失败,从而影响最终结果的准确性。

数值传热学

数值传热学

数值传热学数值传热学(numerical heat transfer)数值传热学,又称计算传热学,是指对描写流动与传热问题的控制方程采用数值方法,通过计算机求解的一门传热学与数值方法相结合的交叉学科。

数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场(如速度场,温度场,浓度场等),用一系列有限个离散点上的值的集合来代替,通过一定的原则建立起这些离散点变量值之间关系的代数方程(称为离散方程)。

求解所建立起来的代数方程已获得求解变量的近似值。

数值传热学(numerical heat transfer)数值传热学,又称计算传热学,是指对描写流动与传热问题的控制方程采用数值方法,通过计算机求解的一门传热学与数值方法相结合的交叉学科。

数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场(如速度场,温度场,浓度场等),用一系列有限个离散点上的值的集合来代替,通过一定的原则建立起这些离散点变量值之间关系的代数方程(称为离散方程)。

求解所建立起来的代数方程已获得求解变量的近似值。

数值传热学常用的数值方法1.有限差分法历史上最早采用的数值方法,对简单几何形状中的流动与换热问题最容易实施的数值方法。

其基本点是:将求解区域中用于坐标轴平行的一系列网格的交点所组成的点的集合来代替,在每个节点上,将控制方程中每一个导数用相应的差分表达式来代替,从而在每个节点上,形成一个代数方程,每个方程中包括了本节点及其附近一些节点上的未知值,求解这些代数方程就获得了所需的数值解。

2.有限容积法将所计算的区域划分成一系列控制容积划分为一系列控制容积,每个控制容积都有一个节点做代表。

通过将守恒型的控制方程对控制容积坐积分导出离散方程。

在导出过程中,需要对界面上的被求函数本身及其一阶导数的构成做出假定,是目前流动与换热问题的数值计算中应用最广的一种方法。

3.有限元法把计算区域划分为一系列原题(在二维情况下,元体多为三角形或四边形),由每个元体上去数个点作为节点,然后通过对控制方程做积分来获得离散方程。

传热学第五章

传热学第五章
' u1 " u1
例2:流体外掠平板对流换热边界层温度场相似问题 温度沿 x、y 方向变化 如果在空间 对应点上: 几何相 似倍数
' x1 " x1
=
' r2 " r2
=
r3'
" r3
= .... =
' um " um
R'
几何相 似倍数
=
' u2 " u2
=
' u3 " u3
= .... =
' u max
∂x
∂y
∂y
Cu Cl Ca
” ∂t” ” ∂t” ” ∂ 2 t” =a u +v 2 ∂x” ∂y” ∂y”
hl — — 努谢尔特数( Nusslet ) λ ρul ul = Re = — 雷诺数 ( Reynolds ) η ν Nu = ν a ∆p Eu = ρu 2 Pr = — — 普朗特数 (Prandtl) — — 欧拉数 (Euler) ul — 贝克利数 (Peclet)20 a
17
∂u” ∂v” + =0 ∂x” ∂y”
Cu Cl u 'l ' υ " =1 ⇒ =1 Re ' = Re" Cυ υ ' u "l " C∆p Cl ∆p ' u 'l ' ρ "u "2 υ " = 1 ⇒ ' '2 ' =1 C ρ Cu Cυ ρ u υ ∆p" u "l "
Eu ' Re ' = Eu " Re "

数值传热学

数值传热学

数值传热学数值传热学是一门研究如何采用计算机技术模拟传热过程的学科。

它的出现,使得传热学在进行理论分析和数值计算方面更加具有实际意义。

以前要进行传热问题的求解,需要有丰富经验的工程师去积累相应的数据,而且也只能做到相对比较好的效果,这对于我们来说也并非难事,但是就目前情况来看,人工建立一个数值模型对于复杂的传热过程进行求解将会变得更加困难。

而且传统的工业过程模拟中大多都是依靠经验,但是这种经验往往又是片面的、偶然的。

因此,建立一个完整的数值传热学体系就成了当务之急。

然而数值传热学作为一门年轻的学科,它的许多思想都源于传热学的实际应用。

因此,对于传热学基础理论知识的掌握以及综合应用能力的培养对于该课程学习至关重要。

在教学过程中,除了注重学生自身素质的培养外,还应该结合教学内容,创新教学方法,在充分调动学生主观能动性的基础上,发挥他们的创造性思维,让学生参与其中。

另外,还可以借助多媒体教学等现代化教学手段,提高课堂效率,增强教学效果。

从而促进学生对于课程的理解,提升教学质量。

传统的传热学模型很难完全适用于现代化传热分析与设计。

针对于传热模型方面,首先,需要增加数值传热学模型,在原有的基础上引入新的概念和规则;其次,模型的编制需要精确考虑每一个单元模块之间的关联,在保证各个模块都能够单独准确地计算出结果的同时,还必须将他们联系起来。

例如,对于燃烧室内传热问题,由于温度的分布情况是非常复杂的,因此我们需要构造适当的网格进行相应的计算,将所有网格划分成细小的区域,再逐步地建立起燃烧室内温度场的整体结构,进而达到一个有效的、清晰的计算结果。

除此之外,计算结果的收敛速度和计算的精确度也是影响分析效率的两个主要因素。

但是随着对于这些相关领域研究人员的不断增多,一些技术已经可以得到大幅度的改善,甚至部分可以直接用于商业用途。

而计算流体力学( CFD)就是在计算机运算能力不断增强的基础上逐渐形成的一门新兴学科,它可以通过在计算机上建立一些专门的数学模型,利用计算机仿真软件进行求解,最终获取相应的物理图像或者曲线,帮助工程师快速地找到解决方案。

化工原理(第五章传热第五节)

化工原理(第五章传热第五节)

吉 首 大 学
流体在管束外横掠流动
化 工 原 理 由于各排的给热系数不同,则整个管束的平均给热系数应按 下式求出: a1A1+ a2A2 + a3A3 + … am = A1+A2 + A3 + … 式中:A1、A2、A3……分别为第一排,第二排,第三排…… 的传热面积; a1 、 a2 、 a3……分别为第一排,第二排,第三排…… 的传热系数。
d A2 π d 2dl d 2 = = d Am π d mdl d m
吉 首 大 学
1 = d2 + b d2 + 1 a2 a1 d1 λ dm K2
当间壁为平壁,或管壁很薄或管径较大时,dA1 、dA2 、dAm 和 dA 相等或近似相等,则: 1 = 1 + b + 1 a2 a1 λ K2
Q QR
Q A + QR + QD = Q Q A QR QD + + =1 Q Q Q
QA
吉 首 大 学
QD
A+ R + D =1
A、R 和 D 分别为物体吸收率、反射率和透过率。 单色吸收率、反射率和透过率
a(λ , T ) + r (λ , T ) + d (λ , T ) = 1
基本概念
化 工 原 理 黑体(绝对黑体):能将辐射能全部吸收的物体,即 A=1, R=D=0。自然界中并不存在绝对黑体,例如没有光泽的黑墨 表面,其吸收率 A=0.96~0.98,定义黑体的目的是为了在计 算中确定一个比较的标准。 镜体(绝对白体):能将辐射能全部反射的物体,即 R=1, A=D=0。自然界中也不存在绝对镜体,例如表面抛光的铜, 其反射率 R=0.97。 0.97 透热体:辐射能全部透过的物体,即D=1, A=R=0。例如对 称双原子气体 O2、N2、H2 等都是透热体。 灰体:能够以相等的吸收率吸收所有波长辐射能的物体。灰 体也是理想物体,其特点为:吸收率 A 与波长无关;为不透 热体 (A+R=1)。工业上常见的固体材料均可视为灰体。

传热学-第五章

传热学-第五章
第五章 对流换热 2
2 对流传热的特点 ① 导热与热对流同时存在的复杂热传递过程。
② ③
必须有直接接触(流体与壁面)和宏观运动;也必须有温差。 由于流体的粘性和受壁面摩擦阻力的影响,紧贴壁面处会形 成速度梯度很大的边界层。 3 对流传热的基本计算式
牛顿冷却式:
Φ hA(t w t ) W
对流传热微分方程组:(常物性、无内热源、二维、不可 压缩牛顿流体)
u v 0 x y
u u u p 2u 2u ( u v ) Fx ( 2 2 ) x y x x y v v v p 2v 2v ( u v ) Fy ( 2 2 ) x y y x y
Q = E + W
Q — Q导热 Q对流 Q内热源
E — U 热力学能 U K(动能) U P(位能)
W — 体积力(重力)作的功、表面力作的功 假设:(1)流体的热物性均为常量,流体不做功 (2)流体不可压缩 (3)一般工程问题流速低
W=0
(4)无化学反应等内热源
第五章 对流换热
3 密度 [kg m ] 热导率 [ W (m C) ] 2 比热容 c [J (kg C) ] 动力粘度 [ N s m ] 2 运动粘度 [m s] 体胀系数 [1 K ]

1 v 1 v T p T p
UK=0 Q内热源=0
21
Q导热 + Q对流 = U热力学能
Q导热
2t 2t 2 dxdy+ 2 dxdy x y
单位时间内、 沿 x 方向热对流传递到微元体的净热量:
" " " Qx Qx (ut) " " " Qx Qx dx Qx Qx dx dx c p dxdy x x x

数值传热第五章课件2陶文铨

数值传热第五章课件2陶文铨

主讲陶文铨西安交通大学能源与动力工程学院热流中心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 数。

第五章传热(本专业)(1)精品PPT课件

第五章传热(本专业)(1)精品PPT课件
0K,均存在辐射传热; ▪ 不需要任何中介; ▪ 传热过程中伴随能量 形式的转换。
三种传热方式的比较:
传导 对流 辐射
注:三种传热方式往往共存
5.1.3 工程上常用的换热方式
⒈混合式换热
冷热两种流体直接接触换热,如凉水塔,湿式混 合冷凝器。 优点:传热速度快、效率高、设备简单等。
2. 蓄热式换热
计算:
厚度为b 的无限大平壁,壁 面两侧温度t1、t2 ,t1>t2 , 取厚度为dx 的薄层,由傅
立叶定律:
q dt
dx dt q dx
对上式积分,积分限为:
t : t1 t2
x:0
λ取一平均值,视为常 数,积分得:
q
t1
t2
t
Q
qA
A
t1
t2
说明:
①将上式写成速率方程的一般形式为:
分率、分子量及导热系数。
气体的 导热系数:
1-水蒸气;2-氧;3-二氧化碳;4-空气; 5-氮;6-氩
5.2.4 平壁的稳定热传导
㈠单层平壁的稳定热传导
平壁模型:
▪ 平壁材质均匀,λ可视为常数;
▪ 平壁内只有一维温度梯度,导热方向垂直于壁面 ─等温面为平行于侧面的平面;
▪ 导热平壁的长和宽>>壁厚b ,忽略边缘热损失。
─等温面为与圆筒同心的圆筒面;
▪ 筒壁材质均匀,λ视为常数。
计算:
内、外半径r1、r2 , 内外壁温度t1、t2(t1>t2), A=2πrl,导热系数λ,由傅
立叶定律:
QArq2rLddrt
分离变量: dt Q dr
2L r
积分: d t2 t r2 Q dr
t1
r1 2L r

数值传热学第5章作业答案

数值传热学第5章作业答案

第5章作业答案5-2对于5种三点格式来说,一维对流扩散方程都是可以写成下列通用离散形式:P P E E W Wa a a φφφ=+ 其中: [](){}()[]{}()w e W E P w w w W e e e E F F a a a P P A D a P P A D a -++=+=-+=∆∆∆∆0,0,5种三点格式的()∆P A格式()∆P A迎风差分 1混合格式 []|5.01,0|∆-P 指数格式 ()()1exp -∆∆P P对网格Peclet 数为5,10的情形,应该得出如下图的结果,FUD 与混合格式没有振荡,而CD 和QUICK 均有,而且CD 比QUICK 更为严重。

5-3不同网格∆P 数下各系数计算结果如下∆P E aW a 0P a P a 0.1 28.53 31.53 2 62.05910 0 3255-5 四个节点之值如下一阶迎风 混合格式 乘方格式 二阶迎风(边界一阶) 二阶迎风(边界二阶)1φ 94.26 73.96 79.01 58.57 91.122φ 147.61 91.10 115.13 76.65 144.19 3φ 82.14 72.40 74.19 69.33 81.34 4φ 126.99 85.31 102.70 87.38 124.505-7不计扩散项,采用QUICK 离散i 控制容积的非稳态与对流项得:12117338n nn n n ni i i i i i x utφφφφφφ+--+--++∆=-∆ ((0)u >采用离散扰动分析法,对i+1得到扰动为78n i u t ρε∆,对i-1 得到扰动为38ni u t xε∆-∆,符号不变原则要求:0832≥∆Γ∆+∆∆-ninin i x t x t u εερερ,由此得:38≤=Γ∆∆P xu ρ5-9根据三阶迎风格式的定义:⎪⎩⎪⎨⎧<∆--+->∆+-+=∂∂-++--+0,62360,6632112211u x u xx i i i i i i i i φφφφφφφφφ仿照QUICK 格式,令三阶迎风格式的控制容积右界面上的值的形式为:⎪⎩⎪⎨⎧<+--+>+--+=0,220,22u a u a EEE P E P WP E E P e φφφφφφφφφφφ同理可以写出w φ的计算式。

数值传热学答案范文

数值传热学答案范文

数值传热学是热力学的重要分支之一,研究物质中热量的传递和分布规律。

与传统的实验方法相比,数值传热学采用计算机模拟技术,通过数学模型和计算实验方法,能够更加深入、系统地研究热传递现象的规律和特性,为工程设计和实际生产提供重要的技术支持。

数值传热学的本质是热传递方程的数值求解。

热传递方程是描述物质中热量传递和分布的方程,它包含了热传导、热对流和热辐射三种传热方式。

热传导是指热量沿着物质内部的温度梯度传递,主要发生在固体和液体中;热对流是指热量随物质的流动而传递,主要发生在液体和气体中;热辐射是指热量通过辐射传递,主要发生在光学和辐射热转换材料中。

通过数值方法求解热传递方程,可以得到物体的温度分布、热传递速率和热流密度等参数,为材料和工程设计提供准确的数据支持。

数值传热学的核心是数值方法,主要包括有限差分、有限元和边界元等方法。

有限差分法是一种利用离散化方法求解微分方程的数值方法,它将微分方程中的连续变量离散化,将求解微分方程转化为求解线性方程组。

有限元法是一种利用有限元逼近方法解决偏微分方程的数值方法,采用对物体进行简单的几何划分,将问题离散化,通过数学建模来表示物体的温度分布和热流密度分布。

边界元法是一种较新的有限元法补充,它能够快速解决边界值问题,并且可以减少问题的维数。

数值传热学的应用范围广泛,包括热工和物理问题的研究、能源系统分析和设计、建筑工程中的热传递和能源效率研究等。

例如,在太阳能发电系统设计中,数值传热学可以帮助设计人员确定集热器表面温度和吸收率等参数,提高太阳能效率并减少系统成本。

在建筑工程中,数值传热学可以帮助设计师分析建筑物的保温性能,合理评估保温材料的性能和使用效果,确保建筑节能和环保。

在机械加工领域中,数值传热学可以帮助工程师分析材料切削过程中的热量和温度分布,挑选适合材料和刀具的加工工艺,提高机械切削效率。

数值传热学是现代科学技术的重要分支之一,是研究物质中热传递和分布规律的重要工具。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

v 0, 0 v 0, 0
0
Fs W S ! Fw Fs , aPP aW W aSS , aP aW aS , P 2
12/60
三个对流问题的归纳 1) 一维稳态对流扩散问题 2) 一维非稳态对流 问题(Noye,1976) 此种误差称为流向假 扩散 (Streamwise
如型线上凹,则 ( P 2W WW ) 0 修正为正; 如型线下凹,则
(P 2W WW ) 0 修正为负。
23/60
(2) FVM定义-向上游取两点定义界面插值
w
1.5W 0.5WW , u 0 1.5P 0.5E , u 0
两种定义等价性的说明:
所以
w
E
0 0
来流速度为 U
P W !
11/60
上游温度一直保持到下游。
2. 设x-y坐标分别与来流呈45度倾角 仍采用一阶迎风格式,对P控制容积,此时:
aE De Fe ,0
2 U, 来流速度为 U u v 2 u 0, 0
0
w
aW DW FW ,0 u 0, 0 F aN Dn Fn ,0 aS Ds Fs ,0
说明:这里讨论的是二阶假扩散;即在二阶截断误差 的意义上,该格式实际模拟的是一个扩散-对流过程。 一般工程计算中讨论二阶假扩散已经能够满足计算精度 的要求。
6/60
x
2. 扩充的含义
现有文献中常常将较大的计算误差都称为假扩 散,大致有以下几项原因: (1) 一阶导数的一阶截差格式; (2) 流动方向与网格线呈倾斜交叉; (3) 离散格式未计及非常数源项的影响。 5.5.2 一阶截差格式引起严重假扩散举例 1. 一维稳态对流扩散问题 对流项用FUD,扩散项用CD,当Pe较大时,数值 计算结果严重偏离精确解。
所以上述一阶截差的离散方程在二阶截差的意义上 相当于模拟了一个对流扩散方程,而不是纯对流方程。
5/60
只有当 1 u t 0 时才没有这部分的计算误差。
ut 称为Courant数,以纪念德国数学家Courant。 x ux ut 2 )i ,n u )i ,n [ (1 )] 2 )i ,n O(x 2 , t 2 ) 2 t x x x
ui u )i (2i 1 3i 6i 1 i 2 ), u 0 x 6x
特点:利用了下游一个节点,提高了格式精度,但减 弱了稳定性。 (2) FVM定义-向上游取两 点,向下游取一个节点来定 义界面插值。
25/60
3. 采用QUICK格式 1) FVM定义-界面的插值在中心差分基础上考虑曲 率修正
1 E ) 0 Cur 使插值微减! 8 1 E ) 0 Cur 使插值微增! 8
27/60
怎样相邻的三点? (2) 为增加格式的稳定性要引入迎风思想:对e-界面 当界面流速ue大于零时,取 当界面流速ue小于零时,取
W , P ,E P ,E ,EE
ue大于零时,取
aE De A( Pe ) Fe ,0 , A( Pe ) 0,1 0.5 Pe
设网格Peclet数的改变通过改变扩散系数来实现,则
P e 2, aE 之值均相同,使得所有 P e 2 数值解均同。
16/60
5.5.5 两个名例 1. Smith-Hutton问题(1982)已知流场计算温度场
u 2 y (1 x 2 ), v 2 x(1 y 3 )
已知的流场
已知的进口温度分布
Tin ( x) 1 tanh[ (1 2 x)]

称为陡度系数
17/60
用二维对流扩散方程求解,对流项采用所研究格式:
T 参 考 解
乘方格式
x
QUICK格式
QUICK格式用20X20的节点数就可达到乘方格 式80X80的精度。
1 e w (1.5P 0.5W ) (1.5W 0.5WW ) dx x w x x x
e
3P 4W WW 2x
FVM定义的是计算区域中的 积分中值;FDM定义的是一 点的离散值。
24/60
2.采用三阶迎风 (TUD) (1) Taylor展开定义-一阶导数的三阶偏差分格式
20/60
5.6 可以克服或减轻假扩散的格式与方法 5.6.1 采用高阶格式克服流向假扩散 1. 采用二阶迎风 2.采用三阶迎风 3. 采用QUICK格式 4. 采用SGSD格式 5.2.2 克服、减轻交叉假扩散的方法 1. 采用有效扩散系数 2.采用自适应网格
21/60
5.6 可以克服或减轻假扩散的格式与方法 5.6.1 采用高阶格式克服流向假扩散 1. 采用二阶迎风-向上游两个点获取信息构成导数 (1) Taylor展开定义-一阶导数的二阶偏差分格式
数值传热学
第五章 对流扩散方程的离散格式(2)
主讲 陶文铨
西安交通大学能源与动力工程学院 热流科学与工程教育部重点实验室 2012年10月17日, 西安
1/60
5.5 关于假扩散的讨论 5.5.1 假扩散的含义与成因
1.本来的含义 2.扩充的含义 3.Taylor展开法的分析
5.5.2 一阶截差格式引起严重假扩散举例 5.5.3 网格倾斜交叉引起的计算误差 5.5.4 非常数源项引起的假扩散 5.5.5 两个名例
18/60ຫໍສະໝຸດ 2) Leonard问题(1996) 细高方腔中的自然对流换热
H / L 33
Gr gL3T
2
9500, Pr 0.71
32 129 4128
低阶格式无法分辨出小涡。
19/60
QUICK scheme PWL scheme 节点数增加到一定数目时乘方格式也能分辨出小 涡。格式间的差别随着网格的细化而逐渐消失
3/60
以一维非稳态纯对流过程为例俩分析,其中有两 个一阶导数项,均采用一阶截差的格式。
u t x
将 ,
n i 1 n 1 i
一阶截差格式离散
in1 in
t
u>0
u
in in 1
x
对点(i,n)做Taylor展开,代入上式:
2 1 2 n i t )i ,n t 2 )i ,n ... in t 2 t t
2 1 2 n n ) ... i [i x )i ,n x 2 i ,n 2 x x u x
4/60
t 2 ux 2 2 2 )i ,n u )i ,n ) ) ( , ) x O t 2 i ,n 2 i ,n 2 t 2 x t x 其中关于时间的二阶导数项可做如下变化:
纯对流传递
C 0.8
9/60
当Courant数小于1时,产生了严重的扩散作 用,将尖峰逐渐抹平。此种误差称为流向假扩散 (Streamwise false diffusion)。 5.5.3 网格倾斜交叉引起的计算误差 两股速度相同温 度不同的气流相遇, 设气流的扩散系数为 零,当流速与网格线 夹角倾斜时,数值计 算结果误差很大。 气流
0
0
10/60
1. 设x-y坐标分别与来流平行与垂直
A( P ) 1 对P控制容积,有 采用一阶迎风格式,
aE De Fe ,0 U 0, 00 aW DW FW ,0U 0, 0 F aN Dn Fn ,0 V 0, 0 aS Ds Fs ,0 V 0, 0
d ( u ) d d ( ) S , S 非常数,分布给定 dx dx dx
x 0, 0 ; x L, L
15/60
对于这种具有非常 数源项的情形,无论五 种三点格式中的哪一 种,在一定条件下都会 产生较大的计算误差。 以混合格式为例,网格Peclet数小于2时,数值 解与分析解符合甚好,但大于2 后则偏差甚大。
e
E P
2
1 Cur 8
中心差分插值
(方位定义) 界面曲 率修正
界面的 实际值
26/60
如何确定曲率修正CUR?需要满足两个条件: (1)能自动反应型线凹向对CD插值的正确修正:相邻 三点间二阶中心差分的结构可以反应型线的凹向 型线上凹 (W 2 P 型线下凹 (W 2 P
8/60
由离散方程:
in1 in
t
得:
u
in in 1
x
纯对流传递

n 1 i
ut ut n (1 )( )i 1 x x
n i
ut 当 1 时 x

n 1 i

n i 1
C 0.8
此时只有对流,没有扩散! 与物理问题相一致! 此值不为1时则有严重假扩散!
7/60
FUS Physically plausible solution
FUS严重的误差
CD振荡的数 值解 2. 一维非稳态对流问题(Noye,1976)
u , 0 x 1, u 0.1, (0, t ) (1, t ) 0 t x
x [0,0.1] 范围内 ,初始分布为三角形,其余为零。
Jin WW, He YL, TaoWQ. How many secondary flows are in Leonard’s vertical slot? Progress in Computational Fluid Dynamics, 2009, 9(3/4):283-291
节点数=316404
2 2 ( ) ( u ) u ( ) u ( u ) u 2 x x x 2 t x x t t t t
相关文档
最新文档