传热学导热数值计算共44页
传热学导热问题的数值解法
导热问题的数值解法1 、重点内容：① 掌握导热问题数值解法的基本思路；② 利用热平衡法和泰勒级数展开法建立节点的离散方程。
2 、掌握内容：数值解法的实质。
3 、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。
由前述3 可知，求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。
但是，对于工程中几何形状及定解条件比较复杂的导热问题，从数学上目前无法得出其分析解。
随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，并得到广泛应用，并形成为传热学的一个分支——计算传热学（数值传热学），这些数值解法主要有以下几种：（１）有限差分法（ 2 ）有限元方法（ 3 ）边界元方法数值解法能解决的问题原则上是一切导热问题，特别是分析解方法无法解决的问题。
如：几何形状、边界条件复杂、物性不均、多维导热问题。
分析解法与数值解法的异同点：相同点：根本目的是相同的，即确定① t=f(x ，y ，z) ；②。
不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。
§4-1 导热问题数值求解的基本思想及内节点离散方程的建立实质对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。
该方法称为数值解法。
这些离散点上被求物理量值的集合称为该物理量的数值解。
2 、基本思路：数值解法的求解过程可用框图4-1 表示。
由此可见：1 ）物理模型简化成数学模型是基础；2 ）建立节点离散方程是关键；3 ）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。
一数值求解的步骤如图4-2 （a ），二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下：1 建立控制方程及定解条件控制方程：是指描写物理问题的微分方程针对图示的导热问题，它的控制方程（即导热微分方程）为：（a ）边界条件：x=0 时，x=H 时，当y=0 时，当y=W 时，区域离散化（确立节点）用一系列与坐标轴平行的网格线把求解区域划分成若干个子区域，用网格线的交点作为需要确定温度值的空间位置，称为节点( 结点) ，节点的位置用该节点在两个方向上的标号m ，n 表示。
《传热学》第四章 导热数值解法基础
边界
2.第二类边界条件：
Байду номын сангаас
Δx＝Δy时简化为：
绝热边界：
3.第三类边界条件：
Δx＝Δy时简化为：
其他情况的节点方程 ——见教材表4－1
外拐角与内拐角节点
对流边界内部拐角节点热平衡：
节点方程式推导实例 ——对流边界外部拐角节点
Δx＝Δy时简化为：
数值导热离散方程组＝内节点离散方程＋边界节点离散方程
二、常用计算软件
1.MATLAB——矩阵计算软件
matlab软件主界面
2.FLUENT——流体流动通用数值计算软件
3. FLUENT AIRPAK ——人工环境系统分析软件，暖通空调专业和传热学领域必备软件
AIRPAK模拟温度场
第四章重点： 1.有限差分方程的建立 2.高斯－赛德尔迭代方法
谢谢观看
《传热学》
第四章 导热数值解法基础
本章研究的目的 ——利用计算机求解难以用 分析解求解的导热问题 基本思想 ——把原来在时间、空间坐 标系中连续的物理量的场， 用有限个离散点的值的集合 来代替，通过求解按一定方 法建立起来的关于这些值的 代数方程，来获得离散点 上被求物理量的值。 研究手段——有限差分法
数值导热离散方程组内节点离散方程边界节点离散方程三节点离散方程组的求解迭代法迭代法的原理离散方程组的求解方法消元法方程过多时计算机内存不足迭代法假定初值根据假定的初值求新值并重复此步骤若干次两次计算值足够接近认为达到真实值简单迭代法每次迭代时使用上次迭代的结果允许误差简单迭代法的缺点由于每次迭代中使用与真实值偏差较大的上次迭代的旧值使运算过程接近真实值的时间增加高斯赛德尔迭代法将本次迭代的最新结果立刻代入本次迭代过程计算其他未知值高斯赛德尔迭代法的优点由于每次迭代中使用与真实值偏差较小的本次迭代的新值使运算过程接近真实值的时间缩短第三节非稳态导热的数值计算一显式差分格式研究对象一维非稳态导热问题一维非稳态导热内节点差分方程
传热系数与导热系数换算公式
传热系数与导热系数换算公式
传热系数与导热系数之间存在换算关系，具体如下：
热传导率 = 导热系数 / (物质的密度× 物质的比热容)
根据这个公式，我们可以将导热系数和传热系数进行相互换算。
例如，假设某物质的导热系数为W/(m·K)，密度为1000 kg/m³，比热容为1000
J/(kg·K)，我们可以先计算出该物质的热传导率：
热传导率= / (1000 × 1000) = 5 × 10^-7 m²·K/W
然后，通过热传导率可以计算出该物质的传热系数：
传热系数 = 1 / 热传导率= 1 / (5 × 10^-7) = 2 × 10^6 W/(m²·K)
通过以上计算，我们得知了该物质的传热系数为2 × 10^6 W/(m²·K)。
以上内容仅供参考，建议查阅传热学或物理学书籍获取更全面和准确的信息。
传热学数值计算
Fe aE De 2
Thermal
Fw Fe aP Dw De aW aE Fe Fw 2 2
2、对方程的几点说明 由于连续性，Fe=Fw， aP aW aE（只是在流 场满足连续性条件时才具有这一性质）； 方程 aP P aE E aW W 隐含着分段线性分布的含
讨论只有对流项和扩散项存在时的一维稳态问题，控 制方程为：
u j ( ) S x j x j x j
d d d u ( ) dx dx dx
d 连续方程： u 0 dx
u const
任务：导出相应方程的离散化形式
义，也是熟知的中心差分格式（用左右节点值表示
界面上的值以及界面上的导数值）； 方程必须遵守四项基本法则，否则会产生灾难性的 结果。
2018-11-24
太 原 理 工 大 学
7 /70
Thermal
例如：设 De Dw 1， Fe Fw 4 若E、W给定，即可由离散方程求得P 。
即， F 2D时，有可能使 aE 或aW 为负 产生不切实际的结果
2018-11-24
太 原 理 工 大 学
8 /70
Thermal
这就是中心差分格式求解对流换热问题时仅限于低
Fw Fe Re(低的F/D)的原因 . aP Dw De aW aE Fe Fw 2 2
原通用方程可改写为
u j ( )S x j x j x j
对于已知的ρ、uj、Γ及S（常量）的分布，任何解及
+c 将同时满足方程，故系数和的法则仍然适用。
2018-11-24
太 原 理 工 大 学
传热学-第4章-导热数值解法基础
环境与能源工程学院 ( SEEE )
二阶导数的中心差分表达式 [1]-[2]相加
在表示温度对时间的一阶导数时只采用向前或向后差分 表达式——温度对时间的中心差分表达式求解非稳态导 表达式 中心差分 热问题将导致数值解的不稳定，参考：《工程传热学》 导致数值解的不稳定 以常物性、无内热源、二维稳态导热为例 P(I,j)
2ΔxqW 1 = (2ti −1, j + ti , j +1 + ti , j −1 + ) λ 4
第三类边界条件
i,j+1 i1,j i,j
Beijing University of Civil Engineering and Architecture ( BUCEA )
传热学 ( Heat transfer )
环境与能源工程学院：School of Environment and Energy Engineering -SEEE
第四章 导热数值解法基础 4-1 建立离散方程的方法 4-2 稳态导热的数值计算 4-3 非稳态导热的数值计算
ti +1, j ⎛ ∂ 2t ⎞ Δx 2 ⎛ ∂ 3t ⎞ Δx 3 ⎛ ∂t ⎞ = ti , j + ⎜ ⎟ Δx + ⎜ 2 ⎟ ⎜ ∂x ⎟ 2! + ⎜ ∂x 3 ⎟ 3! + LL ⎜ ⎟ ⎝ ∂x ⎠i , j ⎝ ⎠i , j ⎝ ⎠i , j
节点（i, j）一阶导数的向后差分表达式（一阶截差公式）；
P(i.j)
q ∂t ∂ 2t ∂ 2t = a( 2 + 2 ) + v ρc ∂τ ∂x ∂y
二维导热问题；网格；沿x、 y方向的间距为Δx、Δy； 网格单元 （Nodal point） 节点：网格线的交点；p(i,j) 节点
传热学 数值计算
数值计算
4、一无限大平壁厚度为 0.3m，其导热系数为λ =36.4 W/ (m·K)。平壁两侧表面均给定 为第三类边界条件，即 h 1 =60W/ (m 2 ·K)，t f 1 =25℃；h 2 =300W/ (m 2 ·K)，t f 2 =215℃。 当平壁中具有均匀的内热源 q v =2 10 W / m 时，试计算沿平壁厚度的稳态温度分布。
k
0
k k 1
11
k 2Fo (t 10 t f Bi ) (1 2Fo 2Bi Fo) t11 t11
0
四、
计算过程
⑴ 设定初值：
t (1~11) 35 ℃;
36W / m k ;
Bi h / ;
根据不同的 Fo 计算Δ τ ：
qv x 2
)/2
h2 x
h2 x TRB
qv x 2
) /(1
)
|T[i]-t[i]|<=EPS
NO
IT=IT+1
YES
打印“t[i]” ， “IT”
YES
IT>K
NO
停止
⑷ 程序与计算结果 #include"iostream.h" #include"iomanip.h" #include"math.h" #include"stdio.h" #define N 16 void main() { int M,i,IT,flag; //定义节点个数
计算结果：
各节点温度： 节点 1 2 411.24 10 422.65 3 420.33 11 414.23 4 427.22 12 403.62 5
传热学课件第四章 导热问题数值解法基础
i , j
t x
t i 1 , j t i , j x
0 x
2.一阶导级的向后差分表达式：舍去<2>式△x2后各项，则有：
i , j
t x
t i , j t i 1 , j x
0 x
第一节 建立离散方程的方法
二、泰勒级数展开法（有限差分法）
k 2 k 1
对 流 h t f t1 A
k k
显式
△x
C.内能增量△u：
u c
x 2
A t1
k
k 1
t1 /
k
△x/2
k hx
据热平衡A+B=C并整理得：
k f
t 2 t1
k
t
t1
k
1 2
c
x
2
t1
k 1
LP
△y
t i 1 , j t i , j x
t i , j 1 t i , j y
y 2
x 2
1
BP
1
x 2
y 2
EP h t f t i , j
△x
1
FP h t f t i , j
t x
t
2
2
x i , j 2!
2
t x
3
x i , j 3!
3
3.一阶导级的中心差分表达式：<1>-<2>式且忽略后项，则有：
i , j
t x
传热学第四章导热问题的数值解法
三.灰体表面间的辐射换热
由于灰体表面存在着多次的反射和吸收，计算起来比 黑体复杂得多，为了计算方便我们引进了两个概念： 1.投入辐射G和有效辐射J G—单位时间内投射到表面单位面积上的总辐射能。 J—单位时间内离开该表面单位面积的总辐射能。
G(投入辐射) (1-α)G(反射辐射) αG (吸收辐射) J(有效辐射) E= ε Eb (本身辐射)
定义：表面1发出的辐射能直接落到表面2上的百分数， 称表面1对表面2的角系数，记为X1,2，同理有X2,1。 角系数是一个无量纲能量百分比。引入角系数是为了说 明两个表面之间的辐射换热量与它们之间的相对位置有 很大关系。 角系数为几何因子，其值取决于物体的几何特性（形状、 尺寸及物体的相对位置）而与物体的种类和温度无关。
i =1 n
表面为凸面或平面时，有何性质？
3.可加性
从表面1上发出落到表面2上的总能量，等于落到表 面2各部分的能量之和。于是有：
X 1， 2 = X 1， 2 a + + X 1， 2 n = ∑ X 1， 2i
i =a n
设表面由a、b两部分组成，写出其可加性表 达式。
三.角系数的计算
或写成：
Eb − J q = ， 1− ε
ε
1− ε 称为表面辐射热阻 εA
2.灰体表面间的辐射换热
考虑两个等温的漫灰表面组成的二维封闭系统，表 面1、2间的Φ为： φ1, 2 = A1 J 1 X 1, 2 − A2 J 2 X 2,1
J 1 A1 = A1 Eb1 − (1 / ε 1 − 1)φ1, 2 J 2 A2 = A2 Eb 2 − (1 / ε 2 − 1)φ 2,1
本节讨论的是被透热介质隔开的两固体表面间的辐射换 热。透热介质是指不参与热辐射的介质，最常见的是空 气。