传热学-第四章导热问题的数值解法
合集下载
传热学导热问题的数值解法
导热问题的数值解法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 表示。
四章节导热问题数值解法
O(h2)
(h)
由式(b)和式(d)消去f (x) 得：
f (x)
f (x)
f
(
x
2h) h2
2
f
(x
h)
O(h2
)
(i)
由式(a)和式(b)消去f (x) 得： f (x) f (x h) f (x h) 2 f (x) O(h3) (j) h2
由(e)式～(j)式分别略去 h 、h2 及 h3 以上各项得一阶、二阶
导数向前、向后及中心差分公式为：
、
一阶导数向前差分：
f (x) f (x h) f (x)
h
一阶导数向后差分： f (x) f (x) f (x h) h
一阶导数中心差分：
f (x) f (x h) f (x h) 2h
3 三种方法的特点 (1) 分析法
a 能获得所研究问题的精确解，可以为实验和数值计算提供 比较依据；
b 局限性很大，对复杂的问题无法求解； c 分析解具有普遍性，各种情况的影响清晰可见。
(2) 数值法
在很大程度上弥补了分析法的缺点，适应性 强，特别对于 复杂问题更显其优越性；与实验法相比成本低。
(3) 实验法
f (x)
fi ,
f (x h)
f i 1 ,
f (x h)
fi
……
1
x
函数 f（x）在点 x 的一、二阶导数的有限差分表达式分别为：
一阶导数向前差分：fi '
fi1 h
fi
一阶导数向后差分：fi '
fi fi1 h
一阶导数中心差分：fi '
《传热学》第四章 导热数值解法基础
——以右边界为例 1.第一类边界条件：
边界
2.第二类边界条件：
Байду номын сангаас
Δx＝Δy时简化为：
绝热边界：
3.第三类边界条件：
Δx＝Δy时简化为：
其他情况的节点方程 ——见教材表4－1
外拐角与内拐角节点
对流边界内部拐角节点热平衡：
节点方程式推导实例 ——对流边界外部拐角节点
Δx＝Δy时简化为：
数值导热离散方程组＝内节点离散方程＋边界节点离散方程
二、常用计算软件
1.MATLAB——矩阵计算软件
matlab软件主界面
2.FLUENT——流体流动通用数值计算软件
3. FLUENT AIRPAK ——人工环境系统分析软件，暖通空调专业和传热学领域必备软件
AIRPAK模拟温度场
第四章重点： 1.有限差分方程的建立 2.高斯－赛德尔迭代方法
谢谢观看
《传热学》
第四章 导热数值解法基础
本章研究的目的 ——利用计算机求解难以用 分析解求解的导热问题 基本思想 ——把原来在时间、空间坐 标系中连续的物理量的场， 用有限个离散点的值的集合 来代替，通过求解按一定方 法建立起来的关于这些值的 代数方程，来获得离散点 上被求物理量的值。 研究手段——有限差分法
数值导热离散方程组内节点离散方程边界节点离散方程三节点离散方程组的求解迭代法迭代法的原理离散方程组的求解方法消元法方程过多时计算机内存不足迭代法假定初值根据假定的初值求新值并重复此步骤若干次两次计算值足够接近认为达到真实值简单迭代法每次迭代时使用上次迭代的结果允许误差简单迭代法的缺点由于每次迭代中使用与真实值偏差较大的上次迭代的旧值使运算过程接近真实值的时间增加高斯赛德尔迭代法将本次迭代的最新结果立刻代入本次迭代过程计算其他未知值高斯赛德尔迭代法的优点由于每次迭代中使用与真实值偏差较小的本次迭代的新值使运算过程接近真实值的时间缩短第三节非稳态导热的数值计算一显式差分格式研究对象一维非稳态导热问题一维非稳态导热内节点差分方程
边界
2.第二类边界条件：
Байду номын сангаас
Δx＝Δy时简化为：
绝热边界：
3.第三类边界条件：
Δx＝Δy时简化为：
其他情况的节点方程 ——见教材表4－1
外拐角与内拐角节点
对流边界内部拐角节点热平衡：
节点方程式推导实例 ——对流边界外部拐角节点
Δx＝Δy时简化为：
数值导热离散方程组＝内节点离散方程＋边界节点离散方程
二、常用计算软件
1.MATLAB——矩阵计算软件
matlab软件主界面
2.FLUENT——流体流动通用数值计算软件
3. FLUENT AIRPAK ——人工环境系统分析软件，暖通空调专业和传热学领域必备软件
AIRPAK模拟温度场
第四章重点： 1.有限差分方程的建立 2.高斯－赛德尔迭代方法
谢谢观看
《传热学》
第四章 导热数值解法基础
本章研究的目的 ——利用计算机求解难以用 分析解求解的导热问题 基本思想 ——把原来在时间、空间坐 标系中连续的物理量的场， 用有限个离散点的值的集合 来代替，通过求解按一定方 法建立起来的关于这些值的 代数方程，来获得离散点 上被求物理量的值。 研究手段——有限差分法
数值导热离散方程组内节点离散方程边界节点离散方程三节点离散方程组的求解迭代法迭代法的原理离散方程组的求解方法消元法方程过多时计算机内存不足迭代法假定初值根据假定的初值求新值并重复此步骤若干次两次计算值足够接近认为达到真实值简单迭代法每次迭代时使用上次迭代的结果允许误差简单迭代法的缺点由于每次迭代中使用与真实值偏差较大的上次迭代的旧值使运算过程接近真实值的时间增加高斯赛德尔迭代法将本次迭代的最新结果立刻代入本次迭代过程计算其他未知值高斯赛德尔迭代法的优点由于每次迭代中使用与真实值偏差较小的本次迭代的新值使运算过程接近真实值的时间缩短第三节非稳态导热的数值计算一显式差分格式研究对象一维非稳态导热问题一维非稳态导热内节点差分方程
第四章导热问题数值解法基础_传热学
y o
∆x
(i, j) P (i,j-1) S ∆x
(i+1,j) E
x
QV为单位时间控制体内热源 的发热量； 的发热量；ΔΕ为控制体单
位时间内热能的增加量。 位时间内热能的增加量。
由导热傅立叶定律
QW = QE = QS = QN =
N (i,j+1) (i-1,j) W P (i,j-1) S ∆x ∆x (i, j)
∆x
S
建立离散方程的方法一 建立离散方程的方法一： 泰勒级数展开 方法一： 泰勒级数： 泰勒级数：
ti +1, j ( ∆x ) 2 ( ∆x ) 3 ∂t ∂ 2t ∂ 3t = t i , j + ( ) i , j ∆x + ( 2 ) i , j + ( 3 )i , j + ••• 2! 3! ∂x ∂x ∂x
近代发展（1985年-至今） 近代发展（1985年 至今） Singhal 在“Numerical Heat Transfer “撰文 指出了促使数值传热学应用于实际应解决的问 题 前后处理软件的快速发展 巨型计算机的发展促使了并行算法和紊流直接 数值模拟（DNS)和大涡模拟 LES)的发展 和大涡模拟（ 数值模拟（DNS)和大涡模拟（LES)的发展 大型商业软件投放市场 1993年 PHOENICS对中国的禁运被解除 对中国的禁运被解除， 1993年， PHOENICS对中国的禁运被解除，中国 科技大学火灾实验室首先买进了使用权
走向工业应用阶段（1975-1984年 走向工业应用阶段（1975-1984年） 1979年 国际杂志Numerical 1979年，国际杂志Numerical Heat Transfer 创刊，分为Application 创刊，分为Application 和 Fundamentals 1979年 1979年，大型通用软件 PHOENICS(Parabolic,Hyberbolic,Elliptic Series)问世 Numerical Integration Code Series)问世 1979年 Lconard创建了优于中心差分格式的 1979年，Lconard创建了优于中心差分格式的 QUICK格式 格式( QUICK格式(精度高和稳定性好 ） 1980年 Patanker教授的名著 教授的名著“ 1980年，Patanker教授的名著“Numerical Flow” Heat Transfer and Fluid Flow”出版 随后，许多商用软件如FLUENT,Star FLUENT,Star随后，许多商用软件如FLUENT,Star-CD, CFX 问世
∆x
(i, j) P (i,j-1) S ∆x
(i+1,j) E
x
QV为单位时间控制体内热源 的发热量； 的发热量；ΔΕ为控制体单
位时间内热能的增加量。 位时间内热能的增加量。
由导热傅立叶定律
QW = QE = QS = QN =
N (i,j+1) (i-1,j) W P (i,j-1) S ∆x ∆x (i, j)
∆x
S
建立离散方程的方法一 建立离散方程的方法一： 泰勒级数展开 方法一： 泰勒级数： 泰勒级数：
ti +1, j ( ∆x ) 2 ( ∆x ) 3 ∂t ∂ 2t ∂ 3t = t i , j + ( ) i , j ∆x + ( 2 ) i , j + ( 3 )i , j + ••• 2! 3! ∂x ∂x ∂x
近代发展（1985年-至今） 近代发展（1985年 至今） Singhal 在“Numerical Heat Transfer “撰文 指出了促使数值传热学应用于实际应解决的问 题 前后处理软件的快速发展 巨型计算机的发展促使了并行算法和紊流直接 数值模拟（DNS)和大涡模拟 LES)的发展 和大涡模拟（ 数值模拟（DNS)和大涡模拟（LES)的发展 大型商业软件投放市场 1993年 PHOENICS对中国的禁运被解除 对中国的禁运被解除， 1993年， PHOENICS对中国的禁运被解除，中国 科技大学火灾实验室首先买进了使用权
走向工业应用阶段（1975-1984年 走向工业应用阶段（1975-1984年） 1979年 国际杂志Numerical 1979年，国际杂志Numerical Heat Transfer 创刊，分为Application 创刊，分为Application 和 Fundamentals 1979年 1979年，大型通用软件 PHOENICS(Parabolic,Hyberbolic,Elliptic Series)问世 Numerical Integration Code Series)问世 1979年 Lconard创建了优于中心差分格式的 1979年，Lconard创建了优于中心差分格式的 QUICK格式 格式( QUICK格式(精度高和稳定性好 ） 1980年 Patanker教授的名著 教授的名著“ 1980年，Patanker教授的名著“Numerical Flow” Heat Transfer and Fluid Flow”出版 随后，许多商用软件如FLUENT,Star FLUENT,Star随后，许多商用软件如FLUENT,Star-CD, CFX 问世
《传热学》第4章-导热问题的数值解法
3
4-2. 节点温度差分方程组的求解方法
导热物体所有内部节点和边界节点温度的差分方程都是线性代 数方程。 n个未知节点温度，n个代数方程式：
a11t1 + a12t2 + L + a1 jt j + L + a1ntn = b1
a21t1 + a22t2 + L + a2 jt j + L + t2ntn = b2
空间步长
4
2) 节点温度差分方程的建立
控制 容积
（1）内部节点温度差分方程
对于常物性、无内热源的无限大平壁 的一维非稳态导热问题
热平衡：在k时刻，单位时间内从相邻控制
容积i-1与i+1分别导入的热流量与之和等于该 控制容积热力学能的增加
Φλ′ + Φλ′′ = dU
节点i 的温度对时间的变化率采用向前差分
≤ε
k及k+1表示迭代次数；
t
(k) max
—第k次迭代得到的最大值
当有接近于零的t时，第三个较好
有时还要同时考虑热流密度收敛
4-3. 非稳态导热问题的数值解法
非稳态导热与稳态导热的主要区别：控制方程中多一个非稳 态项；温度随空间和时间变化
∂t ∂τ
=
a(
∂2t ∂x 2
+
∂2t ∂y 2
)
能量平衡关系：网格单元不仅与相邻的网格单元之间有热量的 导入或导出，网格单元本身的热力学能也随时间发生变化
t t 在用第二个方程计算节M点温度
1 2 时，直接将
依a此n1类t1 推+ an2t2 + L + anjt j + L + anntn = bn
传热学—第4章 热传导问题的数值解法
(k ) t max
⎧a11t1 + a12 t2 + a13t3 = b1 ⎪ ⎨a21t1 + a22 t2 + a23t3 = b2 ⎪a t + a t + a t = b 33 3 3 ⎩ 31 1 32 2
假定初场
⎧ (1) ⎪t1 = ⎪ ⎪ Jacobi ⎨t(1) = 2 ⎪ ⎪ (1) ⎪t3 = ⎩
4.1.1 4 1 1 基本思想 把原来在时间、空间坐标系中连续的物理量的场， 用有限个离散点上的值的集合来代替，通过求解按 定方 建 起来 关 值 代数方程 来获 一定方法建立起来的关于这些值的代数方程，来获 得离散点上被求物理量的值。 这些离散点上被求物理量值的集合称为该物理量 的数值解。
4.1.1 基本思想
λ Δy
Δx = Δy 时： tm −1,n
+ tm+1,n + tm,n+1 + tm,n−1 − 4tm,n = 0
tm ,n
1 = ( tm−1,n + tm+1,n + tm,n+1 + tm ,n−1 ) 4
与Taylor级数法相比，热平衡法物理意义明显。
4.3.1 边界节点离散方程的建立
4-2 内部节点离散方程的建立
4.2.1 4 2 1 Taylor级数展开法
4-2 内部节点离散方程的建立 内部节点离散方程的建
∂ 2t ∂x 2
=
m ,n
tm+1 n − 2tm ,n + tm −1 n 1, 1, Δx 2
控制方程
∂ 2t ∂ 2t + =0 ∂x 2 ∂y 2
∂ 2t ∂y 2
⎧a11t1 + a12 t2 + a13t3 = b1 ⎪ ⎨a21t1 + a22 t2 + a23t3 = b2 ⎪a t + a t + a t = b 33 3 3 ⎩ 31 1 32 2
假定初场
⎧ (1) ⎪t1 = ⎪ ⎪ Jacobi ⎨t(1) = 2 ⎪ ⎪ (1) ⎪t3 = ⎩
4.1.1 4 1 1 基本思想 把原来在时间、空间坐标系中连续的物理量的场， 用有限个离散点上的值的集合来代替，通过求解按 定方 建 起来 关 值 代数方程 来获 一定方法建立起来的关于这些值的代数方程，来获 得离散点上被求物理量的值。 这些离散点上被求物理量值的集合称为该物理量 的数值解。
4.1.1 基本思想
λ Δy
Δx = Δy 时： tm −1,n
+ tm+1,n + tm,n+1 + tm,n−1 − 4tm,n = 0
tm ,n
1 = ( tm−1,n + tm+1,n + tm,n+1 + tm ,n−1 ) 4
与Taylor级数法相比，热平衡法物理意义明显。
4.3.1 边界节点离散方程的建立
4-2 内部节点离散方程的建立
4.2.1 4 2 1 Taylor级数展开法
4-2 内部节点离散方程的建立 内部节点离散方程的建
∂ 2t ∂x 2
=
m ,n
tm+1 n − 2tm ,n + tm −1 n 1, 1, Δx 2
控制方程
∂ 2t ∂ 2t + =0 ∂x 2 ∂y 2
∂ 2t ∂y 2
传热学课件第四章 导热问题数值解法基础
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
传热学 第4章-导热问题的数值解法
第四章 导热问题的数值解法1、重点内容： ① 掌握导热问题数值解法的基本思路；② 利用热平衡法和泰勒级数展开法建立节点的离散方程。
2、掌握内容：数值解法的实质。
3、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。
由前述3可知，求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。
但是，对于工程中几何形状及定解条件比较复杂的导热问题，从数学上目前无法得出其分析解。
随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，并得到广泛应用，并形成为传热学的一个分支——计算传热学（数值传热学），这些数值解法主要有以下几种：（1） 有限差分法 （2）有限元方法 （3）边界元方法数值解法能解决的问题原则上是一切导热问题，特别是分析解方法无法解决的问题。
如：几何形状、边界条件复杂、物性不均、多维导热问题。
分析解法与数值解法的异同点：1、 相同点：根本目的是相同的，即确定① t=f(x ，y ，z)；② ),,,(τz y x g Q =。
2、 不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。
§4—1 数值求解的基本思路及稳态导热内节点离散方程的建立一、 解法的基本概念1、 实质对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。
该方法称为数值解法。
这些离散点上被求物理量值的集合称为该物理量的数值解。
2、基本思路：数值解法的求解过程可用框图4-1表示。
由此可见：1）物理模型简化成数学模型是基础； 2）建立节点离散方程是关键；3）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。
大学传热学第四章 导热问题的数值解法1
对解进行分析
获得物体中的温度分布常常不是工程问题的最终目的，所 得出的温度场可能进一步用于计算热流量或计算设备、零 部件的热应力及热变形等。对于数值计算所获得的温度场 及所需要的一些其他物理量应作仔细的分析，以获得定性 或定量上的一些新的结论。
建立节点方程的泰勒级数展开法
• 函数的泰勒级数展开式为
二维稳态导热内部节点方程式的建立
ydy m,n+1
m-1,n
m,n
x
y
x m,n-1
m+1,n
xdx
y
• 从左面进入微元体的热量
t t
y m1,n
m ,n
x
x
• 从右面进入微元体的热量
t t
y m1,n
m ,n
xdx
x
• 从下面进入微元体的热量
t t
x m,n1
m ,n
y
y
• 从上面进入微元体的热量
t t
x m,n1
m ,n
ydy
y
二维稳态导热内节点方程
• 当物体内没有内热源时，根据能量守恒定律，从各个方向 进入微元体的热量之和为零
0
x
xdx
y
ydy
• 将各热量计算表达式带入，整理后得到
t 1 t t t t
m ,n
4 m1,n
m1,n
m ,n1
y
y
从上面进入微元体的热量
t t
m,n1 m,n x / 2
ydy
y
二维稳态导热平直边界上节点方程
• 当物体内没有内热源时，根据能量守恒定律，从各个方向 进入微元体的热量之和为零
0
m,n
m,n
传热学第四章导热问题的数值解法
三.灰体表面间的辐射换热
由于灰体表面存在着多次的反射和吸收，计算起来比 黑体复杂得多，为了计算方便我们引进了两个概念： 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
本节讨论的是被透热介质隔开的两固体表面间的辐射换 热。透热介质是指不参与热辐射的介质，最常见的是空 气。
传热学第四章
k及k+1表示迭代次数；
t
(k) max
—第k次迭代得到的最大值
华北电力大学
刘彦丰
传热学 Heat Transfer
4-3 一维非稳态导热问题的数值求解
在非稳态导热问题中，不但需要对空间区域进 行离散，还需要对时间变量进行离散，接下来以一 个一维非稳态导热问题为例，重点介绍对非稳态项 的离散方法，以及不同离散方法对计算带来的影响 等。
n
y
m
x
华北电力大学
刘彦丰
传热学 Heat Transfer 2.建立节点物理量的代数方程
每一个节点都与它周围相邻的节点存在一定的 关系，通过相应的物理定律，可建立它们之间的关 系式（属于代数方程式），此关系式又称作节点的 离散方程。
(m,n+1)
(m,n)
(m-1,n)
(m+1,n)
华北电力大学
(m,n-1)
h,tf
0 q=0 H x
华北电力大学
刘彦丰
传热学 Heat Transfer
二、数学描述
华北电力大学
∂2t ∂x 2
+
∂2t ∂y 2
+
Φ& λ
=
0
x=0 x=H y=0
−
∂t ∂x
=
0
−
λ
∂t ∂x
=
h(t
−tf
)
−
∂t ∂y
=
0
y =W
−
λ
∂t ∂y
=
h(t
−tf
)
刘彦丰
传热学 Heat Transfer
tm,n+1 −tm,n ∆y
(m,n+1)
传热学-第四章-热传导问题的数值解法
23
判断迭代是否收敛的准则：
迭代次数，表示第k次迭代
Monday, March 30, 2020
表示第k次迭代所得计算域内的最大值 当有温度t接近于零的时，选此准则较好
24
例题：
Monday, March 30, 2020
25
Monday, March 30, 20day, March 30, 2020
27
1. 一维非稳态导热的数值求解： 第三类边界条件下，常物性、无内热源无 限大平壁的一维非稳态导热问题为例。
1) 求解域的离散
2) 节点温度差分方程的建立
运用热平衡法可以建立非稳态导热物体内部节点和 边界节点温度差分方程。
Monday, March 30, 2020
29
➢ 两点结论：
(a) 任意一个内部节点n在(i+1)时刻的温度都可以由该节点及 其相邻节点(n-1) 、(n+1)在i 时刻的温度由上式直接求出，不必联 立求解方程组，这是显式差分格式的优点。这样就可以从初始温 度出发依次求出各时刻的节点温度；
(b) 必须满足显式差分格式的稳定性条件，即
物理意义:
15
§4-3 边界节点离散方程的建立及代数方程的求解
第一类边界条件：已知全部边界的温度，作为已知值加入到内节点的离散方程中， 组成封闭的代数方程组，直接求解。
n=N
封闭
(m,n+1)
第二类边界条件或第三类边界 条件：部分边界温度未知。
不封闭
w (m-1,n)
n e
(m,n) s
(m,n-1)
(m+1,n)
y
n=1
m=1
m
x
m=M
Monday, March 30, 2020
传热学-第4章-热传导问题的数值解珐
若步长∆x=∆y，有： ， 若步长
t m ,n = 1 ( 2 t m −1 , n + t m , n + 1 + t m , n −1 + 4 ∆2 x Φ m , n
λ
+
2 ∆ xq w
λ
)
2. 外部角点 控制容积的热平衡为： 控制容积的热平衡为：
∆y tm−1,n − tm,n ∆x tm,n−1 − tm,n ∆x∆y ∆x + ∆y λ +λ + Φ m, n + qw = 0 ∆x 2 2 ∆y 4 2
4. 边界热流密度的三种情况
q （1）绝热边界： w = 0 ）绝热边界：
（2） qw 值不为零：代入给定的 qw 值。 ） 值不为零： （3）对流边界：qw = h(t f ）对流边界： 平直边界节点： 平直边界节点：
2( h∆x
− t m n = 2 t m − 1 , n + t m , n + 1 + t m , n −1 +
第一类边界条件 — 边界温度已知 m-1,n 第二类边界条件 需建立边界节点温度 ∆y 第三类边界条件 的差分方程 n 1. 位于平直边界上的节点
λ∆y
tm−1,n − tm,n ∆x +λ
m m,n+1
qw
m,n m,n-1
∆x
∆x tm,n+1 − tm,n ∆x tm,n−1 − tm,n ∆x∆y +λ + Φm,n + ∆yqw = 0 2 ∆y 2 ∆y 2
若步长∆x=∆y，有： ， 若步长
t m ,n = 1 ( t m −1 , n + t m , n −1 + 2
传热学第四章-导热问题的数值解法-2
1. 节点离散方程的建立：
(1)内部节点
相邻节点导入控制单元体的热流量= 单元体内能量增量
i-1
i
i+1
A ti( k 1 ) ti(k )A ti( k 1 ) ti(k )c x A ti(k 1 ) ti(k )
x
x
整理，得：
x
x
ti(k 1 )[12 ( a x )2]ti(k)(a x)2[ti( k1 )ti( k1 )]
取103 ~ 106
k及k+1表示迭代次数； tm(ka)x —第k次迭代得到的最大值
当有接近于零的t 时，第三个较好
4-3 非稳态导热问题的数值解法
非稳态导热问题与稳态导热问题的区别是，温度分布不仅 与空间坐标有关，还与时间有关。 本节要求掌握一维非稳态导热问题的数值解法，能够写出 内部节点和边界节点的有限差分方程，掌握显式差分方程 的稳定性条件。
2.节点方程组的求解： 步骤：
1）选择坐标和时间的步长，按选定的坐标步长划分节点网 格，并将节点按位置编号。
2）按节点的情况（位置和具体边界条件）写出各节点的差
分方程，并检查是否符合稳定性 条 件。
3）从初始条件出发，逐点计算 时刻各节点的温度，然后
再逐点计算 2,3,...... 时刻各节点的温度，直到指定
例 如 t03 t1 3， 但 t0 4<t1 4。
i0 1 2 3 4 5 6 7
t
n
0
100 100 100 100 60
148 -109.6 550
1
100 100 100 80
104 19.2 220.2 -328.9
2
100 100 80
84
传热学-学习课件-4-1 导热问题数值求解基本思想
二、主要的数值解法
有限差分法（FDM），主要介绍方法 有限元法（FEM） 边界元法（BEM）
传热学 Heat Transfer
三、数值解法的基本步骤
传热学 Heat Transfer
以一个二维稳态导热问题为例，介绍数值求解导热问题的具体 过程，重点是节点离散方程的建立和代数方程组的迭代求解。
（1）物理问题 二维矩形域稳态无内热源，常物性的导热问题 y
传热学 Heat Transfer
4-1 导热问题数值求解的基本思想
一、数值解法的本质
数值解法是用物理问题所 涉及的空间和时间区域内有 限个离散点（称为节点）的 物理量近似值来代替物体内 实际连续的物理量分布，将 连续物理量分布函数的求解 问题转化为各节点物理量值 的求解问题。
传热学 Heat Transfer
传热学 Heat Transfer
传热学 Heat Transfer
主讲老师：王舫 适用专业：能源与动力工程专业
传热学 Heat Transfer
第四章 热传导问题的数值解法
§4-0 引言
1 求解导热问题的三种基本方法： (1) 理论分析 (2)实验 (3)数值计算
2 三种方法的特点 3 三种方法的基本求解过程
传热学heattransfer?本章教学内容41导热问题数值求解基本思想42内节点离散方程的建立43边界节点离散方程的建立及代数方程的求解44非稳态导热问题的数值解法传热学heattransfer41导热问题数值求解的基本思想一数值解法的本质数值解法是用物理问题所涉及的空间和时间区域内有限个离散点称为节点的物理量近似值来代替物体内实际连续的物理量分布将连续物理量分布函数的求解问题转化为各节点物理量值的求解问题
传热学 Heat Transfer
有限差分法（FDM），主要介绍方法 有限元法（FEM） 边界元法（BEM）
传热学 Heat Transfer
三、数值解法的基本步骤
传热学 Heat Transfer
以一个二维稳态导热问题为例，介绍数值求解导热问题的具体 过程，重点是节点离散方程的建立和代数方程组的迭代求解。
（1）物理问题 二维矩形域稳态无内热源，常物性的导热问题 y
传热学 Heat Transfer
4-1 导热问题数值求解的基本思想
一、数值解法的本质
数值解法是用物理问题所 涉及的空间和时间区域内有 限个离散点（称为节点）的 物理量近似值来代替物体内 实际连续的物理量分布，将 连续物理量分布函数的求解 问题转化为各节点物理量值 的求解问题。
传热学 Heat Transfer
传热学 Heat Transfer
传热学 Heat Transfer
主讲老师：王舫 适用专业：能源与动力工程专业
传热学 Heat Transfer
第四章 热传导问题的数值解法
§4-0 引言
1 求解导热问题的三种基本方法： (1) 理论分析 (2)实验 (3)数值计算
2 三种方法的特点 3 三种方法的基本求解过程
传热学heattransfer?本章教学内容41导热问题数值求解基本思想42内节点离散方程的建立43边界节点离散方程的建立及代数方程的求解44非稳态导热问题的数值解法传热学heattransfer41导热问题数值求解的基本思想一数值解法的本质数值解法是用物理问题所涉及的空间和时间区域内有限个离散点称为节点的物理量近似值来代替物体内实际连续的物理量分布将连续物理量分布函数的求解问题转化为各节点物理量值的求解问题
传热学 Heat Transfer
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
x
∆x
m
M
第四章 导热问题的数值解法
7
4 建立离散方程的常用方法： 建立离散方程的常用方法： (1) Taylor（泰勒）级数展开法； （泰勒）级数展开法； (2) 多项式拟合法； 多项式拟合法； (3) 控制容积积分法； 控制容积积分法； (4) 控制容积平衡法 也称为热平衡法 控制容积平衡法(也称为热平衡法 也称为热平衡法)
∂ 2t ∂x 2
m ,n
t m +1,n − 2 t m ,n + t m −1,n = + o(∆x 2 ) ∆x 2
同样可得： 同样可得： 截断误差
∂t tm,n+1 − 2tm,n + tm,n−1 = + o(∆y 2 ) ∂y 2 m,n ∆y 2
2
未明确写出的级数余项 的最低阶数为2 中的ΔX的最低阶数为2
第四章 导热问题的数值解法
3
(2) 数值法：在很大程度上弥补了分析法的缺点，适应性 数值法：在很大程度上弥补了分析法的缺点， 强，特别对于复杂问题更显其优越性；与实 特别对于复杂问题更显其优越性； 验法相比成本低 (3) 实验法 是传热学的基本研究方法，a 适应性不好； 实验法: 是传热学的基本研究方法， 适应性不好； b 费用昂贵 数值解法：有限差分法（ 数值解法：有限差分法（finite-difference）、 ）、 有限元法（ 有限元法（finite-element） 、 ） 边界元法（ 边界元法（boundary- element）、 ）、 分子动力学模拟（ 分子动力学模拟（MD） ）
第四章 导热问题的数值解法 2
(3) 实验法 就是在传热学基本理论的指导下，采用对所 就是在传热学基本理论的指导下， 研究对象的传热过程所求量的方法
3 三种方法的特点 (1) 分析法 a 能获得所研究问题的精确解，可以为实验和数值计算 能获得所研究问题的精确解， 提供比较依据； 提供比较依据； b 局限性很大，对复杂的问题无法求解； 局限性很大，对复杂的问题无法求解； c 分析解具有普遍性，各种情况的影响清晰可见 分析解具有普遍性，
第四章
导热问题的数值解法
§4-0 引言
1 求解导热问题的三种基本方法：(1) 理论分析法；(2) 数 求解导热问题的三种基本方法： 理论分析法； 值计算 法；(3) 实验法 2 三种方法的基本求解过程 (1) 所谓理论分析方法，就是在理论分析的基础上，直接 所谓理论分析方法，就是在理论分析的基础上， 对微分方程在给定的定解条件下进行积分，这样获得的解 对微分方程在给定的定解条件下进行积分， 称之为分析解，或叫理论解； 称之为分析解，或叫理论解； (2) 数值计算法，把原来在时间和空间连续的物理量的场， 数值计算法，把原来在时间和空间连续的物理量的场， 用有限个离散点上的值的集合来代替， 用有限个离散点上的值的集合来代替，通过求解按一定方 法建立起来的关于这些值的代数方程， 法建立起来的关于这些值的代数方程，从而获得离散点上 被求物理量的值；并称之为数值解； 被求物理量的值；并称之为数值解；
第四章 导热问Leabharlann 的数值解法4§4-1 导热问题数值求解的基本思想 及内部节点离散方程的建立
1 物 理 问 题 的 数 值 求 解 过 程
建立控制方程及定解条件 确定节点（区域离散化）
设立温度场的迭代初值
建立节点物理量的代数方程
求解代数方程
改进初场
否 是否收敛 是 解的分析
第四章 导热问题的数值解法 5
2 例题条件
y
h3 t f
二维矩形域内 稳态无内热源， 稳态无内热源， 常物性的导热 问题
t0
h2 t f
h1t f
x
第四章 导热问题的数值解法 6
3 基本概念：控制容积、网格线、节点、界面线、步长 基本概念：控制容积、网格线、节点、界面线、
N
(m,n)
n
∆y
y
二维矩形 域内稳态 无内热源， 无内热源， 常物性的 导热问题
m ,n
m ,n
(i-1,j)的 用节点(i,j)的温度ti,j来表示节点(i-1,j)的
温度ti-1,j
t m −1,n = t m ,n
∂t ∂ 2t − ∆x + 2 ∂x m ,n ∂x
m ,n
∆ x 2 ∂ 3t − 3 2! ∂x
m ,n
∆x 3 +L 3!
9
第四章 导热问题的数值解法
若取上面式右边的前三项，并将式①和式③ 若取上面式右边的前三项，并将式①和式③相加 移项整理即得二阶导数的中心差分： 移项整理即得二阶导数的中心差分：
第四章 导热问题的数值解法
10
第四章 导热问题的数值解法
8
(1) 泰勒级数展开法 根据泰勒级数展开式，用节点( 根据泰勒级数展开式，用节点(i,j)的温度ti,j 来表示节点( 来表示节点(i+1,j)而温度ti+1,j
t m +1,n = t m ,n ∂t ∂ 2t + ∆x + 2 ∂ x m ,n ∂x ∆ x 2 ∂ 3t + 3 2! ∂ x ∆x 3 +L 3!