集总参数法( Bi 0 )
-
Upload
chandler-wiley -
Category
Documents
-
view
258 -
download
2
description
Transcript of 集总参数法( Bi 0 )
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
集总参数法( Bi 0 )
0
exp( )V VBi FoV Ve Bi Fo
2
( / ),
( / )V V
h V A aBi Fo
V A
( )V
h V ABi
0.05
0.1
0.033Bi ≈ 0?
平板
圆柱
球
Vc
hA
etttt
00
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
0 0 0
( , ) m
m
t x t
t t
( 3-30 )
2 20
( , )( , , ) ( , , ) ( , , )n
x a x a h x xf f f Fo Bi
1 120
1( , ) ( , )m af f Foh Bi
2 2
1( , ) ( , )
m
x xf fh Bi
图( 3-7 )
图( 3-8 )
22
20
( ) ( )Q h a
f FoBi fQ
图( 3-9 )
第四章
热传导问题的数值解法
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
1 、重点内容: ① 掌握稳态导热问题数值解法的基本思路; ② 利用热平衡法建立节点的离散方程。2 、掌握内容:稳态导热问题的离散方程的建立。 3 、了解内容:了解非稳态导热问题的两种差分格式及其稳定性。
基本要求
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
分析解法与数值解法的异同点: • 相同点:根本目的是相同的,即确定 ① t=f(x , y, z) ;②导热量 。 • 不同点:数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场;分析解法求解的是连续的温度场的分布特征,而不是分散点的数值。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
• 数值解法的实质 对物理问题进行数值解法的基本思路可以概括为:把原来在时间、空间坐标系中连续的物理量的场,如导热物体的温度场等,用有限个离散点上的值的集合来代替,通过求解按一定方法建立起来的关于这些值的代数方程,来获得离散点上被求物理量的值。该方法称为数值解法。 这些离散点上被求物理量值的集合称为该物理量的数值解。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
§4-1 导热问题数值求解的基本思想
§4-2 内节点离散方程的建立方法
§4-3 边界节点离散方程的建立 及代数方程的求解
§4-4 非稳态导热问题的数值解法
以二维稳态导热为例
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
一、问题提出
1 、对于一维稳态导热,可用理论法求解2、若H 不满足,二维导热,如图
3 、二维稳态无内热源的导热微分方程式2 2
2 20
t t
x y
( a)
•此式求理论解困难,边界复杂时则不可能
•用数值法求解 区域离散化 建立离散方程
求离散点的温度值 温度场(分布)
§4-1 导热问题数值求解的基本思想2 2 2
2 2 2( )
t t t t
c x y z c
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
建立控制方程及定解条件 确定节点(区域离散化)
建立节点物理量的代数方程设立温度场的迭代初值
求解代数方程组
是否收敛
解的分析
改进初场
是
否
二、导热问题数值求解基本步骤
0t
y
f3th
f2th
f1th x
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
0t
y
f3th
f2th
f1th x
二维矩形域内稳态无内热源,常物性的导热问题
例题条件
( a )
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
如图( a)所示二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下:
( 1)建立控制方程及定解条件 针对图示的导热问题,它的控制方程(即导热微分方程)为: 2 2
2 20
t t
x y
0t
y
f3th
f2th
f1th x
( a)
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
( b )x
y
x
y
n
m
(m,n)
M
N
基本概念:网格线、节点、步长、控制容积(元体)
二维矩形域内稳态无内热源,常物性的导热问题
( 2)区域离散化(确立节点)
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
用一系列与坐标轴平行的网格线把求解区域划分成若干个子区域,用网格线的交点作为需要确定温度值的空间位置,称为节点 ( 结点 ) ,节点的位置用该节点在两个方向上的标号 m , n 表示。
相邻两节点间的距离称步长,△ x,△ y。
如图 (b) 所示。 x
y
x
y
n
m
(m,n)
M
N
控制容积:节点代表的区域
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
( 3)建立节点物理量的代数方程(离散方程) 节点上物理量的代数方程称离散方程。其过程如下: • 首先划分各节点的类型; • 其次,建立节点离散方程; • 最后,代数方程组的形成。 对节点 (m,n) 的代数方程,当 △ x=△y 时,有:
, 1, 1, , 1 , 1
1( )
4m n m n m n m n m nt t t t t ( b )
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
( 4) 设立迭代初场 代数方程组的求解方法有直接解法与迭代解法,传热问题的有限差分法中主要采用迭代法。采用迭代法求解时,需对被求的温度场预先设定一个解,这个解称为初场,并在求解过程中不断改进。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
( 5) 求解代数方程组
求解时遇到的问题: ① 线性; ② 非线性; ③ 收敛性等。 如图 ( b ),除 m=1 的左边界上各节点的温度已知外,其余 (M-1)N 个节点均需建立离散方程,共有 (M-1)N 个方程,则构成一个封闭的代数方程组。 1 )线性代数方程组:代数方程一经建立,其中各项系数在整个求解过程中不再变化;
x
yx
y
n
m
(m,n)
M
N
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
2 )非线性代数方程组:代数方程一经建立,其中各项系数在整个求解过程中不断更新。
3 )是否收敛判断:是指用迭代法求解代数方程是否收敛,即本次迭代计算所得之解与上一次迭代计算所得之解的偏差是否小于允许值。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
( 6 ) 解的分析
通过求解代数方程,获得物体中的温度分布,根据温度场应进一步计算通过的热流量,热应力及热变形等。因此,对于数值分析计算所得的温度场及其它物理量应作详细分析,以获得定性或定量上的结论。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
§4-2 内节点离散方程的建立方法这是导热问题数值计算的关键一步。
要得出节点的离散方程,首先要了解
该节点是哪种类型。
如图所示共给出了 6 种不同的节点:
( 1 )具有对流边界条件的外角顶;( 2 )具有对流边界条件的平直边界节点;( 3 )具有对流边界条件和对称绝热角顶;( 4 )具有绝热边界条件的平直边界节点;( 5 )具有对流边界条件的内角顶;( 6 )内部节点。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
1、 Taylor(泰勒)级数展开法:
2、 控制容积平衡法(也称为热平衡法) 能量守恒
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
一、 泰勒级数展开法
根据泰勒级数展开式,用节点 (m,n) 的温度 tm,n
来表示节点 (m+1,n) 的温度 tm+1,n
用节点 (m,n) 的温度 tm,n来表示节点 (m- 1,n)的温度 tm-1,n
2 2 3 3 4 4
1, , 2 3 4, , ,
( )2 6 24m n m n
m n m n m n
t x t x t x tt t x a
x x x x
2 2 3 3 4
1, , 2 3 4, , ,
4( )
2 6 24m n m nm n m n m n
t x t x t x tt t x b
x x x x
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
将以上两式相加可得2 4 4
21, 1, , 2 4
,
2 ( )12m n m n m n
m n
t x tt t t x c
x x
2
2
,m n
t
x
将上式改写成 的表达式,有
)(2
22
,1,,1
,2
2
xox
ttt
x
t nmnmnm
nm
)(2
22
1,,1,
,2
2
yoy
ttt
y
t nmnmnm
nm
同样可得: 表示未明确写出的级数余项中的ΔX 的最低阶数为 2
2 2 3 3 4 4
1, , 2 3 4, , ,
( )2 6 24m n m n
m n m n m n
t x t x t x tt t x a
x x x x
2 2 3 3 4
1, , 2 3 4, , ,
4( )
2 6 24m n m nm n m n m n
t x t x t x tt t x b
x x x x
( 4-1a )
( 4-1b )
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
根据导热问题的控制方程 ( 导热微分方程 )2 2
2 20
t t
x y
1, , 1, , 1 , , 1
2 2
2 20m n m n m n m n m n m nt t t t t t
x y
, 1, 1, , 1 , 1
1( )
4m n m n m n m n m nt t t t t 若 △ x=△y 则有
得 ( 4-2 )
P164
( b )
, 1, 1, , 1 , 12 2
2 2
1 1 1[ ( ) ( )]
2 2 ( ) ( )( ) ( )
m n m n m n m n m nt t t t tx y
x y
1, , 1, , , 1 , , 1 ,( ) ( ) ( ) ( ) 0m n m n m n m n m n m n m n m nt t t t t t t t
)(2
22
,1,,1
,2
2
xox
ttt
x
t nmnmnm
nm
)(2
22
1,,1,
,2
2
yoy
ttt
y
t nmnmnm
nm
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
1, , 1, , , 1 , , 1 ,( ) ( ) ( ) ( ) 0m n m n m n m n m n m n m n m nt t t t t t t t
1, ,*1( )m n m ny t t
x
在稳态下,流向任何节点的热量总和必定为 0
1, ,
, 1 , , 1 ,
*1( )
*1( ) *1( )0
m n m n
m n m n m n m n
y t t
xx t t x t t
y y
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
二、 控制容积平衡法 ( 热平衡法 )
基本思想:是傅里叶导热定律和能量守恒定律的体现。对每个元体,可用傅里叶导热定律写出其能量守恒的表达式。如图所示, 从节点 (m-1,n) 通过界面 w 传导到节点 (m,n) 的热流量:
1, ,m n m nw
t ty
x
同理:通过界面 e,n,s 传导给节点( m,n )的热流量也可求得(省略)
(垂直纸面方向取单位长度)
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
0w e n s
其中,规定:导入元体( m,n )的热流量为正;导出元体( m,n )的热流量为负。在未知温度高低的情况下一律以周围温度都高于该节点温度来列方程。
对元体 (m,n). 根据能量守恒定律可知:
1, ,m n m nt ty
x
1, ,m n m nt t
yx
, 1 ,m n m nt t
xy
, 1 ,m n m nt tx
y
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
说明:① 上述分析与推导是在笛卡儿坐标系中进行的; ② 热平衡法概念清晰,过程简捷; ③ 热平衡法与建立微分方程的思路与过程一致,但不同的是前者是有限大小的元体,后者是微元体。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
§ 4-3 边界节点离散方程的建立 及代数方程的求解
对于第一类边界条件的热传导问题,处理比较简单,因为已知边界的温度,可将其以数值的形式加入到内节点的离散方程中,组成封闭的代数方程组,直接求解。
而对于第二类或第三类边界条件的导热问题,所有内节点的离散方程组成的代数方程组是不封闭的,因未知边界温度,因而应对位于该边界上的节点补充相应的代数方程,才能使方程组封闭,以便求解。
为了求解方便,这里我们将第二类边界条件及第三类边界条件合并起来考虑,用 qw表示边界上的热流密度或热流密度表达式。为使结果更具一般性,假设物体具有内热源 ( 不必均匀分布 )。
三类边界条
件的回顾
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
如图所示 边界节点 (m,n) 只能代表半个元体,若边界上有向该元体传递的热流密度为 qw ,该元体具有内热源 ,据能量守恒定律对该元体有:
1. 边界节点离散方程的建立:(1) 平直边界上的节点
1, ,m n m nt ty
x
yx2
,, 1, , 1 , 1
212
4m n w
m n m n m n m n
x Φ xqt t t t
( 4-4a )
( 4-4b )
Φ
, 1 ,
2m n m nt t x
y
, 1 ,
2m n m nt tx
y
, 2m n
xΦ y
0wyq Φ
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
(2) 外部角顶点
2,
, 1, , 1
21
2 2m n w
m n m n m n
x Φ xqt t t
1, ,
2m n m nt t y
x
yx
如图所示,二维墙角计算区域中 (A.B.C.D.E点),该节点外角点仅代表 1/4 个以 为边长的元体。假设边界上有向该元体传递的热流密度为 ,则据能量守恒定律得其热平衡式为:
x y 、wq
( 4-5b )
( 4-5a )
, 1 ,
2m n m nt t x
y
,
4m nx yΦ
02 w
x yq
Φ
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
(3) 内部角点
22
, 1, , 1 , 1 1,
21 3(2 2 )
6 2w
m n m n m n m n m n
x qxt t t t t
1, ,m n m nt ty
x
yx
如图所示内部角点 F代表了 3/4 个元体,在同样的假设条件下有
( 4-6b )
( 4-6a )
, 1 ,m n m nt tx
y
, 1 ,
2m n m nt t x
y
1, ,
2m n m nt t y
x
,
3
4 m n
x yΦ
0
2 w
x yq
Φ
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
讨论关于边界热流密度的三种情况: ( 1)绝热边界
即令上述各式 即可。 0wq
( 2) 值不为零wq
流入元体, 取正,流出元体, 取负使用上述公式
wq wq
)( ,nmfw tthq ( 3)对流边界此时 ,将此表达式代入上述各方程,并将此项中的 与等号前的 合并。对于 的情形有
,m nt ,m nt x y
2
, 1, , 1 , 1 1,
1 3(2 2 )
6 2m n m n m n m n m n
xt t t t t
22
, 1, , 1 , 1 1,
21 3(2 2 )
6 2w
m n m n m n m n m n
x qxt t t t t
( 4-6b )
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
( a)平直边界
( b)外部角点
( c)内部角点
2
,, 1, , 1 , 1
22 2 2 m nm n m n m n m n f
h x x h xt t t t t
2,
, 1, , 1
22 1
2m n
m n m n m n f
h x x h xt t t t
2
,
, 1, , 1 1, , 1
3 22 3 2
2m n
m n m n m n m n m n f
h x x h xt t t t t t
( 4-7 )
( 4-8 )
( 4-9 )
1, , , 1 ,
,
2 2
04 2
m n m n m n m n
m nw
t t t ty x
x y
x yΦ x yq
,( )f m nh y t t
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
2 代数方程的求解方法
2 ) 迭代法:先对要计算的场作出假设(设定初场),在迭代计算中不断予以改进,直到计算前的假定值与计算结果相差小于允许误差为止的方法,称迭代计算收敛。
1) 直接解法:通过有限次运算获得精确解的方法,如:矩阵求解,高斯消元法。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
迭代法目前应用较多的是:
1 )高斯——赛德尔迭代法:每次迭代计算,均是使用节点温度的最新值。 2 )用雅可比迭代法:每次迭代计算,均用上一次迭代计算出的值。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
雅可比迭代法 高斯——赛德尔迭代法
11 1 12 2 1 1
21 1 22 2 2 2
1 1 2 2
n n
n n
n n nn n n
a t a t a t b
a t a t a t b
a t a t a t b
( 1) ( ) ( ) ( )1 1 12 2 13 3 1
11
( 1) ( ) ( ) ( )2 2 21 1 23 3 2
22
( 1) ( ) ( ) ( )1 1 2 2 , 1 1
1
1
1
k k k kn n
k k k kn n
k k k kn n n n n n n
nn
t b a t a t a ta
t b a t a t a ta
t b a t a t a ta
( 1) ( ) ( ) ( )1 1 12 2 13 3 1
11
( 1) ( 1) ( ) ( )2 2 21 1 23 3 2
22
( 1) ( 1) ( 1) ( 1)1 1 2 2 , 1 1
1
1
1
k k k kn n
k k k kn n
k k k kn n n n n n n
nn
t b a t a t a ta
t b a t a t a ta
t b a t a t a ta
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
设有一三元方程组:
11 1 12 2 13 3 1
21 1 22 2 23 3 2
31 1 32 2 33 3 3
a t a t a t b
a t a t a t b
a t a t a t b
其中 ( i=1,2,3 ; j=1,2,3 )及 是已知的系数(均不为零)及常数。
,i ja ib
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
采用高斯——赛德尔迭代法的步骤: ( 1)将三元方程变形为迭代方程:
1 1 12 2 13 311
2 2 21 1 23 322
3 3 31 1 32 233
1( )
1( )
1( )
t b a t a ta
t b a t a ta
t b a t a ta
11 1 12 2 13 3 1
21 1 22 2 23 3 2
31 1 32 2 33 3 3
a t a t a t b
a t a t a t b
a t a t a t b
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
( 2)假设一组解(迭代初场),记为:
并代入迭代方程求得第一 次解 每次计算均用最新值代入。
(0) (0) (0)1 2 3t t t、 、
(1) (1) (1)1 2 3t t t、 、
( 3)以新的初场 重复计算,直到相邻两次迭代值之差小于允许值,则称迭代收敛,计算终止。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
判断迭代是否收敛的准则:
)(max
)()1(
)(
)()1(
)()1(
max
max
max
k
ki
ki
ki
ki
ki
ki
ki
t
tt
t
tt
tt
k及 k+1 表示迭代次数; —第 k次迭代得到的最大值
(k)maxt
当有接近于零的 t 时,第三公式个较好
3 6
10 ~ 10
— 允许的偏差;相对偏差 值一般
取
(4-10)
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
说明:
1 )对于一个代数方程组,若选用的迭代方程式不合适,有可能导致发散,即称迭代过程发散;
2 )对于常物性导热问题,组成的差分方程组,迭代公式的选择应使一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值的代数和,此时,结果一定收敛。
11 1 12 2 13 3 1a t a t a t b 12 13
11
1a a
a
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
12 13 21 23 31 32
11 22 33
1 1 1a a a a a a
a a a
, ,
这一条件数学上称主对角线占优(对角占优);
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
3 )采用热平衡法导出差分方程时,若每一个方程都选用导出该方程中心节点的温度作为迭代变量,则上述条件必满足,迭代一定收敛。
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
02
20)(
222 4,164,163,164,164,15
ttx
hx
y
tty
x
tt
之所以除以 2 是因为1mm 宽的电热带在该元体上只有一半
河海大学常州校区热能与动力工程系—传热学 23/4/20 15:37 杨祥花
作业 P188 4-9